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Abstract 


We make an updated review and a systematic and comprehensive analysis of the decays of Higgs 
bosons in the Standard Model (SM) and its three well-defined prototype extensions such as the 
complex singlet extension of the SM (cxSM), the four types of two Higgs-doublet models (2HDMs) 
without tree-level Higgs-mediated flavor-changing neutral current (FCNC) and the minimal super- 
symmetric extension of the SM (MSSM). We summarize the theoretical predictions for the decay 
widths of the SM Higgs boson and those of Higgs bosons appearing in its extensions taking account 
of all possible decay modes. We incorporate them to study and analyze decay patterns of CP-even, 
CP-odd, and CP-mixed neutral Higgs bosons and charged ones. We put special focus on the 
properties of a neutral Higgs boson with mass about 125 GeV discovered at the LHC and present 
constraints obtained from precision analysis of it. This review is intended to be self-contained and 
consolidated by coherently integrating relevant physics information for studying decays of Higgs 
bosons in the SM and beyond. 
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1 Introduction 


Since the discovery of a resonance with a mass of approximately 125 GeV at the Large Hadron Collider 
(LHC) in 2012 [I| P], the substantial subsequent studies of its properties have been carried out with 
the data set collected during the LHC Run 1 period from 2009 to 2012 and the LHC Run 2 period 
from 2015 to 2018. They have firmly confirmed the compatibility of the resonance with the spin-zero 
and parity-even SM Higgs boson which appears in the spontaneously broken gauge theory where the 
electroweak interactions are governed by the SU(2);, x U(1),- gauge symmetry group [3} 4 5]. |] At the 
present time we are on a watershed peak for exploring a new territory of particle physics through the 
Higgs landscape. 


The total and differential rate measurements of all the possible production and decay channels of the 
resonance state so far are consistent with those predicted in the SM within experimental and theoretical 
uncertainties [7] [8]. The mass of the Higgs boson has been measured at the per-mille precision level, 
mainly through the high-resolution decay modes with four-lepton and di-photon final states [9] [10]. 
Furthermore, the couplings of the Higgs boson to the gauge bosons and the charged fermions of the 
third generation [11] [2 {13} [14] (15) and, recently, to the muons [16} [7] were established independently 
and unambiguously. Based on the observational facts, we call the discovered resonance particle as the 
SM-like Higgs boson wherever appropriate in the following. 


Nevertheless, the couplings of the SM-like Higgs boson to the electrons and lighter quarks of the first 
and second generations and its cubic and quartic self-couplings defining the profile of the Higgs poten- 
tial are yet to be established and measured independently. Furthermore, more comp Higgs sectors 
associated with additional states have not been ruled out. t | 


structure or not. 
Conceptually, the SM with the Higgs boson could be weakly interacting well above the weak scale 


of v = 246 GeV without violating unitarity and so with no need for new physics. However the Higgs 
boson mass is influenced subtly by the presence of heavy particles and it receives quantum corrections 
destabilizing the weak scale and requiring a delicate fine-tuning of apparently unrelated parameters. 
This so-called naturalness or hierarchy problem has been the key argument for expecting 
new physics to be revealed at the TeV scale. To mention just a few, new theoretical frameworks based 
on a fermion-boson symmetry called supersymmetry [21] 22 [23], a collective symmetry between the 
SM particles and heavier partners as in Little Higgs theories [24] or an effective reduction of 
the Planck scale to the TeV scale as in extra-dimension models have been proposed and 
intensively investigated. 


‘Note that the gauge and fermion sectors of the SM have been already well probed with great precision both theoret- 
ically and experimentally as can be checked with the particle physics reference book Ref. [6]. 


After the successful completion of Run 1 and Run 2, the LHC is presently in the second long shut 
down period while undergoing important upgrades for its high luminosity phase. Much larger data 
sets are to be collected during the Run 3 period and, ultimately, during the operation period of the 
high-luminosity LHC (HL-LHC) and they will enable us to explore new physics beyond the SM (BSM) 
by performing more challenging as well as more precise measurements. In light of such promising 
experimental prospects at the LHC and at other future high-energy and high-precision experiments 
[38] [39] it is quite timely and worthwhile to perform a systematic and comprehensive review and 
analysis of the decays of Higgs bosons including the theoretical calculations known up to now, not only 
in the SM but also in various BSM scenarios with unique features in their extended Higgs sectors. 


Certainly it is formidable to review all theoretical and experimental aspects of Higgs sectors in all 
BSM models proposed so far in a single report with limited space. Unavoidably, we restrict ourselves 
in this review to the SM and the three well-defined prototype BSM models with extended Higgs sectors 
possessing their own characteristic features and broad implications. Specifically, in addition to the 
SM, we consider the following representative examples: (i) the complex singlet extension of the SM 
(cxSM) KOL [42] [44] [46] [47], [48], (ii) the four types of two Higgs doublet models (2HDMs) 
[49| with natural flavor conservation at the tree level and the 
so-called p parameter close to unity abiding by the stringent experimental constraint on it, and (iii) the 


minimal supersymmetric extension of the SM (MSSM) [60] [61] (62) {63} [64 (65) (66) [67] 168, [69]. For Higgs 
sectors in BSM models beyond cxSM, 2HDMs and MSSM, see, for example, Refs. [70] [71 [721 [73} [74] [75]. 


Several related previous reviews on Higgs physics in the SM and the MSSM can be found in Refs. [76] 
Kári! [79] [88]. A few tailor-made sophisticated computational packages 
have been developed for the mass spectra and decay widths of neutral and charged Higgs bosons in the 
SM and the MSSM with real parameters and with explicit CP violation (94) [95] [96]. 
This review updates the previous works substantially by including two popular BSM models in addition 
to the MSSM and also by allowing for complex parameters leading to CP-violating phenomena [80} [97]. 
We perform a systematic and comprehensive analysis for the decays of neutral and charged Higgs bosons 
in those three prototype BSM models as well as in the SM. We take account of all possible decay modes 
of these models including those into non-SM particles among which some are invisible and/or 
exotic [100]. We anticipate more complete reviews on the Higgs sectors of many other BSM scenarios 
to come out timely in step with more advanced experimental developments. 


In this review, we try to contain all the relevant information needed to implement the up-to-date 
theoretical calculations of Higgs decays. We aim to make it be of pedagogical and practical use especially 
for incorporating corrections beyond the leading order (LO). We elaborate on how the partial decay 
widths of neutral and charged Higgs bosons are calculated at LO and how we treat QCD and electroweak 
(ELW) corrections in each decay mode. |*| We intend to make this review stand-alone, self-contained, 
and consolidated by integrating relevant physics information for Higgs decays coherently. Incidentally, 
we try to make it be as model-independent and analytic as possible in order for our approach to be 
useful and easily applicable in studying Higgs decays even in the BSM models not explicitly mentioned 
in this review. 


This review is organized as follows. Section 2 is devoted to reviewing the Higgs sectors of the SM 
and three extended scenarios - cxSM, 2HDMs and MSSM - without imposing any constraints on the 
model parameters. We work out the analytic structure of the Higgs potential and mixing. Also worked 
out are the Higgs interactions with gauge bosons, the SM fermions, and new scalars and fermions as 
well as the Higgs-boson self interactions. We review and update the decays of neutral Higgs bosons in 


?The ELW corrections considered in this review are mostly SM ones since the BSM ELW corrections, compared to 
universal QCD corrections, are still subleading, complicated, and strongly dependent on specific BSM models under 
consideration. 


Section [3] and those of charged Higgs bosons in Section [4] We provide explicit analytical expressions of 
the individual partial decay widths as precisely as possible by including the state-of-the-art theoretical 
calculations. In Section |5| we present the constraints on the couplings of the SM-like Higgs boson 
weighing about 125 GeV obtained from global fits to the LHC precision Higgs data. Conclusions are 
made in Section [6| To make this review self-contained, various supplemental materials are provided in 
six appendices. Appendix{A]includes a summary of the SM parameters used for the numerical estimates 
of the Higgs decay widths and a description of the running of the strong coupling constant and heavy 
quark masses. In Appendix |B| the supersymmetric contributions to the loop-induced couplings of the 
Higgs boson to two gluons, two photons and Zy are presented and Appendix |C] is devoted to the 
presentation of the QCD corrections to the partial width of the Higgs-boson decay to two photons. In 
Appendix [D] we present expressions for the most general 2HDM potential parameters in terms of the 
masses of charged and neutral Higgs bosons and the elements of the orthogonal matrix describing the 
mixing among neutral Higgs bosons and, in Appendix [E] we apply them for deriving cubic Higgs-boson 
self-couplings. Finally, Appendix |F|is added as a guide to numerical packages for calculating precise 
SM and full BSM-dependent ELW corrections. 


2 Standard Model and Beyond 


In this section, we derive and describe the basic form of Higgs boson masses and mixing as well as 
their interactions in the SM, cxSM, 2HDMs and MSSM. The derived analytical results are utilized 
comprehensively in the sequential sections for the detailed review of the decays of neutral and charged 
Higgs bosons and also for the model-independent precision study of the SM-like neutral Higgs boson 
which has been extensively probed at the LHC since its discovery in 2012. 


We T a eae are eh ae as a pa mo [102] [03], Sarah [104, [105] {106} [107], and 
atic gener: of k the Feynman rules in 


the o Te as PAA T: as in P ET T se ras EE using these codes, we start by presenting 
interaction Lagrangians in order for the readers to work out independently the analytic structure and 
parametric dependence of the partial decay widths of neutral and charged Higgs bosons and to more 


deeply understand the theoretical and phenomenological aspects of Higgs physics in the SM and beyond. 


2.1 Standard Model 


The self-interactions of the SM Higgs boson and its interactions with the massive vector bosons are 
derived from the Higgs Lagrangian: 


Luiss = (DB)! (Dy) — VD), (1) 


where ® denotes a complex SU(2);, doublet Higgs field with hypercharge Y = 1/2 and its covariant 
derivative is defined as 


. Ta a . yl 
D,® Q = ig W, — ig TH ® 


Grae Re (2) 


—2(WitiW2) 0,4 $(gW2 - 9'B, 


in terms of the SU(2); and U(1)y gauge couplings g and 9’, respectively, the three SU(2);, gauge bosons 
W>, and the single U(1)y gauge boson B, with the usual three 2 x 2 Pauli matrices 


welt Doll D4 AY o 


Vem(®) = p’ (d'8) + A(S), 


Taking ® = (0,v + H)T/V2 with the vacuum expectation value (vev) v = \/—2/X and the real 
scalar field H after rotating away three Goldstone modes and using Wy (Wy; = iWe)/ V2 and 


Z, = (gW2 — 9'B,)/\/9? + 9”, we can render the kinetic term of the Higgs Lagrangian in Eq. (1) into 
the form expanded as 


(D"®)' (D,®) = 5 (OuHH)(O"H) + M?,W IW" + E M32," (5) 


1 1 M3 
+ gMw (wiw - zyz") VE (mywiw- -— az ze) H?, 


in the unitary gauge. We use the abbreviation sw = sin Ôw for the sine of the weak mixing angle Ow 
and cw = cos Ow, tw = sin 0w / cos Ow, etc. The masses of the massive gauge bosons W and Z are given 
by Mw = gv/2 and Mz = Mw /cw with v = (vV2Gr) 1/2 O 46 GeV fixed by the Fermi constant Gp, 
which is determined with a precision of 0.6 ppm from muon decay measurements (175) 116]. Incidentally, 
the SU(2)z and U(1)y gauge couplings are g = e/sw and g' = gtw = e/cw, respectively, where the 
magnitude of the electron electric charge e = 2,/7a with a being the fine structure constant. On the 
other hand, the SM Higgs potential takes the form of 


L 1 > 1 (3M2 s- de BME -a 
Vsm(H) = = vMi +3Ma H +3 ( 2 v H +7 = a, (6) 
which is completely fixed in terms of v and the Higgs mass My with the replacements of u? = —Av? 


and À = MẸ, /20". 


The Higgs interactions with the SM fermions 
interactions 


are derived by considering the following Yukawa 


-Ly = Uru QT (in) ® — DrbUUQ? (im) ® — Eghe LË (ir) ® + h.c., (7) 


where ® = inð* and QT = (UŁ, Dz) and LT = (v, Ez) with U and D standing for the three 
up- and down-type quarks, respectively, and v and E for the three neutrinos and charged leptons, 
respectively, in the weak eigenstate basis. And the 3 x 3 Yukawa matrices are denoted by h,¢,-. Taking 
® = (0,v + H)"/V2 again, we have 


—Lajp= >) fT, (8) 


f=u,d,c,s,t,b,e,u,7 
with the masses my = hy vu/ V2 in the fermion mass eigenstate basis diagonalizing the Higgs-fermion 


interactions. 


2.2 Complex Singlet Extension of the SM 


In this subsection, as the first BSM example, we consider a model in which the SM is extended by 
adding a complex SU(2);x U(1)y singlet (cxSM). 


2.2.1 Potential and mixing 
When a complex scalar singlet field S is added to the SM Higgs sector [40] [43} [44] [45] [46] [4°71 [48], 


the most general renormalizable scalar potential takes the form [40] 
V(®,S) = (ip) + ASD)? 

+ (5,8'S+cc.) + PtP |S? + (6,516 S? + h.c.) 

+ (aS+ec.) + (bS? +ec.) + be Is? + (Si + cc) + (eS IS)? + c.c.) 

+ (d,S'+ec.) + d |S|* + (d3 S° |S]? + c.c.) . (9) 
Imposing a global U(1) symmetry eliminates all terms containing complex coefficients. One may allow a 
soft U(1)-breaking b; term to avoid a massless CP-odd Goldstone boson which is not phenomenologically 
viable. And then, in order to avoid the cosmological domain wall problem caused by the presence of the 


bı term, one may additionally include the linear a, term which breaks the global U(1) and a discrete 
Zə symmetry under S > —S. The resulting cxSM scalar potential takes the form 


V(®,S) = pP (DS) + ASID)? + iPP |S|? + bo |S]? + d2 [Sf + (a1 S + b1 S? + c.c.) , (10) 
in terms of the original couplings in Eq. (9), or, alternatively [I7], 
V(®,S) = u(t) + u? [S| + A, (G1)? + A2 |S|* + Az" |S]? + (a1 S +0, S? + cc.) , (11) 


in terms of a more systematic parameter set of 5 real parameters of Hio and à1,2,3 and 2 complex massive 
parameters of a, = |a;|e’%+ and bı = |b, |e’*. 


By parameterizing the SU(2); doublet ® and singlet S as 


Gt e's 
B= (4 wrati ) Sg are er ally (12) 


we obtain the following three tadpole conditions for minimizing the potential: 


1 
Ç + Ayu + A v=0, 


1 ; i A 
È + A202, + 543 v +2 Rebre”) Vs — 2 3m(bie”®) va + V2 Relace!) = 0, 


1 , , : 
Ç + A202 + 53 ye —2 Bebe”) va — 2Sm(bye”) v, — V2 Im(aye*) = 0, (13) 
with the abbreviation v2, = v? + v2. The mass terms of the scalar states are given by 
1 $ 
VexsM „mass — ae s a’) Me s! ’ (14) 
a’ 


in terms of a real and symmetric 3 x 3 mass-squared matrix Mĉ decomposed into the two parts: 


Xi +20? AZ3UUsq, 0 


M = AZUV sq Xo + 2r2v?, 0 
0 0 Xo 
0 0 0 
+ 2] 0 Re(bye*)cy, — Sm(bie” 59, —Re(bye*5) so, — Sm(dye"*)co, |. (15) 


0 —Re(bye?"*) so, — Sm(bie*5)co, —Re(bie”)co, + Sm(bie*) 52, 


T 


The two parameters of X; 2 appearing in the diagonal components of the first term are defined by 
La 2,1, 2 = 2 2,1, 2 
Xi = uí + Àv + 5A Vsa» Xo = u3 + Àz Ua + 5r3¥ $ (16) 


In Eq. (14), the primed scalar fields s’ and a’ are related to the original scalar fields s and a through 


the rotation 
s _ Cy “By s (17) 
a =S: Gy a)? 


with c, = cosy = Us/Usq and sy = siny = U,g/Vsa. Note that Xı = 0 always to have the non-zero 
vev of v, as can be checked with the first tadpole condition in Eq. (13). On the other hand, only in 
the U(1)-conserving case with both a; = 0 and bı = 0, X> = 0 to have the non-vanishing vevs of va 
and v, giving rise to the massless Goldstone boson a’. In this case, the singlet vacuum takes the U(1) 
symmetric vev of Us, = \/v2 + v2 while each of the vevs remains undetermined. 


In some cases, instead of the discrete Zə symmetry, a different discrete Z4 symmetry under the 
interchange S + S* is imposed. In this case, the U(1)-breaking part of the scalar potential reads 


(aS + bS? + c.c) 


z, = 2 [Re(a1) S + Re(b) (S? — A?)] , (18) 


with S = S + iA = [(v, + s) + i(va +a)] / V2. Note that, if the potential has the Z, symmetry, just two 
real parameters are sufficient for parameterizing the U(1)-breaking part of the potential. Assuming a, 
and bı to be real with no loss of generality, the tadpole conditions of the scalar potential become 


Xıv=0, (Xə + 2b1) vs + V2a, =0, (Xə — 2b1) va = 0. (19) 


Assuming v Æ 0, which forces X, = 0, the mass-squared matrix is simplified into the form 


2)? AsV Usa 0 
Me = AZUUsq Xo + 2A202, + 2b1 C24 —2bj 827 š (20) 
0 —2b1 S2 Xo = 2b1 C27 


When vs # 0 and va ¢ 0 guaranteeing so, # 0, the CP symmetry is spontaneously broken and all 
the three states mix. In this case, one may parameterize the potential with a set of 7 parameters of 
{A1, A2, A3, V, Usa; tan y, b1} with the relations 


1 
A = — À v? = 548 v? Xo = 2b1, ay = —2V2 biv, (21) 


sa? 


where the second relation is solved to give u3 = 2b) — A2 v2, — $A3 0’. 


When v, Æ 0 and va = 0, the angle y = 0, i.e. s’ = s and a’ = a and the scalar mixing occurs only 
between the two states of @ and s with the pseudoscalar mass-squared 


— 4b. (22) 


In this 2-state mixing case, the scalar potential can be parameterized with a set of 7 parameters of 
{A1, A2, A3, V, Us, a1, b1} with the relations 


2 
2 X, = —¥2% oh, (23) 
Us 


1 
B = rv? - z% 


sa? 


where the second relation is solved to give 13 = —V/2a;/v,— 2b, — A2 v2, — $A3 v*. Incidentally, if vs = 0 
in addition to va = 0, the parameter a, should vanish due to the second tadpole condition in Eq. 
and Xə = uô + A3v?/2 giving the squares of three masses as 


Mj = 27, M; = p5+A3u?/2+2b,, M2 = p3 + A3v?/2 — 2b, (24) 
and the 6 parameters {\;, A2,A3, V, u2, 0, } can be employed for describing the scalar potential. Various 


types of vacua in the cxSM are summarized in Table 


Table 1: Vacua in the cxSM imposing Z}: the parameters a; and bı are real and v Æ 0 is taken. From 


Ref. [117]. 


vacua Xə a possible set of inputs miscellaneous relations 
Va #0 & vs #0 2b; £1, 2, A3, U, Usa, tan y, bi} ay = —2/2byv, 
Va #0 & v = 0 2b; {A1, A2, A3; V, Va, b1} aı = l; Vsa = Va, Cy = 0,5, =1 
Vg =0 & vs #0 | —vV2a1/Vvs —2b1  {A1, Ao, Az, V, Ve, a1, b1} Usa = Vz = 16,20 
Ua = 0 & v, = 0 u2 + Agv?/2 {X1, A2, Az, V, HS, b1} ai = l; Usa =0,c, 91,8, 30 


Finally, without loss of generality, the orthogonal 3 x 3 mixing matrix O diagonalizing the real and 
symmetric mass-squared matrix in Eq. is defined through 


(d, Bl Ve z Ooi( Hi, H2, Hs); , (25) 
such that OT MO = diag(Mj,,Mz,,.Mj,) with the increasing ordering of My, < My, < Muy. 


2.2.2 Higgs-boson interactions 


The interactions of the three Higgs bosons with the SM fermions and the massive vector bosons in the 
cxSM are given by 


Lip = B ox l ghg HEE 


f=u,d,c,s,t,b,e,4,T 


3 
Luvv = g Mw (wiw- aoe a Z" ) Da Sue Hi 


1 
Lunvy = 3 (mewy + "7 z) 5 Uaa n (26) 
i, j=1 
with the normalized dimensionless couplings simply given by 
one = avv = Ogi » Guajvv T Ogi09; - (27) 


And the cubic and quartic couplings are given by the self-interaction term of the scalar potential: 
—Lsat = MVE? + 2dA3veae"s' + 2Agud(s” + a”) + Arov.8' (S? + a’) 
1 
de pe 2 d3¢°(s” +a’) As do(s + 28"a” +a’) 


3 3 
= > JH, H; Hy HH; H, + > JH, H; HH; HH; HH, (28) 
izj>zk=1 izj2k2l=1 


where the normalized cubic and quartic couplings of the three Higgs mass eigenstates are 


3 


ra = 2, {OciOpjOon} Gapy 
a<p<y=1 
3 
JH, H; HE = pp {0:0 OykOs1} Jaßyô » (29) 
a<B<y<ő=1 


with 7, J, k,l = 1,2,3, the cubic weak-eigenstate couplings 


Usa 
Good = Ms Goes! = 2A3——, 
Usa 
Joss = Yda'a’ = 2A3 ; Ys!'s's! = Ys'a'a’ = 4X» cs ; 
(30) 
and the quartic weak-eigenstate couplings 
Josep = U/4,  Gbds's' = Idda'a’ = AB, 
Js's's'sh 5 Gatsigig | 2 = Ga'a'a'a’ = A2 : (31) 


In Eq. (29), the expressions within the curly brackets {---} need to be symmetrized with respect to the 
indices 2, j, k,l and divided by the corresponding symmetry factors in cases where two or more indices 
are the same. For example, {O.;08;O0,,} can explicitly be evaluated as follows: 


1 
{0x08 On} = N; (OaiOsjOrn + Ogi0 grO75 + OajOBiOyk + OajOpkOyi 
+On¢.O0 9104; + OarOgjOri) ’ (32) 


with Ns = 3! = 6 when i = j = k, Ns = 1 when (i,j,k) = (3,2,1), and Ng = 2! = 2 in all the other 
cases. 


2.3 Two Higgs Doublet Models 


In this subsection, we give a detailed description of the models in which the SM is extended by adding 
one more SU(2);, doublet (2HDMs) while taking the same gauge group SU(3)c¢ x SU(2), x U(1),, as in 
the SM. 


2.3.1 Potential and mixing 


The general 2HDM scalar potential containing two complex SU(2); doublets of ®; and ®, with the 
same hypercharge Y = 1/2 may be given by 


Vnom = Li (®191) + 13(G)82) + mi (DiD) + miaa) 
+A (DiD)? + A2(B}H2)? + As(H]H1)(@}H2) + Ay(G} Hy) (8181) 
+A5(®1 2)? + AZ (O51)? + Ag(@]H1) (DİD) + AZDI 1) (O)4,) 
+A7(&}B2)(B} b2) + Az (LD2) (8181) , (33) 


3Here, the indices a, 6, y, and 6 count the Higgs weak eigenstates of ¢, s’, and a’ and the inequalities among them 
imply that the cubic and quartic terms in the Higgs potential are ordered in the weak eigenstates. 


10 


in terms of 2 real and 1 complex dimensionful quadratic couplings and 4 real and 3 complex dimensionless 
quartic couplings. With the parameterization of two scalar doublets ®; 9 as 


_ oT gp. aif 3 
b= ( seth tin) )} ae Ge sa oe 


and denoting vı = v cos 8 = veg and vz = vsin 8 = vsg with v = yv? + v3, one may remove 7, u2, 
and Sm(m7,e"4) from the 2HDM potential using three tadpole conditions: 


1 ; 
R = —v hg + 5388 + easg Relace) + s3Mye, 
1 À 
R = —v ost + 5a + esa Relare)] + Mis , 
2 
Sm(m?2,e%) = a [2 cgsgSm(Ase"*) + c3S3m(Ase®) + s3Sm(Aze)] , (35) 
with the square of the charged Higgs-boson mass 
R 2 aig 2 . . : 
M? = — A a ee dacgsg +2 cgsgRe(Ase) + c2Re(Age”) + s2Re(Aze*)] . 36 
H BB BB B B 
Cg58 2€888 


Then, including the vacuum expectation value v, in general we need the following 13 parameters plus 


1 sign: 


v, tg, [Mmi]; 
Ay, A2, A3, A4; fAs] ’ [Ael ’ |Az|; 
ps F 2E ’ 06 ag € ’ o7 F E ’ sign{cos(¢12 + £) , (37) 


to fully specify the general 2HDM scalar potential in a form given by Eq. (33). Here mî, = |m?,| e1012 
and A567 = |As.6,7|e'?%7 and we note that sin(ġ12 + £) is fixed by the CP-odd tadpole condition 
if the CP phases $5 + 2€,¢6 + € and ¢7 + € are given together with |myo|, |A5s,67|, v, and tg and, 
accordingly, cos(¢12 + £) is determined up to a two-fold ambiguity. One may take the convention with 
€ = 0 corresponding to re-defining the 1 quadratic and 3 quartic complex parameters without loss of 
generality 


The 2HDM Higgs potential includes the mass terms which can be cast into the form consisting of 
two parts 
1 Qı 
V2HDM „mass — M? HTH- F 5 (ht Q2 a) Me Q2 , (38) 
a 


in terms of the charged Higgs boson H~, two neutral scalars ¢1,2, and one neutral pseudoscalar a after 
absorbing the charged and neutral Goldstone bosons G~ and G® in the 2-state mixings of the two 
charged scalars and two neutral pseudoscalars as 


T S SCS). 


And the 3 x 3 real and symmetric mass-squared matrix of the neutral Higgs bosons M@ is given by 


s%  —sgcp 0 
Mj =Mi|{ -sses ch 0 | + Mj, (40) 
0 0 1 


4Note that all the parameters are neither basis independent nor physical. There are only 11 physical degrees of freedom 
in the potential as counted in, for example, Ref. [59]. 


11 


with (reinstating the relative phase € for the sake of generality) 


M2 = Mia + [534 — Re(Aee*) v, (41) 


and the second part expressed in terms of the quartic couplings as 


2102 + We(Ase?"*) 52 2Azacgsg + Re(Age*)c3 + —Sm(Ase?"*)s 
b B pap b b 
S 


+2Re(Age"*) ace +Re(Aze") 5% —Sm(Age"*)eg 
Mi — | 2Azacasg + Re(Age)cZ 2283 + WRe(Ase™)cZ —Sm(Ase?®)cg (42) 
v? +Re(A7e"*) 83 +2Re(Aze") s Cg Sm (Aze"*) 54 
—Jm(Ase”"*) 8.4 —Sm(Ase”"*)eg 0 
—Sm (Age )cg —Jm(Aze)sg 


where the abbreviation A34 = (A3 + A4)/2 and, in passing, we note v = 2My/g, a = —sgaı + cgaz and 
Ht = —sg¢} + cad. We need to specify, therefore, the 13 parameters plus 1 sign listed in Eq. to 
fix the mass-squared matrix. 


Once the real and symmetric mass-squared matrix M@ is given, the orthogonal 3 x 3 mixing matrix 
O is defined through 


(d1,62,a)¢ = Oni(Hi, Ho, Hs); , (43) 
such that OT M20 = diag(Mz,, Mj, Mzi,,) with the increasing ordering of My, < My, < Mps- 


assive vector bosons 
The cubic interactions of the neutral and charged Higgs bosons with the massive gauge bosons Z and 
W~ are described by the three interaction Lagrangians: 


Luvv = gMw (wiw- Pap sar 2") Da Suv Hi 


oO 
Lanz = —— Sa Ou 
2cw i>j : 
oO 
Lage = E Y Ipan Wo" (Hii dy Ht) + hc., (44) 


++ 
respectively, where X 9, Y = X0,Y — (0,X)Y, i,j = 1,2,3 and the normalized couplings Daan 
Gn,n,2 and g, ,+w- are given in terms of the neutral Higgs-boson 3 x 3 mixing matrix O by (note that 
det(O) = +1 for any orthogonal matrix O): 


Guvy = Cada + 8p Ogi, 
Jaa = sign|det(O)| Eijk Iu,vv > 
Guntw- = C8 Odai — 86 Ogri — Oai, (45) 
leading to the following sum rules: 
3 
` lawy = ] and A T Cee ? = 1 for eachi=1,2,3. (46) 
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On the other hand, the quartic interactions of the neutral and charged Higgs bosons with the massive 
gauge bosons Z and W= and massless photons are given by 


1 M2 : 
LHHVV = J2 (mwm T 2,2") ` Inn, vv HiH} , (47) 
ij=l 


with Gu,Hjvv = Oi and 


2 2. 72 
Lutu-vv = (4 wiwe + ZB 207, + AMA, + eg2 chy ANZ HH, 


2 3 
GZ 98 = 
Lusnawe = 2 i (zw "S lana tne) ? 
i=1 
eg : 
Ly+Hawt = = (AAS tanner tt ne) (48) 
i=1 
with lamna =O pct Ogi — Cp Ob. — 10a, Cow = cos 20w, and gz = g/cw = e/(swew). 
2.3.3 Interactions of Higgs bosons with the SM fermions 


Without loss of generality, the Yukawa couplings in 2HDMs could be cast into the form [L18]: 
Ly = h TRR" (ite) ®2 — hadr QT (iT) (nf @, + ng ,) 
— hlr IL" (ite) (ni d + 7, d») + h.c., (49) 
where ®; = inž and QT = (uz,dz) and LT = (v,,lz) with u and d standing for three up- and 
down-type quarks, respectively, and / for three charged leptons. We note that there is a freedom to 


redefine the two linear combinations of ®) and ©; to eliminate the coupling of the up-type quarks to 
®, [119]. The 2HDMs are classified according to the values of 7. and nf, as in Table 


Table 2: Classification of 2HDMs satisfying the Glashow-Weinberg condition which guarantees 
the absence of tree-level Higgs-mediated flavor-changing neutral current (FCNC). 


2HDMI 2HDMIL 2HDMII 2HDMIV 
ng 0 1 0 1 
ne 1 0 1 0 
ni 0 1 1 0 
n 1 0 0 1 


By identifying the couplings in terms of the vev v and the mixing angle 8 as 


h J2m, 1 A Vma 1 h J2m, 1 (50) 
= —— d= ; = , 
v sp v nice + nise v nics + mse 


5Here we take the convention with € = 0 and the couplings hy,a are supposed to be real without loss of generality. 
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we obtain the following Lagrangians 


u O i . 
—Ly,ff = di c (S-i Boun ) | H; 


v 88 
= (Opi t Opi .nisa — nd 
p Ma a (2 gi Ta gai ge F207 | i H; 
v nice + 288 mC + 1288 
B lO, . LO... lsg — n! 
| MI i (2 ou + 2 got i du n Oy 18) 1 H; (51) 
v 1168 + Msg Ca + 1S, 


for the interactions of neutral Higgs bosons with fermion pairs. For the interactions of the charged 
Higgs boson with fermions, 


—Lutiad = — 


v2 (2 

5B 

E V2m ( nise — niceg 

v \ mes + nss 
where Pr = (1 + 5)/2 and Pr = (1 — 45) /2. 


) ūP,d H+ = V/2ma (a — n§ce 


ti Prd Ht 
v nice ae) 


) vPrRlH* + h.c., (52) 


2.3.4 Higgs-boson self-interactions 


Given the orthogonal mixing matrix O diagonalizing the mass-squared matrix of the neutral Higgs 
bosons, the cubic and quartic Higgs-boson self-couplings are given in terms of the Higgs mass eigenstates 


by [121] (122) (123) (124): 


3 3 
-Ln = v > (ee E ee H,H+H-, (53) 
i>j>k=1 i=l 
3 3 
Li = ` Taani HoH Hey + >D Donans eg A 
i>j>k>l=1 izja 
iO I ee (54) 


where the normalized cubic and quartic weak-eigenstate couplings are ý 


3 3 
Gun Hy, = D {0GO0 gO} Japy > Ia H+H- = X Ox U ayd (55) 
a<p<y=1 a=1 
3 
JHH; HH = Lp {Oxi O8; OykOsi} JaBy6 > 
ax<p<y<b=1 
3 
J a,n 1+ = ` {OaiOp;} Jaßut a= - (56) 
a<ß=1 


We note again that, in Eqs. and (56), the expressions within the curly brackets {---} need to 
be symmetrized fully with respect to the indices 7, j, k,l and divided by the corresponding symmetry 
factors in cases where two or more indices are the same as in, for example, Eq. (32). 


Here, the indices a, 8, and y count the Higgs weak eigenstates of ¢,, #2, and a and the inequalities among them 
imply that the cubic and quartic terms in the Higgs potential are ordered in the weak eigenstate basis. 
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For the sake of completeness, we present all the effective cubic and quartic Higgs-boson self—couplings 
of the Higgs weak eigenstates. The cubic self-couplings of the neutral Higgs bosons are given by [] 


1 
Cài + 5 sp Rede, 


Isiti = 
3 
Yoreidg Z Sp A34 + sg Reds + 5 C2 Rede, 
3 
Is; 60¢2 = CB A34 + C8 Reds + 3 8 Redz, 
1 
VEA = SgA2 + 9 Cg Redz, 
1 
Ióigja = —sSgep SmA5 E 9 (1 + 2c3) Småge j 
Joiosa = —23MAs — Sgcg Sm (Ag + Az), 
1 
Dba = —SgCB SmAs5 = z (1 + 28%) SmAz , 


1 1 
Ipaa = sgcgd1 + ch Asa — ce(1 +53) Reàs + 5 80(83 — 2c) Reds + 58863 Rerz , 


1 1 
Genoa = 83C3A2 + s3 A34 — sg(1 + c3) Reds + 59808 Reds + 5 olch — 283) Redz, 
1 1 
Jaaa = Seca SMAs — 5 88 SmAg — 58 JmA7 , (57) 
with the abbreviation A34 = i (A3 + Aa). The effective cubic couplings J autu- read: 
Dn = 2s3CBÀ1 + ce3 — SgCpàa — 28%cg Reds + splsg — 2c3) Rede 
+ sgh Redr , 
a Qsgcgr2 + 83A3 — 8gc3Aq — 2sgch Reds + s3cg Rede 
-+ calch — 253) Redz, 
I yru- Z 28eCa SMÀs — s% SmAg — c3 Småà7. (58) 


On the other hand, the quartic couplings for the neutral Higgs bosons are 


1 1 1 1 
Jsigigidi = 4 Àr ’ Doi dvd = z Rerg , Joididaga = 2 A34 + 2 Reds , 
1 1 
Is 62020 = 2 Rez ? I 42600260 = 4 A2 , 
1 1 
Dinne =O Småge , Issa = -Ce Ms — 5 S8 SmA¢ , 
1 1 
Doi dotaa = —§~8 SmA5 a 2 Cg SmAz ’ Dbsdo daa = T3 SB SmA7 3 
1 2 1 2 1 2 1 
9$,¢,00 = 3 S$ At + 9 8 A34 E 3B Reds z z SPCL Rede, 
1 1 
Terogaa = 258C Reds + 5 s% Res + 5 c3 Redz, 
l 2 1 9 1 2 1 
Téscdauts = 3B A2 + 3 8 À34 = z 88 Reds = 5 56Ce Redz, 
1 1 
G,aaa E 58C3 SmAs5 D 2 S3C SmAg ag 2 ch SmA7 5 


"The relative phase £ could be reinstated by replacing As with Ase7"§ and Ag.7 with A¢,7e"s, if necessary. 
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1. 1 
= 53g Små; — 5 88 SJmAg — 5 86C8 JmAz , 


1 


O aaa = J Jat nH n+ H- ? 


U izana 
(59) 
together with the quartic coupling of the charged Higgs bosons given by 

Oe ShA + CBA + s3chlÀs + As) + 2s3chReàs — s}cpRerg — 2sgcReàz. (60) 


Finally, the remaining quartic couplings involving the charged Higgs boson pairs, gagu+n-, are given by 


1 
= 834 + 5 C88 _ 53Cg Rede , 


$19, HTH- 
suta- — SACP Aa — 2sgcg Reds + s% Reds + c3 Rerz , 
o gutam T c3 Ag + 1 Ag — sgeg Redz, 
Orare ~ 25c SMAs — s3cgSmAg — cZSmAz, 
Ianin- = 255Cp SMAs — s3 SmAg — sgcgSmz, 
Janta- T Gutu-ntu- * (61) 


2.4 Minimal Supersymmetric Extension of the SM 


In this subsection, we give a brief review of the minimal supersymmetric extension of the SM (MSSM) 
with our particular focus on the MSSM Higgs sector. For general reviews on the MSSM, see Refs. 


[61] (621 (63) [64 [65] [66} [67 [68] [69]. 


2.4.1 Potential and mixing 
potential of the MSSM 


The super is written as 


Wyussm = UCh,Q- A, + D°hafa Q + Ehf- L + pfu Ba, (62) 


where Aaa are the two Higgs chiral superfields, and Q, L, U a D? and E® are the left-handed doublet 
and right-handed singlet superfields related to up- and down-type quarks and charged leptons. The 
Yukawa couplings hua, are in general 3 x 3 complex matrices describing the charged-lepton and quark 
masses and their mixings. The superpotential contains one supersymmetry preserving mass parameter, 
the u parameter that mixes the two Higgs supermultiplets, which has to be of the electroweak order 
for a natural realization of the electroweak symmetry breaking mechanism without any significant fine 
tuning. See Table [B] for the full particle contents appearing in the MSSM. 


In an unconstrained version of the MSSM, there are a large number of different mass parameters 
present in the soft SUSY-breaking Lagrangian 


1 ae ates ee ee eee 
Lot = 5 (m BB + M WIW’ + Myg7g? + h.c.) + M20 + ÜM + OtM2,o 
+ DIM2D + EME + M2 HiH, + M3, HtHa + (Bu Hu: Ha F hie.) 


$ (Ua.Q- Hy Dar OQ ti arel h.c.) (63) 


8Here and in the following, we introduce a product between SU(2); doublets as defined by A- B = ATit3B = ceap A” B? 
with €12 = — €21 = +1. 
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Table 3: Particle contents in the minimal supersymmetric extension of the SM (MSSM). Note the 
relations Q = 73 +Y and W+ = (W! = iW?)/V/2 and W° = W. 
spin-0 spin-5 spin-1 color T T3 Y B L 
(5) e) 3 5 o Ta +3 10 
Ü; (u°), x fof o | -3 | -4 Jo 
Ds Can af 0 0 +3 | -4 | 0 
7._{ % VE 1 +1/2 ae 
in=( 2 | tn=( x | 1 5 a 5 0 1 
D (E°), 1 0 0 +1 | 0 |- 
g g 8 0 0 0 0 0 
W+ wt +1 
w° w° 1 1 0 0 0 0 
W- W- —1 
B B 1 0 0 0 o | 0 
m eo a H? i 1 e? E i i 
HY \ x Hy +1/2 
Hou) = 2 Ao =| 2 1 1 1 0 0 
2(u) ( H? ) 2L ( H? ) 2 ( —1/2 +5 


Here M93 are the soft SUSY-breaking masses associated with one U(1)y gaugino B, three SU(2), 
gauginos W’, and eight SU(3)c gauginos 9°, respectively. In addition, Miz, , and By are the soft 


masses related to the Higgs doublets H, a and their bilinear mixing. Finally, M% z pvp are the 3 x 3 
soft mass-squared matrices of squarks and sleptons, and a,4,- are the corresponding 3 x 3 soft Yukawa 
mass matrices. P Hence, in addition to the u term, the unconstrained CP-violating MSSM contains 109 
real mass parameters including 46 CP phases: 


M23 OX a2 (1) = 6 (3) , 
Mê 1DE 5x 9(3) = 45 (15), 
aude 3 x 18 (9) = 54 (27), 
Min, , 2x 1(0)= 2(0), 
B : 2(1)= 2(1); (64) 


where, in each line, the number of CP phases is separately counted in parentheses for the corresponding 
soft-SUSY breaking parameters. A caution is made in order to avoid misinterpreting that the 109 


° Alternatively, the soft Yukawa mass matrices a4 ae may be defined by the relation: (Quctalig = Muae)iy (Antic): 
where the parameters (Axa,¢)ij are generically of order Msysy in gravity-mediated SUSY breaking models. 
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parameters are all physical. The number of new physical parameters is actually 105. Including them, 
the MSSM contains in total 124 physical parameters, as counted systematically in, for example, Ref. 
65|. 


By identifying H, = ®, and Hg = ©, = inð* = (¢* ,—¢7)T, fo] one may obtain the same form of 
the Higgs potential as in the 2HDMs: 


Vussm = Mili D) + p3(Ghb2) + mi,(G1G2) + miD) + Ar(G}G1)? + A0) 
+ A3(®1®1)(O}by) + AlDID) (DiD) + As(PiD)? + A (Dia)? (65) 
+ \g(b16,)(16,) + rz(016,)(O16,) + \7(O562) (DiD) + A2(O}6,)(G56)) , 


with the potential parameters given in terms of the u and the soft SUSY-breaking parameters as well 
as the SU(2); and U(1)y gauge couplings by [E] 


4 = Mintle, m = Mi, tle’, mi, = —Bu, 
1 1 ; 1 
` = = gig +9"), hg 0): M = =r; 
w = ea’ ees (66) 


Note that the quartic couplings À1,2,3,4 are solely determined by the gauge couplings and A567 are 
vanishing at the tree level. However, the quartic couplings 5.67 receive significant radiative corrections 
from scalar-top and scalar-bottom loops and, especially in the presence of CP-violating phases in the 
soft SUSY-breaking terms, the CP-violating mixing among the three neutral Higgs states are induced 
[125] [126] [1277 [128] [129] [130] [31]. In this case, as in the 2HDMs, the orthogonal 3 x 3 mixing matrix 
O has to be introduced for diagonalizing the 3 x 3 real and symmetric mass-squared matrix of three 
neutral Higgs states through 


(1, G2, a) 4 = Oui( Hi, Hz, Ha); (67) 


with the increasing ordering of My, < My, < Mug. 


2.4.2 Interactions of Higgs bosons with the SM particles and self-interactions 


The interactions of Higgs bosons with massive vector bosons and those among themselves in the MSSM 
are formally the same as in the 2HDMs. 


The tree-level Higgs-boson interactions with the SM fermions are the same as those in the type-II 
2HDM without including the finite loop-induced threshold corrections mediated by the exchange of 
gluinos and charginos [37]. Including the threshold corrections, the couplings 
of neutral and charged Higgs bosons to down-type fermions could be significantly modified for large 
values of tg. More explicitly, by resumming potentially large tan -enhanced effects due to the threshold 
corrections, the down-type quark Yukawa couplings may take the form 


_ Wiig, 1 


h 
UCB 1+ Ay te 


(68) 


KA 


T , T 
10We recall ©; = ( Ta z +o, + ian)) and Po = e (9, Z (vo + do + ia2)) together with a = 
—8ga1 + Cga2. 
Note that Hu -Hg = HT itz Ha = -PT 0% = — 616). 
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where q, = d,s, b and the tg-enhanced threshold corrections enter through the one-loop quantities 


205 * * 
Aa = 3n ” M; I( M3, M5, Mal’), 
2s in 
A; zn V Ms (Ma M5 Mal): 
b 3g K 3 Qa Da? 773 Ir t Ga Da P > 


with the loop function I(a, b,c) defined as 


ab ln(a/b) + bc In(b/c) + acln(c/a) 
(a — b)(b—c)(a—c) 


In the presence of CP-violating mixing in the neutral Higgs sector, the effective Lagrangian describing 
the Higgs interactions with up- and down-type quarks is given by 


3 
Ma J , 
Lug = — 5 p P |a (95.5, F igh 4,75) ql H; 


I(a,b,c) = (70) 


i=1 q=4¢4,9) 
V2 Mmg 

tl- [S ` P + 
uD i ja (03... igh, 75] a Ht +he.p, (71) 

where (qy, q) = (u, d), (c, s), (t,b). At the tree level, as in the type-II 2HDM, 

1/1 m i1 m 

T — Opi SO. gS = ae ea, A =e ee, 72 
I ga diay CRs Iqa BVais Jata 9 + Bl > Jatma 9 is Ma, Bl > ( ) 


for the neutral and charged Higgs bosons, respectively, see Eqs. and (52). While, in the presence 
of tg-enhanced threshold corrections, the couplings g® g? g dg’ are given by [138] 
Gay 


Hyda,? 7 Ayayay’ H+ 
[139] [140]: 


hn, HB) te) ee, 


Ht aya, 


1 + Aate’ €p 1+ Aytg/ cg 1+ Ag, ts 
IH = - Re ( L) Ou F Sm ( e ) pi Sm ( dN ) 2 
1+ Ay ts 1+ Ay tg Cg T- Áj F Ca 
g5 — 1] 1 ds Ma, tg = Ag g? = a} 1 2 Ma, tg = Ai (73) 
Htapay 2|tg mg \1+ Aj ts mtma 2 |ts mg, \1+ As te 


Note that the above couplings approach to the tree-level ones in Eq. in the limit of A,, = 0. The 
size of the resummed threshold corrections could be significant enough to make the down-type quark 
Yukawa coupling hg, as comparably large as the top-quark Yukawa coupling when Re(A,,) < 0 and 
tan 8 is large. Some extreme cases with |,| ~ 1 and |h,| ~ 1 are discussed, for example, in Ref. 
and Ref. [140], respectively, taking account of the CP-violating mixing in neutral Higgs sector. 


!2We note that the one-loop quantities are uncertain at the level of about 10% depending on scale choices for as. 
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2.4.3 Interactions of Higgs bosons with the SUSY particles 


For the sake of completeness and reference, although no serious analyses on them are presented in the 
present review, we fix the convention for the interactions of Higgs bosons with the supersymmetric 
(SUSY) particles such as charginos, neutralinos and sfermions. 


The interactions of neutral Higgs bosons with charginos, which are mixtures of charged gauginos 
and higgsinos, are described by the following Lagrangian: 


2 3 
g =_ [9g ; P f= 
Lmg- = “a yD AX; Ce Fiona) Xj > (74) 


ij=1 k=1 


with the normalized Higgs-chargino-chargino couplings 


date = 5 {lCadalCrnGt + (Cala(CynGel+ to a}, 
dharar = 3 {U(Cr)alCr)eCh + (Cr)a(CrpGel — feo a}, (75) 


where G$! = (Osik — i8gO0ax) and G? = (Opak — icgOax). The two different unitary 2 x 2 chargino 
mixing matrices (Cz)ia and (Cr)ia are required to diagonalize the chargino mass matrix 


M: V2Mw c 
Me= , " ’) (76) 
V2My sg H 


in the (W~, H~), and (W+, H*),, bases with the convention Hy; = Hy; and H} = Hj in such a way 
that 
CrMcCy, = diag{mys, myx}, (77) 


with the increasing ordering of Met < Mg. Explicitly, the mixing matrices relate the electroweak 


eigenstates to the mass eigenstates, via 


We = J 272 >) O 


i=1,2 i=1,2 
We = > Coit, eH >) (Chau (78) 
i=1,2 i=1,2 


Note that the convention Hy, r) = Ay is adopted with the subscripts 1 and 2 being associated with 
the Higgs supermultiplets leading to the tree-level mass generation of the down- and up-type quarks, 
respectively, see Table We recall that we take the following abbreviations throughout this paper: 
sg = sin p, cg = cos Ê, tg = tan p, S28 = sin 26, cog = cos 26, sw = sin Ôw, cw = cos Oy, etc. 


The interactions of three neutral Higgs bosons with neutralinos, which are mixtures of 2 neutral 
gauginos and 2 neutral higgsinos, are described by the following Lagrangian: 


4 3 
g = : Ag 
L Hogogo = =) De Y AR? Gar F ins9hs030) x (79) 


with the normalized Higgs-neutralino-neutralino couplings 


1 * * * 1 * . . 
ARR = z tel(Nj2 = tw N») (NG? = NAGS) +(e j)], 
1 * * * * . . 
IH — ~ 9 Sml(Nj2 = tw N5 (NAGI = NGH) +(i e j)], (80) 
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where i,j = 1-4 for the four neutralino states and k = 1-3 for the three neutral Higgs bosons. One 
unitary 4 x 4 neutralino mixing matrix is required to render the 4 x 4 symmetric neutralino mass matrix 
expressed as 


Mı 0 —Mzcgsw Mzsgsw 
0 Mo Mzcgcw —Mzsgcw 
—Mzcgsw Mzcgcw 0 —L 
Mzsesw —Mzspcw =H 0 


in the (B, W°, H°, 19), basis into a diagonal matrix as 


N*MwyN' = diag (m mgo, Mx, Mz, MZ) , (82) 


with the increasing mass ordering of mgo < mgg < mg xy L mgo. The single neutralino mixing matrix Nig 
relates the left-handed and right- bended slectroweale m ates to the left-handed and right-handed 
mass eigenstates via 


(B,W*, A, A = Naa Xe Xela, and 
(BW HE He ee = Ni PGCE cmc eine (83) 


respectively. 


The interactions of the charged Higgs bosons H~* with charginos and neutralinos are described by 
the following Lagrangian: 


4 2 
L ya xoyt = a ` 5. Hx ra eee + ined rs s09- ) uy + h.c., (84) 
i=1 j=1 
with the normalized couplings of the charged Higgs boson H* with a chargino and a neutralino 
Iiz = ; {ss [V2N5 (Cr) =(N + twNa)(Cr)}o| 
+ Cg [VZN (Cr); + (Nig + twNan)(Cr)ja| \ , 
I+ 30% = - {ss [VZN (Co) — (Ni + tw NA)(Cr)al 
— cg [VZN (Ca); + (Nig + twNan)(Cr)ja| \ ; (85) 


expressed in terms of the chargino and neutralino mixing matrices. 


The neutral Higgs-sfermion—-sfermion interactions can be written in terms of the sfermion mass 


eigenstates as 
3 
Lagp=0 Dd 2 Ing Hifi fe), (86) 


feud i=1 j,k=1,2 


where the couplings of the Higgs bosons with sfermions in the mass eigenstate basis 


Z-— + >» oo. Oot U5, (87) 


a=1,¢2,a B,y=L,R 


expressed in terms of the 2 x 2 scalar-sfermion-sfermion coupling matrix refs in the weak-eigenstate 
basis, of which the explicit form is given later in Eq. (94), and the 3 x 3 Higgs and 2 x 2 sfermion mixing 
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matrices, O and Uf, with the convention of a = (¢1, ¢2,a) = (1,2,3), 8, y = L, R, i = (Hı, Ho, H3) = 
(1,2,3) and j,k = 1,2 chosen properly for the sake of notational convenience. Likewise, the charged 
Higgs-boson interactions with up- and down-type sfermions are given by 


Lpp =t. > S Impp, At F Pa) the. (88) 
(Ff )=(u.d), (v1) j,k=1,2 


where the couplings of the charged Higgs boson H* with sfermions in the mass eigenstate basis 


a HHEP F f! 
B, y=L,R 


expressed in terms of the 2 x 2 charged Higgs-sfermion-sfermion ao matrix TH+] in the weak- 
eigenstate basis, of which the explicit form is given later in Eqs. and , and the 2 x 2 sfermion 
mixing matrix U f with the convention of 6, y= L, Rand j,k = 1, > The unitary 2 x 2 sfermion mixing 
matrix US is obtained by diagonalizing the 2 x 2 sfermion mass matrix Mi for f = t,b and 7 in such 
a way that 7 7 

U M5 US = diag(m2 ,m*),, (90) 
with the increasing mass ordering of me < me. The mixing matrix U I relates the sfermion electroweak 


eigenstates r to the sfermion mass eigenstates fe via 


(ft, fala = Vas (fi, fo) - (91) 
Explicitly, the stop and sbottom mass-squared matrices are written in the (qz, qr) electroweak basis as 


M B Mo, + me a cop Mz (TS — QaSiv) Ayu (A; = pR v2 (92) 
q=t,b A — H *R a) /V2 Mz, + m? + cogMZ Qeaey ’ 


with the third-generation left- and right-sfermion soft SUS Y-breaking mass-squared MS 2 „and Mi, aes 


a cubic soft-breaking term A,, Tt = —T? = 1/2, Qi = 2/3, Qe = —1/3, v = 11, Vi = on R, = tan 8 = 
v2/v1, Ry = cot 8, and the Yukawa coupling h, of the quark q. Similarly, the stau mass-squared matrix 


is written in the (7z, Tp) electroweak basis as 


vom ( M2 +m? + cxgM3 (sły — 1/2) hgui(At — utan B)/V2 ) a 


hu(Ar— u" tanp)/ VI M2 + m2 = cgM3 siy 


derived directly from the sbottom mass-squared mass matrix by replacing b by 7 and Qs, Ra by 
L3, Es, and taking Q, = —1. Incidentally, the mass of the tau sneutrino v; is simply given by 


= v M? + ico M2, as it has no right-handed counterpart in the MSSM unlike the squark and 
3 


charged slepton cases. 


my, 


r 


For the sake of completeness and explicit analytic and numerical calculations, we present the ex- 
plicit form of the Higgs-sfermion-sfermion couplings in the electroweak-interaction basis for the third- 
generation sfermions. The 2 x 2 coupling matrices [°F are given in the (fz, fr) basis with f = t, b, T, v; 
and a = a, ġ1, Q2 by 


E E Bloods +e \ 
V2 —4 ha(sg Ap + Cpt") 0 ? 
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pee _ —|he vce + 4 (9? + $97) veg ship Ay 
-zA =|| veg + gg? veg |’ 
1 1 1p 
pedd -4 is? + 39°) VSB ya 
Bhon” — 49" 08g ’ 
pet _ t 0 i hý (cg A% + sgu) 
2 —1 hi(cgAt + Spit") 0 , 
1 1 i 
pate _ | 4 ( AF 39°) Vp zaan 
Sale —$39' cp 3 
pete = —|h,|?vs, + F (9 = 39°) USB -ah Ai 
— sah =|| vse + 4g°vsg J’ 
parr _ i 0 ih; (spAt + cpu) 
V2 \ —ih,(sgAr + cpp") 0 , 
parr _ —|h,|Pveg + i (g? — g°) veg -aA 
— —gghrA, —|h,[?vcg + $9vcg J’ 
pert _ —4 (9? — 9?) usg alted 
alr —3g“vsg )’ 


a a yr 1 SAD 
[ers = 0, [eres = z (9° Ze g’) UCB, [evr = 


A| = 


The 2 x 2 coupling matrix pH is given in the (wz, dpr) basis by 


pated _ Z (Rul? + [hal]? — 9°) vspeg hå (3344 + cph) 
hu (CgAu + Sept") A h,hġ v 
and the couplings of the charged Higgs boson with a tau sneutrino and a stau given by 
1 TORT * * 
= (|b? —9?)ussee, Te = ht (spA* + emi). 


V2 


peat 


3 Decays of a Generic Neutral Higgs Boson 


(9° + g’) USB. 


(95) 


(96) 


Without loss of generality, the Lagrangian describing the interactions of a generic neutral Higgs boson 
H with two fermions, which is applicable for all the models described in the previous section, can be 


written as mn : 
Luff = ae Af (987 T ishg) Í, 


(97) 


in terms of the normalized scalar and pseudoscalar couplings of g3, F and gh rf with ms denoting the 
fermion mass and v ~ 246 GeV. The Lagrangian describing the interactions of the neutral Higgs boson 


H with massive gauge bosons Z and W= can be written as 


= 1 
Lavv = g Mw (a WEW 5 F Inza zy ZZ") H; 
W 


(98) 


in terms of the normalized couplings of gyww and g,,, with g = e/sw the SU(2)z gauge coupling, 
sw =sinOw, cw = cos Ow, tw = sin 0w / cos Oy, etc. And, if not mentioned otherwise, we set gyww = 


Juzz = Guvy ID the following. 


23 


In the presence of Higgs/scalar bosons y’s lighter than H, the neutral Higgs boson H can decay into 
a lighter Higgs boson and a massive vector boson and also into two lighter Higgs/scalar bosons. The 
interaction Lagrangian describing these types of decays could be cast into the expressions: 


g ae 
Lupz = Dey Juez S (H 0,9), 
g HKG. ot 
LHotwt = = See WF Oye”) + h.c., 
Leat 2-0 S aga HPP) = 0 a TT Oa O Fa H¥R) ++] - (99) 


izj 
Note that the scalar states y; and y; are ordered in the last expression so as to avoid the couplings 
such as gHy,¢, With a wrong ordering of the scalar states. 


It is noteworthy that we are assuming the neutral Higgs boson H to be a general CP-mixed state, 
i.e. a scalar-pseudoscalar mixture. 


3.1 Decays into two fermions: H > ff 


Including the radiative corrections known up to now, the Higgs decay width into fermions can be 
organized as [87]: 


F 7+ BiMu [a s 12 FS sf f 
(A > ff) = Nowe 87 [Bosz (1+ dcp + 6 + Sf nea) (1+ sfa) 


elal (1 + dacp + af?) (100) 


where Bf = 4/1 — 4ks with Kp = M/M} and the color factor NÉ = 3 for quarks and 1 for leptons. 
The lepton pole mass is taken for ms while, for the Higgs decays into quarks, the MS quark mass 
™,(My) is used. In passing, we recall that 1/v? = 2G. Note that, in the scalar part, the QCD 
and electroweak (ELW) corrections are factorized. This factorization is supported by the reduction 
of the mixed corrections by a factor of 3 [141]. In the pseudoscalar part, we neglect the ELW cor- 
rections. A pseudoscalar component appears in BSM models and the corresponding ELW corrections, 
compared to the QCD corrections, are much more complicated and depend on specific BSM models 
under consideration. In this review, we do not consider the ELW corrections for the decay processes 
involving multi-Higgses, CP-odd or charged Higgses sacrificing the precision in the calculation of the 
partial widths for those decays. For the calculation of full BSM-dependent ELW corrections and also 
of more precise SM ones, we provide Appendix |F|in which we make a brief introduction to relevant 
numerical packages. 


The pure QCD corrections for the Higgs decays into a quark pair qq consist of a universal part of 
dacp and two types of flavor- and parity-dependent parts of 67° and 5%" which are given by 
(144) [145] [146] (147) 048, [149] (150) (151) [1521 [153] (154) 155] 


s(M (M 2 
dacd = porta + (35.94 — 1.36NF) (a) 
g T 
(Mp) y 
+(164.14 — 25.77Np + 0.259N2) (at) 


(Mg) \* 
+ (39.34 — 220.9N r + 9.685N7 — 0.0205.N3) (=) , 


13In this review, we denote the pole mass of the fermion f by My and its running mass by m+ (j1). 
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S 2 =, 
qs JHH a.(My)\? 2 Mj, 1, .™m,(Mu) 
= 1.57 — = log l 
Ò; ( = 57 8 T? 9 og M3 , 


IHG 
P 2 2 mye 
: Jun ( %s( Ma) My 1, m (Mn) 
of? = Hu ( ) [3.83 — toe + — log — , (101) 
' Iirg T M? 6 M2 


where Np counts the flavor number of quarks lighter than H. The QCD coupling strength a, and 
the running MS quark mass 7,(Mj) are defined at the scale of the Higgs mass to absorb large mass 
logarithms. 


For the electroweak corrections [156| 158} [159], we adopt the approximation [160| 


Mi G 
nae = — lo +—= at M? + MP |- 5+ 


3 
> los cy — 8 MZ (6v% 75 — Ly A 


sw 

(102) 

where vz 75 = Die Qgrsiy and azz, = If /2 with If denoting the third component of the electroweak 

isospin and Qp the electric charge of the fermion f. We refer to Eq. for the Z couplings with 

fermions. The large logarithm log M?, /M? can be absorbed in the running fermion mass as in the QCD 

corrections. For decays into leptons and light quarks, the coefficient ky = 7 while it is 1 for b and t 

quarks. The electroweak corrections are below the percent level for f = b,c while they are of O(1-5)% 

for f = 7,p. For the more precise evaluations of SM and full BSM-dependent ELW corrections, see 
Appendix 

The mixed corrections evaluated by means of low-energy theorems could be cast into the expressions 


161) (162) 163}: 


GrM? (3 a;(M, 
Omixed = — 3 iat G zg c) aly for light quarks , 
GrM? a,(M, 
Omixed = eg 4 (1+ C2) ay for b and t, (103) 


at next-to-next-to-leading-order (NNLO) with ¢; = 17/6. Also available are the full mixed QCD- 
AN corrections 6° for H — bb which amount to about —0.08% for My = 125.09 GeV 


[164] [141]. 


In the left column of Fig. |1} taking Ira = 1 and ae = 0, we show the decay widths of a Higgs 
boson H into a pair of b quarks, c quarks, tau leptons, and muons for varying My. For the decays 
H — bb and H —> c@, the lower dashed lines are for the decay widths at leading order (LO) while the 
upper (black) solid lines are for those taking full account of the QCD and electroweak (ELW) corrections. 
The decay widths including only the electroweak corrections are denoted by the lower (red) solid lines. 
For the decays H > t*r~ and H —> pt", the dashed lines are for the decay widths at LO and the 
sold lines are for the decay widths including the electroweak corrections. The behavior does not alter 
much for other choices of (eee) as far as Baal + ral = 1 since the QCD correction gcp, 
which is common in the scalar and pseudoscalar contributions to the Higgs decay width into quarks, 
dominates. 


mixed 


In each frame of the right column of Fig. |1| the corresponding full decay widths are shown in the 
low mass region of 120GeV < My < 130GeV for the three choices of (iF oe) = (1,0), (0,1), 
and (1/v2,1/v2) denoted by the solid, dashed and dotted lines, respectively. The pure scalar case 


with Can T = (1,0) has a slightly smaller width compared to the pure pseudoscalar case with 
(97, Fp g} 7 1) = (0,1) due to the kinematical suppression factor of Be. 
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Figure 1: (Left) Decay widths of a neutral Higgs boson with mass My into bb, cé, T+7~, and pty 
from top to bottom taking dar = 1 and Tire = 0. See the main text for details. (Right) In each row, 


the corresponding full decay widths are shown for the three choices of (gf, Fp i 7p) = (1,0) (solid), 


(0,1) (dashed), and (1/2, 1/,/2) (dotted) in the low mass region around My = 125 GeV. The vertical 
lines locate the positions of My = 125.5 GeV. Note that the muonic decay widths are in units of keV. 
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At leading order (LO), taking the consideration of double off-shell effects, the decay width of a 
Higgs boson into a top-quark pair tt, each of which subsequently decays into bW* and bW~, is given 
by 
t mi; (Mu) g*Mu 


T™°(H > t(p;:) t(D) > BWW TW) = Né oe F, (104) 


where the dimensionless quantity F is given by an integrated function of the pole masses of b and t 
quarks, the W-boson mass, and the top-quark total width T; as 


2 =2 252 
= 1/2 1/2; 1/2 2} S |2 2) P |2 Pi + Pe | Ded; 
F = f~i Ar Àz [làs] EA + |Ap| rA | (1+ 2M2, + BF) 
pt+Me-My P+M- Mw 
(p? — M?)? + M?T? (pe — M?)? + M?T? 


dp; dp; , (105) 


with the 5 dimensionless triangle functions 


_ (i)? (Pi)? 


2 =2 22 

Pr Pi PePt 

A = 1 T T = 2 = 2 = 2 
2 M4 


Mā Mā Mā Mọ’ 


Àg = 


[a+ ay? ot- pB] - [a -NPP - t- ema] 2, 


is = 1+ A12) — (p2 — 92)? /M44 12 1— A22 — (p2 — 5?)2 /M4 1/2 2 106 
p= Utir] — We — Br) Ma FUAS = — Be) / Mi /2. (406) 
After integrating over p? and p?, we can recast the LO decay width into 

POH +t? > BW WT) = ghel T39 + onal TEO, (107) 


and, taking account of the radiative corrections as well as double off-shell effects, the Higgs decay width 
into a top quark pair r (H — t*t*) has been estimated as follows 


(H > tP) = |gia| TZO (1+ ögon + 5% + Shixea) (1+ Siw) + oial TEO (1 + daco + ôt”) (108) 


When My > 2M,, taking p? = p? = M? leads to Às = ys = B, = (1—4M?/M?,)'/? and àp = 1 and, 
neglecting the kinematical b-quark mass in the t + bW process, we reach the following factorized form 
for the LO decay width 


a. TLO(t — BW) \? 
LO Fur- _ PLO 
TMH + tt > BWt*W) | og, = OCH > tH) (7) , (109) 
using the narrow-width approximation (NWA) denoted by 
T 1 
S(p? — m?) = lim Z (110) 


r>0 m (p? — m?) + mT?’ 


l4For tw, we apply the SM approximation Eq. (102) both below and above the top-quark-pair threshold. Note that 
the QCD corrections are not valid in the threshold region due to the top-quark mass effects. For them, we refer to [87] 


and references there in. 
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Figure 2: Decay widths of a neutral Higgs boson with mass My into t*t* or [(H — t*t*) taking 
(Gers Gun) = (1,0) (upper) and (934, 947,) = (0,1) (lower). In the left panels, the vertical lines locate 
the top-quark-pair thresholds and the (magenta) dash-dotted lines are for the corresponding full 2-body 
decay widths of T (H — tt). Note that, in the right panels, we switch from [(H — t*t*) to T(H > tt) 
from My = 500 GeV and above. 


and the LO decay widths for H — tt and t > bW given by 


m- (My biMy 2 2 
pioi = a PM Toe sal? a] (111) 
2M, M2 2 M2 
Poe pw) = 22°t(1- 7) (o4 112 
eam = FA - Fe) (+a), (112) 


respectively. 


In Fig. 2 we show the decay width [(H — t*f*) as a function of My, taking (977,947) = (1,0) 
(upper) and (93, gin) = (0,1) (lower), respectively. We take T, = IT O(t — bW) assuming that a 
top quark decays 100% into a b quark and a W boson and T}O(H — tt + bbW+W) converges to 
[4°(H — tt) in the high My limit. Practically, far above the top-quark-pair threshold with My > 500 
GeV, we return to Eq. to suppress the contributions to [(H — t*t*) from the kinetic edge region 
of Vo + VE ~ My, assuming that the intermediate top quarks are reconstructed by requiring on-shell 
conditions of p? ~ M? and p? ~ M?. 


15For example, when My = 1000 GeV, we find that [(H — t*t*)/T(H — tt) takes the values of 1.07 and 1.10 for 
Cpe = (1, 0) and (957) 9 ra) = (0, 1), respectively. 
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Figure 3: Th°(H > FP) / Gee NÉ) as functions of My for the three choices of T |, Cae |) = (1,0) 
(black solid), (l9 pls Meryl) = (1/\/2,1/,/2) (blue dashed), and (19? pls Cal, = (0,1) (red dotted). 
We have taken my = my = 0,100,300 GeV from left to right. Note that there is no dependence on 


the choice of the couplings when mp = my = 0 as far as Ges |? + Seay |? remains the same. 


Finally, for decays into two different fermions such as charginos or neutralinos or for flavor changing 
decays such as H —> bs, one may write the effective interaction as 


Lary = — 979 HF (Gig tiap) + hic. , (113) 
without loss of generality. Then, at LO, the decay width may take a form of 


1/2 
Mir p 
T 


2 
pl TF fr ff’ 1 7 
TCH > ff) = NE EEE i-re) lH o pP) VEIS p- lg 


8 
where s = m4/Mh, K = m4 /Mh and Asp = (1 — & — K’)? — Aen’. The color factor NÉ = 3 for 
quarks and 1 for leptons, charginos, and neutralinos. For the decays into neutralinos E and X}, we 
need to multiply the factor of 4/(1 + jg) with jx = 1 for an identical Majorana neutralino pair. See 
Fig. {| for the My dependence of the normalized LO decay width TO (H > FFG p NE) for various 
choices of the couplings Com and Le eye and the fermion masses with mp = my. When my 4 mp, 
we observe that the LO decay width locates between the pseudoscalar (red dotted) and scalar (black 
solid) cases if the position of the mass threshold ms + my, is the same. In this review, we consider the 
Higgs decays into SUSY particles only at LO. For the QCD and ELW corrections to them, we refer to 


Refs. [168] [169] £70, [171] [172]. 


I?) , (114) 


Hf f! Hf f! Hf f! Hf f! 


3.2  Decays into two massive vector bosons: H > VV with V = Z,W 


Taking the full consideration of double off-shell effects, the Higgs decay width into two massive vector 


bosons is given by 165} {166} 


POH > V*V*) = o3 m evdy reer evda 
© Jo (y—1)? +e Jo (o- 1)? +e, 
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Figure 4: Comparisons of ['(H — V*V*) (solid) and T1O(H — VV*) (dashed). The vertical lines 
correspond to My = 125.5 GeV and g,,,,,, = 1 is taken. 


dvg2.... Gr M3 
VI avy TPH y1/2 GZ) È (1 = ; 4 ) f 22 ; (115) 
16V2n Wy Wy Wy Wy 


wy 
where dw = 2, ôz = 1, wy = 1/Ky = M},/M?2, ey =Ty/My, and X(a, b,c) = (a — b — c)? — 4bc. When 
My < My < 2My as in the case of the 125 GeV Higgs boson, the off-shell effects of one of the two 
vector bosons are negligible and the decay width reads 


rn > VV*) = 


Ovv« i eydx 
T Jo (g=1)? +e, 
VIivvGrMy 1 
16/2 = we 
with yy» = 2. See Fig. [4] for comparisons of ['(H > V*V*) and ['(H > VV*). Incidentally, one 


may neglect all the off-shell effects for a heavy Higgs boson with My > 2My and the decay width takes 
a simple form: 


Mwy, z, 1) [Mwy a, 1) + 12a] , (116) 


vg, GrMy 
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T(H = VV) = By [1 —4ky + 12K] , (117) 


where By = v1 — 4ky with ky = M?/M},. 

Beyond the leading order, including radiative corrections and the double off-shell effects above as 
well as below the gauge-boson-pair thresholds, we estimate the radiatively—corrected decay width into 
two vector bosons by introducing a correction factor as 


T(H + VV > 4f) =T"9(H > V*V*) (1+ bbw) - (118) 


elw 


For the SM electroweak correction factors 6%=“", we use the most recent version of PROPHECY4Fv3.0 


elw 


[91] (176) 177, 178, 179, (180) [81] to calculate the complete O(a) electroweak corrections to the Higgs 


decays into four fermions through intermediate W and Z bosons, supplemented by the corrections orig- 
inating from heavy-Higgs effects and final-state radiation. We note that these electroweak corrections 
are applicable even near and below the gauge-boson-pair thresholds where the narrow-width approxi- 
mation (NWA) is not applicable. The corrections amount to about 3% or less for My = 125 GeV, as 
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Figure 5: The electroweak correction factors 6” obtained from the complete electroweak correc- 
tions of O(a) to the decay processes of H > ZZ > e~et ppt (V = Z) and H > WW > vet WD, 
(V =W). The dashed lines are obtained by combining the results presented in Figs. 7 and 8 of Ref. 
while the solid lines by using the most recent version of PROPHECY4Fv3.0 taking the SM (or “model=0” ) 
in the “inputfile’. For our numerical analyses, we adopt the solid lines by PROPHECY4Fv3.0. 


can be checked in Fig. |5} For the consistent implementation of the ELW corrections using PROPHECY4F, 
we note that one may need to adopt the complex pole masses of the W and Z bosons for the LO decay 
widths which are given by 


My — il 
Mer _ ao — _Mv— uy (119) 


v1+I3/M} 


Using the complex pole masses, we find that the LO decay width into W*W* (Z*Z*) increases by the 
amount of about 0.4% (0.5%) taking My = 125.5 GeV, compared to those obtained using the on-shell 
masses and widths. 


It has been estimated that missing corrections beyond O(a) make the theoretical calculations for 
the inclusive decay rates into four fermions uncertain by the amount of 0.5% [182] [83]. 


In Fig. (6 we show the decay widths of a neutral Higgs boson with mass My into W*W* (upper) 
and into Z*Z* (lower). Note that we include the electroweak corrections shown by the solid lines in 
Fig. |5| We observe that the electroweak corrections grow as My increases and they make the partial 
decay widths into WW and ZZ as large as 350 GeV and 180 GeV, respectively, around My = 1 TeV. 
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Figure 6: Decay widths of a neutral Higgs boson with mass My into W*W* (upper) and into Z*Z* 
(lower) or [(H + VV — 4f) taking g,,,,, = 1. In the left panels, the vertical lines locate the vector- 
boson-pair thresholds and the (magenta) dash-dotted lines are for the corresponding 2-body decay 
widths of (H > VV) = TO (H > VV)(1+4+ ô% ). 


elw 
3.3 Decays into a lighter scalar boson and a vector boson and into two 
lighter scalar bosons: H > pV ,yy 


In the presence of multiple Higgs bosons, the decay of a heavier Higgs boson H into a lighter neutral 
Higgs boson y and a massive gauge boson Z may occur and, considering the case of a virtual Z*, its 


decay width is given by an integral form as [173} [166] [174] 
GoM fay aha 
8/20 A w3r| (x — 1)? + 3] 


with w = Mj,/MzZ and wp = M3/Mz. When My is larger than M, + Mz, using the Z-boson narrow- 
width approximation, it reduces to 


(nS ọZ*) = (120) 


Gr Mẹ} » 
82r HpZ 


where ky = M/M} and kz = Mz/Mj. And, a heavier Higgs boson H might also decay into a lighter 
charged Higgs boson y* and a massive gauge boson WF with its decay width given by an integral form 


as [173] [166] [74 


[oR > yZ) = N/?(1, Ky, Kz), (121) 


7 GrMẹlI atw- |? a ewA3/2(w, w4, £) 


TLO(H > gtW™) = a l an ae (122) 


32 


where M°(H > ptWF*) = TYR > gtW*) = M(H > gp W*) with w = M?/M? and 
w+ =M p /Mj,. When My is larger than M,+ + Mw, in the W-boson narrow-width approximation, 
it reduces to 


Gp M3 
ee cones 


T°(H > y*WT) = AA ko, kW), (123) 


where ky = M2./Mj, and kw = Mj,/Mj,. 
Finally, when a heavier Higgs boson decays into a pair of lighter neutral Higgs bosons of p; and gj 

or into a pair of sfermions, at LO, we have for the decay width 
v"|G/? 

“le 167M y 


pS (H > poni fi) = = AM?(L, ki, K3), (124) 
where (No, G) = (1 + fij, Gueio,) OF (No, G77) and k; = Meg | Mir- The decay width into a pair of 
lighter charged Higgs bosons is given by taking (Nọ, G) = (1,9,,,,,,,-)- Again, the Higgs decays into 
SUSY particles are considered only at LO. In H —> yj; decay, we neglect the off-shell effects assuming 
Ty ~ Tusu with a caution that the off-shell effects can not be neglected when the decay involves a 
transition between two heavy states with wide widths. 


In Fig. [ we show the LO decay width of [¥°(H — ~~ W**) (upper) which is the same as TO (H > 
ytW-*), TO(H — ~Z*) (middle), and T™O(H —> yip;) (lower). For the mass of a lighter scalar boson 
p, we take three values of 50 GeV, 200 GeV, and 300 GeV for H > yV* for illustration. For the decays 
H > ip;, we take M,, = My, = 50 GeV, 100 GeV, 200 GeV, and 300 GeV. All the relevant couplings 
are taken to be 1 in the numerical analyses. 


3.4 Decays into two gluons: H — gg 


By introducing two form factors, without loss of generality, the amplitude for the decay process H —> gg 
can be written as 


_ as(Mp) M3 om 
Ar v 


My = (An) ia 6.) — P Ma) spe eieskoha)} (125) 
where a and b (a,b = 1 to 8) are indices of the eight generators in the SU(3) adjoint representation, 
ky the four momenta of the two gluons and €12 the wave vectors of the corresponding gluons, ef, = 
ef — 2kf (ko - €1)/MF, 5, = & — 2k (ky - €&2)/Mĝ and (eenkike) = €wpo pegki{k. Retaining only the 
dominant contributions from third-generation quarks and introducing AS and APY to parameterize 
contributions from the triangle loops in which non-SM colored particles are running, the scalar and 
pseudoscalar form factors are given by [£] 


S( Mu) = X. gape Fep(ts) +A; PM) = X ggr Fos (ts) + APS, (126) 
f=b,t,c f=b,t,c 


where te = Mj, /4M? is defined by using the pole masses of the bottom and top quarks. The form 
factors F’,, and Fps can be expressed by 


F = r a= F=: (127) 
in terms of a so-called scaling function f(T) which stands for the integral function 
f(T) = of dy In [1 — 4ry(1 — y)] = 1 ee | 72 Pere (128) 
2/0 y -4 [mn (E) - an] ; Tade 


16See Appendix [B] for AS and AP9 in the MSSM. 
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Figure 7: The LO decay widths of a neutral Higgs boson with mass My into p-W** (upper), yZ* 
(middle), and p;p; (lower) taking aipe = l, 94,2 = 1, and N,|G|? = 1, respectively. In the 
upper and middle panels, the (magenta) dash-dotted lines are for the corresponding 2-body decay 
widths and we are taking M,/GeV = 50,200, and 300 from left to right. While in the lower panel, 
M,/GeV = 50,100, 200, and 300 are taken. 


It is clear that the imaginary parts of the form factors appear for Higgs-boson masses greater than twice 
the mass of the colored particle running in the loop, i.e., 7 > 1. In the limit 7 > 0, F,,(0) = 2/3 and 
F,(0) = 1, see Fig. 
The decay width of the process H —> gg may be cast into the form 
Mj,0% 2 “8 S 2 :P :P 

M(H > g9) = 9% [I9 (1+ éo) (1+ 688) + P Ma) (1+ 686p) (1+ 682) | (129) 
including the QCD and electroweak corrections. The two corrections are factorized as dictated by the 
universal infrared and collinear behavior of QCD corrections and the universality of the dominant part 
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Figure 8: Behavior of the real (solid) and imaginary (dashed) parts of the four form factors of Fsp, 
F,,-, Fı and Fo versus T = s/4m? with ys = My and m being the mass of the particle running in the 
triangle loops. The form factors are defined explicitly in the main text where they first appear, see 


Eqs. (127), (136), and (B.3). The vertical lines denote the mass threshold above which My > 2m. 


of ELW ae os The QCD correction of bgcp is known up to NLO including the full quark mass 
dependence [183] and up to the NLO in the limit of heavy top quarks 185} [186]: 


ae 95 7 as\ oS?) (Ma) 
docp = a -gNr + An — 


+ 


2 
M3 (Np) M 
[37020 — 47.19N# + 0.902(NZ)? + (2.375 + 0.667 N£) log =] ( (Mu) 


+ 


1583.10 — 1062.82N% + 52.62(NF)* — 0.5378( NZ)? 


M3 
+(66.66 + 14.60 NE — 0.6887(N£)°) log — 


M? 
3 
M2, (Ne) M 
+(6.53 + 1.44NE — 0.111(N£)?) log? = ( : n) ; (130) 


with NE = 5 counting the flavor number of light quarks and A%° ~ 0.7 F] for the NLO quark mass 
effects from the top, bottom and charm quarks [183]. Taking My = 125.5 GeV, we find 5 = = 
0.67 +0.20 + 0.02 for the NLO, NNLO, and N3LO corrections. On the other hand, the QCD corrections 
to the pseudoscalar part develop a Coulomb singularity at the top-quark-pair threshold which can be 


1TWe take the value from Fig. 7(b) in Ref. [183]. 
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Figure 9: Behavior of the two-loop SM electroweak corrections of Oe for [(H — gg) versus My. The 


vertical lines locate My = 125.5 GeV and the WW, ZZ, and tt thresholds. The line width represents 
uncertainties taking account of M, = 172.64+1.58 GeV. From the authors of Refs. [198} [197] [199} [200]. 


See also Refs. [201] [196] [194]. 


regularized by including the finite top decay width [187], [88] [189] [90]. The correction SLE Qcp is known 
up to NNLO in the limit of heavy top quarks [191] while the NLO corrections are known oil 133}: 


ge, = (T_T ye age os” (Mn) 
aom = (I eee o 


2 
Mi, WEN My) 
+ (392.22 — 48.58NE + 0.888(N4)? + NE log i) í = , (131) 


with AFP œ~ 0.04 for tan = 1, for example, in the MSSM. || a Taking My = 125.5 GeV, we find 
ae ~ 0.66 + 0.22 for the NLO a NNLO terms. We note that doosi is derived in the heavy top-quark 
limit and it does not contain the Coulomb singularity. 


The electroweak corrections 5%° of O(Gr M2) are given by [163] 192 


2 
oo (132) 
8/27? 
Note that these O(GpM?) electroweak corrections increase the gluonic decay width only by the amount 
of about 0.3%. It turns out that the full electroweak corrections lead to much 
more significant enhancement of Behe, ~ 5% below the WW threshold, see Fig. lo] On the other hand, 
the other electroweak corrections Ta may be given by 


Gr M? 
gP ay Fitt 
Ona J - (7+ co 8\/2r2 ’ (133) 


18In the MSSM, we note that A9:? significantly depends on tan £, see Fig. 20(b) in Ref. [183]. 
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Figure 10: (Upper) Normalized decay widths of a neutral Higgs boson with mass My into gg taking 
Tone! Grit = 1, AS? = 0, and Ita = AP’ = 0. In the right panel, the low My region is magnified. The 
LO (blue dashed), LO+ELW (lower red solid), and LO+ELW+QCD (upper black solid) results are 
separately shown. (Lower) The same as in the upper panels but for three values of I (oe = 1,10, 
and 20 taking full account of the electroweak and QCD corrections. The vertical lines in the right panels 
locate the position My = Mi. 


where the first factor of 7 counts the contribution of the SU(2); doublet which, in the decoupling limit, 
plays the role of the SM SU(2); doublet including the 125 GeV Higgs boson. And the second factor 
of nae denotes the model-dependent BSM contribution. In the type-I] 2HDM, for example, it is given 
by i. = 10/ tan? 8 in the infinite top quark mass approximation [203]. These electroweak corrections 
reduce the gluonic decay width at the percent level. Note that Eq. is a crude approximation 
obtained by focusing on the top-Yukawa contributions which are subleading in the SM. The same 


arguments apply for Eq. (132) which we are not using for our numerical study though. 


We have addressed all the known QCD and electroweak corrections to the decay width of a neutral 
Higgs boson into two gluons. The theoretical uncertainties due to the unknown higher-order QCD and 
NLO electroweak corrections are estimated as 3% and 1%, respectively [182] [84]. 


In the upper panels of Fig.|10| we show the normalized decay widths '(H — gg)/ (9%)? at LO (blue 
dashed), including only the electroweak corrections (red solid), and including both the electroweak and 
QCD corrections (black solid). We assume that all the pseudo-scalar couplings of Iiiz are vanishing 


and AS¥ = 0. The electroweak corrections o are directly read off from Fig. lo In the right panel, we 


magnify the low My region and locate the position My = M, with a thin vertical line. In the lower 
panels, we show the normalized decay widths including the available electroweak and QCD corrections 
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Figure 11: The same as in Fig. [10] while taking GFtaq = 0 and Tia #0 with AS? = AP? = 0. 


and taking three values of g>,,/9>;, = 1 (black solid), 10 (blue dashed), and 20 (red dotted), considering 
the situation in which the bottom Yukawa coupling is enhanced as in the 2HDM II for large tan 6. Note 
that the black solid lines in the lower panels are the same as those in the upper ones. 


In Fig. the alternative choice is made to show [(H — gqg)/(gj;;,)? assuming all the scalar 
couplings of Gita are vanishing and APY = 0. For the electroweak corrections oe , we take nin, =), 
Compared to the scalar case in which the form factor fsf(T) is involved, the rise near the top-quark-pair 
threshold is sharper and bigger due to the behavior of the real and imaginary parts of the form factor 
fos (7) around 7 = 1, see Fig.|8| In the low My region, we find that the decay widths are larger by about 


the factor of [f,-(0)/fs(0)]? ~ 2 when the b-quark contributions are neglected. As the coupling ratios 


a, / A increase, the contributions from b-quark loops become comparable to and larger than those 
from t-quark loops and, in this case, the decay widths are nearly the same especially around My = 100 


GeV as can be checked with the dashed and dotted lines in the lower-right panels of Figs. [10] and [L1] 


3.5 Decays into t 


a(0)M? S De aie de 
Mng = -ZOME Fst) (cha 81) - PM) ay (etki) } 
A 
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Atv an 


in terms of the two form factors of S% and P”. Here ki. and €12 are the four-momenta and wave 
vectors of the two photons, respectively, as in the decay H — gg. Note that the electromagnetic fine 
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Figure 12: Behavior of the real (solid) and imaginary (dashed) parts of the scaling factors Css(T) (left) 
and Cps(T) (right) for the QCD corrections to the decay of a Higgs boson into two photons. 


structure constant in the coupling should be taken at the scale q? = 0 since the final-state photons are 
real. Retaining only the dominant contributions from third-generation fermions and the charged gauge 
bosons W= and introducing two residual form factors AS’ and AP” to parameterize contributions 
from the triangle loops in which non-SM charged particles are running, the scalar and pseudoscalar 
form factors are given by [E] 


S'(Mz) = 2 X. NEQ oR pp Festli) — GuwwFiltw) + AS"; 
f=b,t,7 

P'(My) = 2 5°) NE Q} ofp, Fos(ts) + AP’, (135) 
f=b,t,7 


where NZ = 3 for quarks and NZ = 1 for charged leptons, respectively. For the 7 lepton and the 
W boson, T, = M2 /4M2 and tw = Mj,/4Mj,, respectively. On the other hand, for quarks, T4 is 
defined in terms of the running quark mass at the scale of My/2, i.e. Tg = M},/4[m (Mu/2)|*? where 
Mz, is normalized as m,(M,) = M}. Note that the choice of the scale u = My/2 correctly gives 
Mq(My/2) = M; at the threshold where My = 2M, and makes the full two-loop QCD corrections 
remain small in the entire range of the variable 7, by effectively absorbing all relevant large logarithms 
into the running mass. The form factor F; is given by 


F(t) =2+87* +37 (Q—7 “)f (7), (136) 
which takes the value of 7 in the limit r — 0, see Fig. |8| for the r dependence of the form factor. At 
LO, the decay width of the radiative process is given by 

POC > yy) = HS [> Mu) + PUMP (137) 
25673 v2 = Had 
with the fine structure constant a = a(0) ~ 1/187. 


As a gluon cannot be radiated from the colored quark loop contributing to the Hyy vertex owing to 
charge conjugation invariance and color conservation, the full two-loop QCD corrections [183} [204] 205, 


19See Appendix ne AS? and AP” in the MSSM. 
20For mq(u), see A ppendix |A} 
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206}, 207, [208] [209] (210) [211] (212) [213] (214) [215] [216] can be taken into account by simply introducing 


the scaling factors in the form factors for the b- and t-quark contributions as: 


a) | 


Pala) — Fln) [i+ Castra 


a) . (138) 


Bo —? Foplta) [i+ Cpr) 


The scaling factors C,¢ and Cpr approach —1 and 0, respectively, in the limit 7 — 0, see Fig. 


The asymptotic values of the scalar and pseudoscalar scaling factors in the heavy quark limit, 
Cy ,(0) = —1 and C,,(0) = 0, can also be deduced by means of a general low-energy theorem for ampli- 
tudes involving soft Higgs particles. According to the low-energy theorem, the NLO QCD corrections 
to the scalar coupling to two photons in the heavy quark limit can be obtained from the effective 


Lagrangian 183] 218] 


S Q — BQep/O 


B H 
a A 41 14+ m(as) 


E FY Ew, (139) 


where the QED £ function 64pp/a = 2(a/7)[1 + as/m + ---] with the heavy quark q contribution and 
the anomalous mass dimension Ym = 2a,/7 +--+ [220]. Therefore, the NLO expansion of the effective 
Lagrangian of the scalar coupling to two photons reads 


Le G72, ee Oe E - “4 O(a?) (140) 
127 T v 

which agrees with the asymptotic value of C,;(0) = —1 in the heavy quark limit. On the other hand, 
the pseudoscalar scaling factor is zero in the heavy quark limit, i.e., Cpf(0) = 0, because the QCD 
corrections to the pseudoscalar decay mode vanish in this limit due to the Adler-Bardeen theorem [221]. 
The effective Lagrangian for the pseudoscalar coupling to two photons can be derived from the Adler- 
Bell-Jackiw (ABJ) anomaly in the divergence of the axial vector current and is given to all 
orders of the perturbation theory by 


= 
Lig = Gite Qe oP Fy (141) 


for the pseudoscalar boson, as the lowest order axial-anomaly term is not renormalized at all orders. 


The two-loop electroweak corrections of ie to the scalar part of the decay width via the scalar 
form factor have been calculated in Refs. [198] [224] [201] 200|, see Fig. The electroweak corrections 


to the pseudoscalar part are also available |202! [203]: 


GrM? 
Na t (4 + at 142 
elw 82x2 Mw ( ) 


where, similarly as in the aroa a — gg, the factor 4 counts the contribution of the SM SU(2)z doublet 
in the cco limit and nie the model-dependent BSM contribution. In the type-I] 2HDM, for 
example, me = = 7/ tan? 8 in ita — top quark mass approximation [202| [203], suppressed for large 
tan 6. Note again that Eq. ( is a crude approximation obtained by focusing on the top- Yukawa 
contributions which are va in the SM. 


?1For a detailed description of the scaling factors of C, (7) and Cy ¢(T), see Appendix |C] 
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Figure 13: The two-loop electroweak corrections vii to the part of the decay width of a Higgs boson 


into two photons via the scalar form factor in the full complex-mass scheme which has been designed 
to cure the unphysical infinities of two-loop amplitudes [198]. The line width represents uncertainties 
taking account of M, = 172.64 + 1.58 GeV. From the authors of Refs. [£98] [197] [199] 200]. 
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Figure 14: Behavior of the two-loop QCD corrections to the decay width of a Higgs boson into two 
photons: dQGp (solid) and dap (dashed). We take guww = Jfr = Gon) = ISe = 1 for ôo while 
= gł, = 0 for ion: Both of the residual form factors AS? and AP’ are taken to be vanishing. 
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Incorporating all the QCD and electroweak corrections, one may write 


3 


M(H 4-97) = Lo oA [1a (1+ 8885) (1 + 0X) + 1PM) P (1+ 8860) (1+ 8) | 243) 


Al 


Note that the electroweak corrections are directly from Fig. [13]and Eq. (142) while the QCD corrections 
enter through the scaling factors C',¢(7,) and Cr(7,). For My = 125.5 GeV, the QCD corrections 


n increase the decay width into two photons by about 2%. On the other hand, the electroweak 


S 


corrections 07” decrease the decay width by about 2% , almost canceling the NLO QCD corrections 


to the corresponding part. In Fig. we show the QCD corrections ton (solid) and Slan (dashed) 
with varying My. For the scalar QCD correction Tm we take the SM values of gpww = Jon = 
UFAN = g, = 1 with AS = 0. While, for the pseudoscalar QCD correction Aleni we assume a 
scenario in which the pseudoscalar form factor PY is dominated by the top-quark contribution taking 
Jit = gin = AP = pil At My = 2M,, the pseudoscalar QCD correction aa diverges due to 
the singular property of Cpr at 7 = 1. Around My = 600 GeV where the large cancellation occurs 
between the W-boson and top-quark contributions, the scalar QCD correction Doers is relatively large 
and it could vary between about —0.4 and 0.8. 


In the upper panels of Fig. we show the decay width of a neutral Higgs boson into two photons 
normalized to the gym coupling squared taking gpww/9fp = Uae! Cae = 1 and Ge = AS = 
AP’ = 0. This reduces to the SM decay width of a Higgs particle with mass My when gy = 1. 
For My < 170 GeV, we apply the electroweak corrections ie directly read off from Fig. and we 
simply neglect them above My = 170 GeV, expecting the electroweak corrections to be evaluated with 
high precision for heavy Higgs bosons if necessary. Below the W-boson-pair threshold, the W-loop 
contributions are dominant, leading to the sharp rise as My approaches 2My. Passing My = 2My 
from below, the real part of the W-loop contributions decreases but its imaginary part starts to be 
developed. As a result, the Higgs decay width continues to increase with increasing My until the 
Higgs mass meets the top-quark-pair threshold, My = 2M,. Passing the top-quark-pair threshold, the 
real and newly-developed imaginary parts of the t-quark loop contributions cancel those of the W-loop 
contributions [210]. We observe that, beyond My = 2M, the cancellation between the real and 
imaginary parts of the W-loop and t-loop contributions broadly occurs and it leads to a dip around 
My, = 600 GeV. Specifically we find 


S7(My = 600 GeV) © (=1.593 — 2.5994) Guy + (1.391 + 2.131 i) gS, (144) 


with negligible contributions from the b-quark and 7-lepton loops. In the middle panels of Fig. we 
show the variation depending on g}5,/9h = Gitrr/ Gre Still taking gyww/gfr = 1- The (black) solid 
lines represent the same case as in the upper panels taking full account of the QCD and electroweak 
corrections. And, in the lower panels, we show the results taking g,,,,y, = 0. The last case may apply 
to the heavy neutral Higgs bosons appearing in the 2HDMs and/or MSSM when their couplings to the 
massive vector bosons are naturally suppressed and almost vanishing [225] [226]. 


In Fig. the alternative choice is made to show the normalized decay width ['(H > yy)/(giu)? 
assuming that all the scalar couplings of Ge Ff are vanishing and, again, taking AS’ = AP? = 0. Note 
that, in this case, only the fermion loops are contributing. For the electroweak corrections OF, we 
take mie = 0. At My = 2M;, the decay width diverges because of the singular property of the QCD 
corrections. In the lower panels, we show the dependence on Iii Gin = GHre/ Gy, for the four values 
of 1,10,20, and 50 taking full account of the electroweak and QCD corrections. In the right panels, as 
the same as in Fig. |15| we magnify the low My regions. 
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Figure 15: (Upper) Normalized decay widths of a neutral Higgs boson with a mass My into yy taking 
Cie = Deri ure! TH = ] and Gee = 0. In the right panel, the low My region is magnified. 
The LO (blue dashed), LO+ELW (red solid), and LO+ELW+QCD (black solid) results are separately 
shown. (Middle) The same as in the upper panels but for the four values of g%4,/9%pa = 1,10, 20, 
and 50 taking full account of the electroweak and QCD corrections. gyww/9%g = 1 is taken. (Lower) 
The same as in the middle panels but with g „ww = 0. The vertical lines in the right panels locate the 
position My = 170 GeV beyond which the electroweak corrections are ignored. In all panels, we take 
AS = AP” = 0. 
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Figure 16: Normalized decay widths of a neutral Higgs boson with a mass My into yy via the 
pseudoscalar form factor P7 with AP? = 0. In this case, only the fermion couplings of g} if with 
f =t,6,7 are relevant. 


3.6 Decays into a vector boson Z and a photon: H > Zy 
The amplitude for the decay process H — Zy can be written as [227] 


MzyH = — {87 Mz) [ki - kze- €} — ki- & ke - ef] — P7 (Mp, Mz) (chit) } , (145) 


with the two form factors of $77 and P77. || Here kı 2 are the momenta of the Z boson and the photon 
(we note that 2k, - ko = M?, — M2), €12 are their polarization vectors, and (e:¢9kik2) = Ewapenegkoks. 
Retaining only the dominant contributions from third-generation fermions and W~ and introducing 
two residual form factors AS” and AP% to parameterize contributions from the triangle loops in 
which non-SM particles are running, the scalar and pseudoscalar form factors are given by 


ie = 9520 1 
S?(My,Mz) = 2 X NEQe2——E* gj Galp As) — — Gn Galtw, Aw) + AS, 
f=t.b,r Ww wW 
Zy f If — 25%Qs P Zy 
P”(Myg,Mz) = 2 >` NASE Iarr GorlTp As) + AP”, (146) 
f=t,b,7 


2Tn the general case with arbitrary gyww and Grr couplings, the scaling factors C's(Tq) and Cp f(7q) should be taken 
into account at the amplitude level to incorporate the corresponding QCD corrections properly. 

?3We follow the conventions and notations of Ref. [94] for the form factors. 

24See Appendix [B] for the explicit forms of AS77 and AP” in the MSSM. 
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Figure 17: Behavior of the real (solid) and imaginary (dashed) parts of the form factors G'ss(Te, Az) 
and Gpf(T:, à+) (upper) and that of Gi(7w,Aw) (lower). We recall the relations 7, = Mj7,/4M? and 
A» = M?/4M?. The vertical lines at 7, = 1 denote the mass threshold above which My > 2M,. 


where I” = +1/2, 1° = —1/2, T, = M3 /4m2 and A, = M2/4m2, respectively. |5| The loop functions 
are given by: 


Gep(te, As) = To(te, As) — Lite. As); Gor (te, Àf) = (Ts, Às), (147) 
Gi(tw, Aw) = Cw l [2rw (tiy = 1) + (ty = 5)] Ii (Tw, Aw) + 4(3 = 2) La(tw, aw) j 


where J; are functions of the two variables of 7 and À and they are expressed as 


AG) = geo t aR - FON giy O-A], 


bA = -pl ~ FO], (148) 


in terms of the f(r) function, defined in Eq. (128), and the function g(r) which is defined as 


+ — Jarcsin(/7) : Tl; 
_ V 149 
0) \ Eh (MED) -a 21 a 
aN = aaa i at 


25We take the pole masses of top and bottom quarks for Tib and Azo. 
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18: (Upper) Normalized decay widths of a neutral Higgs boson with a mass My into Zy at LO 
Gece 0,5 = 1 and ae = AS*7 = AP” = 0 for four values of g>5, p -/9frn = 1,10, 20, and 


50. In the right panel, the low My region is magnified. (Lower) The same as in the upper panels but 


taking 
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19: Normalized decay widths of a neutral Higgs boson with a mass My into Zy via the 


pseudoscalar form factor P77 at LO taking AP% = 0. In this case, only the fermion couplings of a if 
with f =t,b,7 are relevant. 


The explicit 7, dependence of the form factors Gss and Gps for x = t and Gi for x = W is shown in 
Fig. clearly exhibiting the development of their imaginary parts beyond My > 2M, with A,w < 1. 
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At LO, the radiative decay width is given by [228] [229] [230] 


a(0) GE May siv M3 


M(H > Zy) = aa 7 


M2\° 

( - 9 (|57 Ma, Mz)|? +|P2(My,Mz)|’) . (150) 
H 

A detailed description of the scalar and pseudoscalar form factors in the framework of the MSSM taken 

as a specific BSM model is given in Appendix |B] 


The QCD corrections turn out to be less than 0.3% [2311232233]. On the other hand, the theoretical 
uncertainties of the electroweak corrections have been estimated as ~ 5% which constitutes the 
largest theoretical uncertainty involved in the experimentally clean and/or dominant SM Higgs decay 
modes into bb, WW ,gg,77,c®,ZZ,yy,Zy7, and up. 


Experimentally, the H —> Zy mode has to be extracted from the Dalitz decays of H > yff 
6] which consist of: (i) the H + Z“)y decay followed by the decay 
Zœ) — ff, (ii) the tree-level process H — ff with a photon radiated from the final-state fermions, 
(iii) the loop-induced process H — 7y* via triangle diagrams followed by the decay y* > ff, (iv) the 
loop-induced process H — yf f via box diagrams and (v) the loop-induced process H > ff via triangle 
diagrams with a photon radiated from the final-state fermions. For a clean separation of the H > Zy 
decay from other processes, appropriate experimental cuts have to be imposed. 


In the upper panels of Fig. we show the LO decay width into a vector boson Z and a photon 
normalized to the gy# coupling squared taking gyww/9fn = 1 and D = AS” = AP” = 0 for the 
four values of Gerber! an = 1, 10, 20, and 50. This reduces to the SM decay width of a Higgs particle 
weighing My when gyywy/9on = Dir Hellan = 1 together with gym = 1. Below the W-boson- 
pair threshold, the W-loop contributions are dominant, leading to the sharp rise as My approaches 
2Mw. Passing the W-pair threshold My = 2My from below, the real part of the W-loop contributions 
decreases while the imaginary part starts to develop leading to another mild rise, see the lower panel 
of Fig. Passing My = 2M,, the t-quark loop contributions start to cancel that of the W-loop 
contributions as shown in the upper-left panel of Fig. In the lower panels, we show the results 
taking 9,yyw = 0 for the four values of Oe as Hor! Girt = 1,10,20,50. This case may apply to the heavy 
neutral Higgs bosons appearing in the 2HDMs and/or MSSM when their couplings to the massive vector 
bosons are suppressed and they are almost vanishing [225] [226]. 


In Fig. {19} the alternative pseudoscalar choice is made to show T(H —> Z7)/(girz)° taking the four 
values of gf5,/9in% = GHrr/ Gn = 1,10,20,50 with AS77 = AP%7 = 0. Here we assume all the scalar 
couplings of re zp are vanishing. Note that, in this case, only the fermion loops are contributing as the 
pseudoscalar state does not couple to gauge bosons at the tree level. In the right panel, as the same as 
in Fig. we magnify the low My regions. Compared to the scalar case shown in Fig. it is clearly 
shown in Fig. [19] that the scalar and pseudoscalar decays exhibit quite distinct patterns, in particular, 
around the t-pair threshold. The distinct patterns between the scalar case with g „ww = 0 (two lower 
frames of Fig. and the pseudoscalar one (Fig. basically come from the differences in sizes and 
behaviors of the form factors Gss and Gps shown in Fig. 


Before moving to the last subsection for numerical results obtained by analyzing the decays of several 
neutral Higgs bosons, we provide Table [4] in which the leading My dependence and the ballpark values 


?6Very recently, the ATLAS collaboration has reported on an evidence for the H — lly process for a Higgs boson with 
a mass of 125.09 GeV and a dilepton invariant mass mee < 30 GeV with £ = e and u [241]. The observed significance 
is 3.20 over the background-only hypothesis, compared to the expected significance of 2.10 for the SM prediction. We 
note that, for the low values of mgp, the Dalitz decay is dominated by the loop-induced H > y*y — €¢y process with 
subleading contributions from the loop-induced H —> Z*y — eey and tree-level H — upy processes for £ = e and p, 
respectively. 


AT 


Table 4: The leading My dependence and the ballpark values of normalized decay widths of the neutral 
Higgs boson H far above the mass threshold of its decay products. Specifically, we have taken My = 1 
TeV. For the radiative loop-induced decays, the scalar (upper) and pseudoscalar (lower) contributions 
are shown separately. For the scalar contributions to H — yy and H — Zy which are dominated 
by t-quark and W-boson loops, we assume 977, = Jyww = 1. See Table [12] for the MS quark masses 
Mg( My) while we refer to Eq. for malu = My/2). 


Decay Mode Leading My dependence [pprmaed [GeV] | Reference Figure 


H f M 40 [Zam] Figs. |1|and |2 
> if H [E 1gs. jl) an 


2 
A Mz p | 


3 | mt(Mz) 
H > WW (ZZ) M3, 350 (180) Fig. 6 
H >V M3 300 
H > pọ 1/Mp 4 
2 
H > gg Mir Ess (ai) TRI 
2 
Mè |For (2) 9 x 107? 
8 M2 M2, \|? 
Taa Mir Foy (mata) = (2) 3x104 
M2 2 7 
Mi [For (mnt) oe 


A M2 M2 M} M2 \(? 4 
mg [0.9 Gss (1 ee) — Gr (H ae 6 x 10 


H > Zy 


2 
3 Mi MZ -5 
Mi, |Gpy (H ane 8 x 10 


of normalized decay widths at My = 1 TeV are shown for all the decay modes elaborated on up to 
this subsection. Table |5} is further provided for a summary of the QCD and electroweak corrections 
considered in Section |3] and Section [4| which are for the decays of H and H™, respectively. 


3.7 Numerical results 


Closing this section dedicated to a detailed study of neutral Higgs boson decays, we present the results 
of two numerical analyses of (i) the decays of a neutral Higgs boson with its mass fixed to 125.5 GeV 
dictated by the LHC discovery and (ii) the decays of heavy neutral Higgs bosons which are mixtures 
of CP-even and CP-odd states. 


In the first numerical analysis, we extend the SM in a somewhat model-independent way by allowing 
for the pseudoscalar as well as scalar couplings of the 125.5 GeV Higgs boson and perform a compre- 
hensive analysis of its decays by estimating all the widths and branching ratios as precise as possible. 
For the second numerical analysis, we have specifically chosen the type-I 2HDM in which the Yukawa 
couplings of the lightest Higgs boson as well as its couplings to a pair of massive vector bosons quickly 
approach the corresponding SM values as the masses of the heavy neutral Higgs bosons increase and 
their decouplings are not delayed (244). In this model, there is no much need of decoupling 
the heavy Higgs bosons to avoid conflicts with the current LHC Higgs precision data. Moreover, though 
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Table 5: 


The QCD and ELW corrections to the decay widths of neutral and charged Higgs bosons 


which are denoted by H and H*, respectively. The abbreviations used are: NC = Not Considered, [S] 
= Scalar part, [P] = Pseudoscalar part. 


H > Tr, utu 


H >tt 


H>WW,ZZ 


Decay Mode Expressions for a partial width QCD corrections ELW corrections 
H > bb,cé Eq. (101) 


[S] : Eq. 


[P] : NC 


H > pZ 


H > yp*Wwt 


H > yy 


Fig. |5| (PROPHECY4Fv3. 0) 


NC 


H —> gg 


H= 


H > Zy 


Ht > cs 


At > vrt ‚vut 


Ht > #5 Eq. for My+ < 500 GeV 
Eq. for My+ > 500 GeV 
Ht > HWt Eq. f167) 


assuming that the lightest neutral Higgs boson is a purely CP-even state, the heavy neutral states 
could still undergo a significant CP-violating mixing in the presence of CP-violating phases in the Higgs 


potential [247]. 


We emphasize that the numerical analyses performed in this subsection are solely based on the 
detailed analytical and numerical results presented in this Section B]and supplemental materials provided 


in Appendices. 


3.7.1 Anatomy of Higgs boson decays with My = 125.5 GeV 


In this subsubsection, we present all the details of calculating the decay widths of a neutral Higgs boson 
by taking My = 125.5 GeV [10], while allowing for nontrivial pseudoscalar as well as scalar couplings 
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of the Higgs boson to fermion pairs. The SM parameters taken in this present analysis are summarized 
in Appendix |A] Note that one may apply the results to the genuine CP-even SM Higgs boson by taking 
are = Jayy = 1 and JHF = AS90427 = APII = 0). 


For the decays of a neutral Higgs boson with My = 125.5 GeV into b and c quarks and muons and 
tau leptons, the relevant radiative corrections are given numerically by 


dacp = +0.203 + 0.037 + 0.002 — 0.001 = +0.241 ; 
Aa = 001LE n G =o, 


Saw = —0.005, Siw = +0.004; 5%, = —0.0097, 54, = —0.0307 ; 
Osea = 0.001, Onivea = —0-000; 
a = 40.244, ô, =+0.265; 67, =—0.0097, 5%, = —0.0307, (151) 


with the rescaling factors defined by €7° = gs i ae Incidentally, we have 62? = 0.018 £"? and 
ôP = 0.030 ¿EP contributing to the pseudoscalar part. 


For the decays H — WW, ZZ, from Fig. |5| we have two electroweak corrections as 


ôW = 0.0308, 64, = 0.0152. (152) 


elw elw 


For the decay H — gg, the scalar and pseudoscalar form factors and the relevant radiative corrections 
are given numerically apart from the scalar and pseudoscalar Higgs—fermion—fermion couplings by 


S9 = 0.688 giz, + (—0.043 + 0.0622) g?g, + (—0.009 + 0.0087) g9 + ASI; 
PI = 1.048 gig, + (—0.049 + 0.063%) giz, + (—0.010 + 0.008 i) gira. + AP! ; 


ol iS :P oP 
ddcp = 0.8850, 6%° = 0.0516; ddcp = 0.8775, ôT- = —0.0218, (153) 


elw elw 


where 6%° is from Fig. [9] and we set n4? = 0 for 6%”. In the SM, we have S9 = 0.636 + 0.070% and 


elw elw’ 


PI = 0. Similarly for the decay H — yy, in terms of the HWW and Higgs-fermion-fermion couplings, 
we obtain the following scalar and pseudoscalar form factors 


SY = —8.341 gyywy + 1.826 g?r + (—0.020 + 0.024 i) gg, + (—0.023 + 0.021 i) g3, + AST; 
PY = 2.772 gry + (—0.022 + 0.024 i) ghg + (—0.025 + 0.021 i) gh, + AP’; 
ódon = 0.016, SN = 0.016; Shap =0028, du = —0.012, (154) 


where the QCD corrections are obtained by implementing the scaling factors Csfpf(Ttb) into the form 
TP we set 7. = 0. For the 


elw ? 


factors Fr p(T») and the electroweak correction 62° is from Fig. For 6 


elw 


numerical estimate of Neem: more precisely, we take the SM values of gyww = Gay Tax = 1. 


While, for the numerical estimate of Ton we assume a scenario in which the pseudoscalar form factor 
P7 is dominated by the top-quark contribution taking Iii = gi, = 0, i.e. neglecting the b- and 
T-loop contributions. In the SM, we have SY = —6.558 + 0.0477 and P” = 0. 


Finally, for the decay H — Zy, the scalar and pseudoscalar form factors are given by 
S47 = 12.372 gww + 0.689 g?r + (—0.019 + 0.011 i) 955, + (—0.0005 + 0.0002 7) gẹ, + AS”; 
P27 = 1.046 gig, + (—0.022 + 0.0114) gi, + (—0.0006 + 0.00027) gi, + AP, (155) 


in terms of the HWW and Higgs-fermion-fermion couplings. In the SM, we have $47 = —11.702 + 
0.0117 and P% = 0. The H —> Zy decay width is estimated at LO. 
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Table 6: Partial and total decay widths of the SM Higgs in MeV taking My = 125.5 GeV. We 
compare our numerical estimates with those presented in Ref. by introducing a quantity dp defined 
by dp = (T This Review — Ter (4) / Dee [B4- Also presented are theoretical uncertainties (THUs) of the 
partial and total decay widths from missing higher orders estimated around My = 125 GeV, see Tables 
178 and 182 in Ref. [84]. 


(fq > bb) I(H AWW) T(H- gg) T(H >rt) T(H > ce) 
This Review 2.367 9.185 x 107} 3.382 x 107} 2.572 x 107} 1.173 x 107! 
Ref. 2.387 9.222 x 107} 3.386 x 107} 2.573 x 107} 1.185 x 107! 
ôr [%] 0.8 -0.4 —0.1 =0.05 -I.0 
THU [%] 0.5 0.5 3.2 0.5 0.5 

A> ZA H>) H> Z) TH >p) Ti 
This Review | 1.139 x 107} 9.405 x 1073 6.531 x 1078 8.913 x 1074 4.129 
Ref. [84] 1.139 x 107} 9.438 x 1078 6.550 x 1078 8.927 x 1074 4.156 
ôr [%] —0.03 —0.4 —0.3 92 —0.65 
THU [%] £0.5 +1.0 +5.0 £0.5 £0.73 


Table 7: Branching ratios (BRs) of the SM Higgs taking My = 125.5 GeV. We compare our numerical 
estimates with those presented in Ref. [84]. Estimation of the total uncertainty THU+PU has been 
done by adding linearly the THU and the total parametric uncertainty (PU) where the latter is obtained 
by adding the individual PUs in quadrature, see Tables 174, 175, 176, 177, and 178 in Ref. [84]. 


BiH > bb) BIH WW) B(H- gg) B(H >rtr) BiH >c) 
This Review | 5.733 x 107} 2.225 x 107} 8.190 x 107? 6.230 x 107? 2.841 x 107? 
Ref. [84] 5.744 x 107}  2.219x 107} 8.147 x 107? 6.192 x 107? 2.852 x 107? 
THU+PU [%] 1.7 2.1 7 2.3 6.6 
B(H > ZZ) B(H >y) B(H > Zy) B(H > up)  TioølļMeV] 
This Review | 2.758 x 107? 2.278 x 1073 1.582 x 1073 2.159 x 1074 4.129 
Ref. 2.741 x 107? 2.271 x 1073 1.576 x 1073 2.148 x 1074 4.156 
THU+PU [%] 2.1 2.9 6.9 2.4 1.9 


In Table g we show the partial and total decay widths of the SM Higgs boson with My = 
125.5 GeV taking Ge = Gyyy = 1 and iF = AS97147 = API747 = 0. For a quantita- 
tive comparison with those presented in Ref. |84], we introduce ôr’s, which are defined by dp = 
(T This Review — TRer [B4)) /T per.[84) for each decay mode, for being contrasted with theoretical uncertain- 
ties (THUs) given in Ref. [84]. |°"| In terms of dp /THU, we find excellent agreement between our analysis 
and that in Ref. except for H — bb and H — cé for which we find 6p/|THU| ~ —2. The largest 
contribution to the discrepancy of Ty, comes from H — bb with the second (third) largest one from 
H — WW (cc). We note that our estimations of the decay widths into quarks are smaller than those 
in Ref. [84]. This might come from our incomplete and rough implementation of the ELW corrections 


of Eq. (102). 
In Table [7| we show the branching ratios and the total decay width of the SM Higgs boson taking 
My = 125.5 GeV. Again comparisons are made with those in Ref. together with the total uncer- 


27Note that we use the same values for all the input parameters as in Ref. [84]. 
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Figure 20: Higgs boson branching ratios and their uncertainties for the mass range around 125 GeV. 
The plot is taken from Ref. [84]. 


tainties. We pick up THUs and PUs from Tables 174, 175, 176, 177, and 178 in Ref. [84]. We note 
that the total uncertainty is about 2 - 3 % for H —> bb, WW, TT, ZZ, yy and uu. While it is about 7 
% for H —> gg,cc, and Zy. The total decay width is determined with about 2% error. In Fig. the 
branching ratios (BRs) are shown in the Higgs-boson mass range between 120 GeV and 130 GeV. For 
each BR line, the band width represents the corresponding total uncertainty. 


3.7.2 Decays of heavy Higgs bosons in CP-violating 2HDMs 


In this subsubsection, we study the decays of heavy neutral Higgs bosons appearing in BSM models. To 
be specific, we choose the type-I 2HDM identifying the lightest neutral Higgs boson as the SM-like 125.5 
GeV one discovered at the LHC. We assume that the lightest Higgs boson is purely CP even while the 
two heavier Higgs bosons exhibit nontrivial CP-violating (CPV) mixing in the presence of the complex 
quartic couplings of Às 6,7 in the Higgs potential. In this scenario, with no much need of decoupling the 
heavier Higgs bosons, all the branching ratios and the total decay width of the lightest Higgs boson 
remain consistent with those of the SM Higgs within the ranges allowed by the current LHC Higgs 
precision data [244]. The decay widths and branching ratios are calculated as summarized in Table 
but the ELW corrections are neglected for consistency. For the full model- and parameter-dependent 
ELW corrections in 2HDMs, see Appendix [F] 


To fix all the relevant couplings of three neutral Higgs bosons, one may start from the orthogonal 
3 x 3 matrix O describing the mixing among them. For CPV scenarios in 2HDMs, the three neutral 
Higgs bosons do not carry definite CP parities and they become mixtures of CP-even and CP-odd 


52 


states. In this case, without loss of generality, the mixing matrix can be parameterized as 


—Sa Ca O Cn 0 s} 1 0 0 
O = Ca Sa O 0 1 0 0 Gy Sy 
0 1 =S O Cn 0 -—s, Cy 
—Soly Calw + SaSnSw CaSw — SaSnCw 
= Can  Salw — CaSqSu SaSw + CaSqCw (156) 


— Sy = Bis CyCw 


introducing a CP-conserving (CPC) mixing angle a and two CPV angles w and 7. We recall that the 
mixing matrix O relates the electroweak eigenstates (%1 ,¢2,a) to the mass eigenstates (Hı , Ha , H3) 
via 

(Qı oy Q2 oy a)? = Oai( Hı 5 Hə 3 H3)? , 
with the ordering of My, < My, < My,. Assuming the lightest Higgs boson is purely CP even or 
taking s, = 0 and c, = 1, the mixing matrix takes the simpler form of 


— Sa CaCy CaSw 
Ca Sawu SasSw |- (157) 
0 -—-S, Gy 


O| 


Sņn=0,cņn=1 aa 


Note that, in the CP-conserving case, one of the heavy Higgs boson is purely CP odd and its coupling 
to a pair of massive gauges bosons is identically vanishing. We observe that H3 is purely CP odd when 
|cu| = 1 while Hz is CP odd when |s,,| = 1. Plugging the above expression of O into Eq. (45), the 
couplings of three neutral Higgs bosons to a pair of massive vector bosons are given by 


Juvy © p-a = V Le, Juvy T CB-alw = 02, Guzvv © CB-a8u = 03, (158) 


with 63 + 63 = e. We note that the two mixing angles are determined as follows 


69 oP) 
Sa — Leen Th a a a 
W W 


N 


Os Oe Og 02 
o= Se =) =e (159) 
ðs +03 € ðs +03 € 

in terms of the couplings 02 = gy, yvy and 03 = 9,,,yy together with tg. And then, the Yukawa couplings 
of the three neutral Higgs bosons are determined by 

ITiu = On. id — Dn = Ogzi/ S6 ; — Gita = Gira = inn = Oai/tg, (160) 
where u and d stand for the up- and down-type quarks, respectively, and £ for three charged leptons. 
To summarize, in the scenario under consideration, all the Yukawa couplings of the two heavy Higgs 
bosons could be fixed by giving their couplings to the massive vector bosons. On the other hand, 
depending on sign[d2/c,,], all the Yukawa couplings of the lightest Higgs boson are determined by 
0401/88 = Co/S3 = V1 — €+ Ve/tg which, especially for large tg, approaches the SM value of 1 as 
quickly as the ga yy = V1—e coupling when e goes to zero. This is the very reason we choose the 
type-I 2HDM for our numerical study avoiding conflicts with the current LHC Higgs precision data [244]. 


?8Here we take the abbreviations such as cosa = Ca, Sina = Sq, etc. 
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Figure 21: The total decay widths and branching ratios of the two heavy Higgs bosons in the type-I 
2HDM taking tg = 5 and |In mm | = |9u2,2,| = 0-1. In the left CPC panels, we take 63 = gp yy = 
2eo (Mh, /Mh,) = € and 63 = Ue = 0. While, in the right CPV panels, we take 63 = ee = 


€o (Mh, /Mh,) and 63 = Ce, = eo (Mh, /Mh,). In the both CPC and CPV cases, we take eo = 0.05. 


For our numerical study, we vary tg but, for d23, we are taking 


Mn, ° Mm i 
Ju vv = do = vy €o (57 ) ’ Juvy = 63 = Vy €o (i ) ; (161) 
3 


Hə 


reflecting the behavior of €e which is suppressed by the quartic powers of the heavy Higgs-boson masses 
at leading order [244]. With the above parameterizations of 623 and taking c, > 0 we have d2/c, = 
Cp—-a = Ve > 0 and s2 = My, /(My, + Mi) leading to a maximal CPV mixing between the two heavy 
Higgs bosons when they are degenerate. For €9 or the largest possible value of 65 35 we choose a value 
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Figure 22: The same as in Fig. [21] but for tg = 30. 


which is a little bit larger than the lower lo error of C, in the CPC4 fit: 
€o = 0.05, (162) 


having in mind the relation gp yy = y1 — e > V1 — 2e ~ 1 —e. For the masses of Higgs bosons, we 
take 
My, = 125.5 GeV, My, = Mp, — 50 GeV, My, < Mps ~ 1 TeV, (163) 


with My, varied. This choice may result in the simpler decay pattern by forbidding or suppressing the 
decay channels of Hy; > H*W**, Ho3 > H*H*, H; + H, Ho, H3 + HoH, etc. By assuming very 
heavy charged Higgs boson, also neglected are the contributions from the charged-Higgs-boson loops to 
the decay processes of the heavy neutral Higgs bosons into yy and Zy. Finally, for H23 > Hy A, 


29 See E 2 ie a eee Note Cy = Cw = Cz, see Eq. (175). 


3°For the details of the contributions from the charged-Higgs-boson loops to the neutral Higgs boson decays into two 
photons in the 2HDM and the MSSM, see Appendix [E] 
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we take |Jan | = |\9u,2,1,| = 0-1. For the rigorous treatment of the cubic HH, Hı and H3H, A; 
self-couplings expressed in terms of the masses of charged and neutral Higgs bosons and the elements 
of the mixing matrix O, see Appendix [E] 


In Fig. 21] we show the decay widths and branching ratios of the le heavy Higgs bosons in the 
type-I 2HDM taking tg = 5. For CPC, H; is taken to be CP odd with a. a = 0 and, accordingly, the 
decays of Hz into WW, ZZ, and HH, are forbidden. Incidentally, we note that H> — HZ decay is 
also forbidden since a m g, yy: On the other hand, for CPV, there are no forbidden decay modes 
as long as they are kinematically allowed. In the CPC case, the total decay widths of CP-even H and 
CP-odd H; are largely enhanced at the Hı Hı and tt thresholds, respectively. While, in the CPV case, 
both of the thresholds contribute to the total decay widths as shown in the upper panels of Fig. 


For CPC, we further observe that the couplings of the CP-even state of H3 t oe are identically 
vanishing when vee = 0. It does happen at sa = 0 or e = 1/(1+1%), see Eq. . This explains why 
there are ask for the fermionic decay modes of the CP-even Hə state at My, = a + t3))'/4Mn, ~ 
70 \/ts GeV }*'|as found in the middle-left panels of Fig. |21| and Fig. [23] 22| around My, = 160 GeV and 
385 GeV, respectively. We note that the branching ratios of fermionic decay modes H33 — tt, bb, TT, pp 
are smaller for the larger value of tg since the corresponding decay widths are suppressed by the factor 
of ~ 1/t3. For large values of My,,, the numerical results are consistent with the observation that 
the decay width of fermionic decay modes is proportional to My, while that of bosonic decay ones is 
inversely proportional to My, , especially with the parameterization of Eq. for Hy3 > VV , AZ. 
On the other hand, the decay width [(H3 + HZ) is not suppressed by the heavy mass My, because 
of the relation g, “a es yy: im our numerical study, it is suppressed since we have taken the small 
mass difference between Hs and H of My, — My, = 50GeV < Mz. Otherwise, it may increase in 
proportion to Mẹ,- 


4 Decays of a Charged Higgs Boson 


The effective couplings of the charged Higgs boson H* to quarks and leptons are described by the 
interaction Lagrangian: 


23 m m 
Lutfrfy = V2 Ht, (an, Pu +E gp Pr) fi + h.c., (164) 


where Prr = (1 F ys)/2 and (fy, fi) = (t, b), (c, 8), (r, T), (pu; H), etc. The masses of the up- and 
down-type fermions are denoted by my, and m,,, respectively. On the other hand, the interaction of 
the charged Higgs boson with a massive gauge boson W and a neutral Higgs boson H is given by 

g 


oO 
Lunswt = —59uqtw-W "(Hi dy H*) + he, (165) 


with the convention X Bu Y = X(0,Y)—(0,X)Y. With these effective interactions, in this section, we 
study the charged Higgs decays at LO except for the QCD corrections considered in the decay modes 
into quarks, concentrating on mainstream instead of being complete. For comprehensive studies of 
decays of a charged Higgs boson in BSM models, we refer to, for example, Refs. B51]. 


When a charged Higgs boson decays into quarks, the decay width is given by 


ff, Mae 
ria > ffi) = Nọ t Srv? A (1, K fyo Kf) l (1 TRAT Kf) (mo, + m39) 


31In CPC, note that we take 3 = 2eo (Mh, /Mh,) = € with 63 = 0. 
32For H+ — tb, note that the QCD corrections are not valid in the threshold region due to the top-quark mass 
effects. For them, we refer to [87] and references there in. 
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—4 EAR fy mms snan } (1 + dacp) ; (166) 


including the QCD correction factor dacp with Ky, 7, = Mj}. p,/Mj+ and nit * = 3. Note that the MS 
quark masses such as mM c(My+) and ™},5(My+) are taken for ms, and mp, respectively. While, for the 


decays into leptons, the charged lepton pole masses are used for my, together with ms, = 0, Ne hh =, 
and dacp = 0. A charged Higgs boson may decay into a lighter neutral Higgs boson H snd a massive 
gauge boson W and the LO decay width is given by [173} [166] 174] 


Porat > HW**) = 


WE—-/w)? 
GrMF2|Ganey—l vw) ia ewà? (ws, w, £) (167) 


827 0 wm[(a — 1)? + ey] ’ 


with w = M / MẸ, and w+ = M? /MẸĝ,. When Myz is larger than the sum My + Mw, it reduces to 


Gr M? 
p/n autw- 


where ky = M3 [M= and kw = Mi, /M Wes The electroweak corrections within the 2HDM framework 
have been calculated, specifically for the process of H+ > Wth,W*+H,W=+A 254|. They are 
of moderate size and numerically stable if a process- and gauge-independent renormalization scheme is 
chosen [254]. Note that the package 2HDECAY [255] is available for the calculation of full ELW one-loop 
corrections to both neutral and charged Higgs decays in the CP-conserving 2HDM. 


Ton? => HW?) = ?\?/2(1, KH, KW), (168) 


Finally, at LO, the decay widths of a charged Higgs boson H* into a chargino ra and a neutralino 
X? are given by 


2 Da gee aes 
Lorra ~o — J Matà" (1, Ki, Kj) 
D(AT >X Xi) = T (169) 


x [a- Ki— Kj) (lI zog? + Izol? ) — 2y Fiki (lIr zoz [? — [Izox BE 


with ki = Meo /M74 and kj = me /M?,. and the LO decay widths into a pair of sfermions by 
t j 


~ak „v? zl 
T°(H* > fifi) = NE — N, ki kj), (170) 


where kK; = M /M? +, kj = M P, /M2.., and NÉ ‘= 3 and 1 for squarks and sleptons, respectively. 


For a numerical example, we take the 2HDMs in which the relevant couplings are given by P| 


7 1 1 
Ht tb, : g,=—, gp =- (L,I), tg (II,IV); 
tg tg 
i i 1 
HY —> vr” ,vu™ : gp, =0, gn = -7 (LIV), te (L, ID; 
tg 
H++HWt : |g, Ê =1- g9}, (L, 11,11, IV), (171) 


depending on the 2HDM type as denoted by I, H, HI or IV. For the decay H* + HW+, we set 
My = 125.5 GeV and |g py- = 0.1. For the MS mass of the strange quark, we take the 


33 See Eq. and Table f2} 


34For tha we cand Invy œ 0.95 adopting a little bit larger value than the lower lo error of C, in the CPC4 fit, see 
Section |5.2)and Table|11}there in. 
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approximation M.(u) = Me(u)/11.72 [6]. On the other hand, we neglect the effects of its pole mass 
which is relevant to the H+ — cs decay only through the kinematical factor xs, see Eq. (166), and too 
tiny to influence the numerical results. 


For the charged Higgs-boson decay H+ — tb, we take the contribution of the off-shell top quark 
into account [166} [74]: 


7 7 1 Mp (Mpt -M)? M2 2 2 
aa a E E e Ht f Ae pe (1- ) (2+ r) 


973 
U 2°71 M2, 


(p? — MP)? + MPT? a 


where the kinematical b-quark mass is neglected in the t — bW* decay process and the triangle function 
A+ is given by 


Nyt = 14+ 07 + 0% — 2a; — 2ag — 20:03, (173) 
with a; = p?/Mj,+ and aj = M?/M7,2. When My > M; + Mp, using Eqs. (110) and (112), we have 


LO 
roc > tb = bw) = E > 15) EO) 
t 


(174) 


in a factorized form. 


In the upper panels of Fig. we show the decay widths of a charged Higgs boson into two quarks, 
tb (left) 5| and c3 (right) in the type-I and type-II 2HDMs taking account of the QCD corrections. 
In these models, g1,¢ = —9s,s = 1/tg and the top-quark and charm-quark contributions dominate and 
the decay widths scale as 1/t3. In the lower panels of Fig. we consider the type-II and type- 
IV 2HDMs in which gee = 1/tg and g,, = tg. For low values of tg, the top-quark and charm-quark 
contributions dominate and the decay widths scale as 1/ t. On the other hand, for high values of tg, the 
bottom-quark and strange-quark contributions are enhanced by the factor of t while the top-quark and 
charm-quark contributions are suppressed by the factor of 1/ t3. The bottom-quark and strange-quark 
contributions start to dominate when t, is larger than \/7™,/™», ~ 7.7 and 4/Me/Ms ~ 3.4, respectively. 
When the bottom-quark and strange-quark contributions dominate, the decay widths scale as (t3/60)? 
and (tg/12)?, respectively, compared to the case with tz = 1 and these factors are responsible for the 
significant change of the decay widths for high tg values, as can be checked by comparing the lines with 
tg = 1 and those with tg = 10 and 30. 


In Fig. we show the LO decay widths of a charged Higgs boson with mass Mj+ into rv in 
the type-I/IV 2HDMs (left) and in the type-II/III ones (right) taking tg = 0.5 (blue dashed), 1 (black 
solid), 10 (red dotted), and 30 (magenta dash-dotted). In the type-I/IV and type-II/III 2HDMs, the 
decay widths scale as g? = 1 /t3 and t2, respectively. Compared to the case of H* — rtv, the decay 
widths ['(H*+ — u*v) are simply suppressed by the large factor of m?/ m? ~ 300. 


In Fig. we show the LO decay widths of a charged Higgs boson with a mass My=+ into HW* taking 
My = 125.5 GeV and lI antw- |? = 0.1. In the right panel, we magnify the low My region covering the 
case with a virtual W** and we compare the decay width with the prediction of the two-body decay 
width (dash-dotted magenta line). We note that they are nearly identical for My+ > My + My, as 
expected. 


35At My+ = 500 GeV, we switch from the three-body decay width T(H+ — bbW*) to the two-body decay width 
I'(H* — tb) because of the same reasons as in H — t*t* described previously. 
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Figure 23: (Upper) Decay widths of a charged Higgs boson with a mass My=+ into tb (left) and cs 


(right) in the type-I and type-II] 2HDMs. QCD corrections are taken into account. For tan 3, we take 
four values of 0.5 (blue dashed), 1 (black solid), 10 (red dotted), and 30 (magenta dash-dotted). (Lower) 
The same as in the upper panels but in the type-II and type-IV 2HDMs. 
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Figure 24: (Upper) Decay widths of a charged Higgs boson with a mass My=+ into Tv in the type-I/IV 
2HDMs (left) and in the type-II/III ones (right). For tan 6, we take four values of 0.5 (blue dashed), 1 
(black solid), 10 (red dotted), and 30 (magenta dash-dotted). 


In Fig. we show the total decay widths of a charged Higgs boson for four values of tg: tan 8 = 0.5 
(blue dashed), 1 (black solid), 10 (red dotted), and 30 (magenta dash-dotted). We sum over the 5 decay 
modes of Ht —+ tb, Ht > cs, Ht —> tv, Ht —> pv, and Ht —> HW” with My = 125.5 GeV and 
= 0.1. 


eee |? 
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Figure 25: Decay widths of a charged Higgs boson with a mass My+ into HW* taking My = 125.5 GeV 
and |g paw- |? = 0.1. In the right panel, we magnify the low My region and compare with the two-body 
decay width (dash-dotted magenta line). The vertical line locates the position My+ = My + My. 


Firstly, we consider the total charged-Higgs decay width in the case of a heavy charged Higgs boson, 
see the left panels of Fig. We note that the top-quark contributions to [(H+ — tb) scale as 1/ th, 
independently of the 2HDM types. As the value of tg grows, the total decay width monotonically 
decreases in the type-I and type-II] 2HDMs. In contrast, in the type-I] and type-IV 2HDMs, the 
decay width reaches the minimum around tg = 8, where the top-quark and bottom-quark contributions 
become comparable, and then it starts to increases as tg grows further. Based on these observations, we 
can work out the behavior of the total decay width of a heavy charged Higgs boson for large values of tg: 
the dominant contribution comes from the Ht — HW* decay mode with the subleading/competing 
contributions from H+ — tb. We identify that there exists an additional subleading contribution from 
H* —> rv for tg = 10 and 30 in the type-II and type-II 2HDMs. On the other hand, for low values 
of tg, the H+ — tb decay mode dominates the total decay width with the subleading contribution 
from Ht — HW* which amounts to about 30 GeV at My+ = 1 TeV taking My = 125.5 GeV and 
= 0.1, see Fig. 


CET |? 


Secondly, we consider the case of a charged Higgs boson lighter than t quark, see the right panels of 
Fig. For large values of tg, the Ht — rv decay mode dominates in the type-II and type-II 2HDMs. 
In the type-I 2HDM, the H* — rv decay mode still dominates with the subleading contributions from 
HY — cs though both of them are suppressed by 1/t3. In the type-IV 2HDM, the H* — cs decay 
mode dominates where the strange-quark contribution is enhanced by t. For small values of tg, the 
H* — rv decay mode mostly dominates, leading to the larger decay widths for tg = 0.5 in the type-I 
and type-IV 2HDMs in which |g,| = 1/tz. 


In Fig. we show the branching ratios of a charged Higgs boson varying its mass between 100 
GeV and 1 TeV in the type-I and type-II] 2HDMs taking tan = 0.5, 1, 10, and 30 from top to 
bottom. We consider 5 decay modes of H+ — tb (red solid), H+ — cS (magenta solid), H+ > rv 
(blue solid), H+ — uv (green solid), and Ht — HW? (black dashed) with My = 125.5 GeV. In the 
type-I and type-III models, the normalized charged Higgs couplings to quarks are the same because 
Jte = |9v,s| = 1/tg while those to leptons are given by ga, = —1/tg (type I) and g,- = tg (type III). 
We observe that, especially for large values of tg, the heavy charged Higgs boson dominantly decays 
into HW*+ (black dashed) since ['(H* — tb) is suppressed as tg grows. And, when the H+ > HW+ 
decay mode dominates, the leptonic branching ratios are enhanced by t in the type-II 2HDM. In the 
type-II 2HDM, we note that B(H*+ — rv) becomes comparable to B(H+ — tb) around tg ~ 10 and it 
becomes larger as tg grows. 
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Figure 26: Total decay widths of a charged Higgs boson with a mass My+ in the type-I, type-II, 
type-II, and type-IV 2HDMs from top to bottom. We are taking tan 8 = 0.5 (blue dashed), 1 (black 
solid), 10 (red dotted), and 30 (magenta dash-dotted). In the right panels, we magnify the low My 
region. 


In Fig. we show the branching ratios of a charged Higgs boson for its mass between 100 GeV 
and 1 TeV in the type-II and type-IV 2HDMs taking tan 6 = 0.5, 1, 10, and 30 from top to bottom. 
We consider the same 5 decay modes as in Fig. In the type-II and type-IV models, the charged 
Higgs couplings to the (right-handed) up-type quarks are again given by gie = 1/tg but those to the 
down-type fermions are by gb, = 9,7 = tg (type II) and gb s = 1/|g,,,7| = tg (type IV). We observe that, 
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Figure 27: Branching ratios of a charged Higgs boson with a mass My+ in the type-I (left) and type-I 
(right) 2HDMs taking tan 6 = 0.5, 1, 10, and 30 from top to bottom. We consider 5 decay modes of 
H* — tb (red solid), H* — cs (magenta solid), H* — rv (blue solid), H+ — uv (green solid), and 


H+ — HW? (black dashed) with My = 125.5 GeV and |g 


=0.1. 


|2 
HH+w- 


especially for large values of tg, the heavy charged Higgs boson dominantly decays into HW* (black 
dashed) eventually. But, compared to the type-I and type-II 2HDMs where the top-quark contributions 
to the H+ — tb decay always dominate and the decay width decreases as tg grows, the dominance of 
the H* — HW? decay mode develops rather slowly because the bottom-quark contributions take over 
the dominance around tg = 8 and the partial width '(H+ — tb) increases as tg grows in the type-II 
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Figure 28: The same as in Fig. [27] but for the type-II (left) and type-IV (right) 2HDMs. 


and type-IV 2HDMs. The leptonic decay widths are also enhanced by t in the type-II 2HDM and we 
find that, specifically for tg = 30, the three decay modes of H+ > tb, Ht > HW*, and H+ — rv are 
competing in the region of 350 £ My+/GeV XS 600. 


For a light charged Higgs boson, it mostly decays into rv and/or c5 before the H+ — t*b decay 
channel opens and starts to dominate. Exceptions occur in the type-I and type-II 2HDMs when tg is 
large, see the lowest panels of Fig. Specifically, for tg = 30 in the type-I model, the charged Higgs 
boson dominantly decays into HW in the narrow region of My=+ between 150 GeV and 180 GeV. 
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Figure 29: Decay widths of a charged Higgs boson with a mass My+ into H'W* for three values of 
My = 300 GeV (blue dashed), 500 GeV (red dotted), and 700 GeV (magenta dash-dotted) taking 
lJ 1w- |° = 1. For comparisons, T1O(H+ > HW*) with My = 125.5 GeV and San = 0.1 is 
also shown in black solid line, see Fig. 


|? 
Htw- 


We confirm that the behavior of the branching ratio of each decay mode versus My+ does not 
depend on the 2HDM type when tg = 1 as it should be: see the tan 8 = 1 panels in Figs. [27] and [28] 


Before closing this section, we address the case in which there exists another neutral Higgs boson H’ 
with its mass smaller than Mj+ and the decay Ht + H'W* is kinematically allowed. Here H’ stands 
for a mixture of CP-even and CP-odd states in general which, in the CP-conserving case, could be either 
of them. In the 2HDM framework, identifying the lightest neutral Higgs boson as the SM-like H with 
My = 125.5 GeV and lI antw- |? =1- Ivy = 0.1 as already considered, the charged Higgs coupling 
to the heavier neutral Higgs boson H’ and a charged vector boson W takes almost the maximum value 
of 1 due to the sum rules given in Eq. (46). Furthermore, the decay width grows by the cubic powers of 
the charged Higgs-boson mass. These combined properties result in a large width for the decay mode 
H= — H'W*= when H= is heavy enough. When My+ = 950 GeV, we observe that the decay width 
[4°(H* — H'W*) is larger than 200 GeV for a H’ boson with My = 300 GeV, see the blue dashed 
line in Fig. Even for H’ as heavy as 700 GeV, we note that the decay width could be comparable to 
[4°(H+ + HW*) when My+ ~ 1 TeV as denoted by the magenta dash-dotted line in the same figure. 


5 Higgcision of the Higgs Boson Discovered at the LHC 


f le u 1 2018, based on the works | 257, [258]. For several TEET 
epee = approaches pith at a own a on aes he global ft fits oft a a couplings, see, for example, 
Refs. [84 88, [271]. For some other early works 
based on REE ER PENTE S S OE we r E m S T 274, 275, 276, 277, 278, 279, 280, 
[281] (2821 283, [284] [285] 286) [287] 283]. 


Assuming generation independence for the normalized Yukawa couplings of Tarp we use the following 
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C notations for the couplings in the global fits: 


S S S S S S 
Ca = GHiu> Cà = Hda: Ci = Guin Cu = Giwwa C= Juggs 
ie P P P P P 
Ca = GHiiu> C] = J Hda > C = DH: (175) 


We further keep the custodial symmetry between the W and Z bosons which leads to the relation 
Cy = Cw = C. We emphasize that one can carry out the precision study of the LHC Higgs data 
using the model-independent approach taken in Section [3] for decays of a generic neutral Higgs boson 
by matching the relevant couplings and fitting parameters as in Eq. (175). 


5.1 Higgs signal strengths 


Each theoretical signal strength can be written in a product form as 
AP, D) ~ uP) AD), (176) 


where P = ggF, VBF, VH, ttH denote the production mechanisms and D = yy, ZZ“), WW“), bb, rtr- 
the decay channels, which are experimentally clean and/or dominant for My ~ 125 GeV. The factoriza- 
tion assumption is valid only when the production and decay processes are well separated like as in the 
resonant s-channel Higgs production in the NWA. By factorizing them, non-resonant and interference 
effects are inevitably neglected. More explicitly, at LO, the production signal strengths are given in 
terms of the relevant form factors and couplings by 


|S°(Mn)| + |P9(Ma)I’ 


a(ggk) = 
|SS Ma) I? 
L(VBF) E u(VH) = ee ’ 
A 2 2 
AttH) = (gre) + (Gnu) > (177) 
and the decay signal strengths by 
= B(H > D) 
D) = ————— 178 
with the branching fraction of each decay mode defined by 
Tr(H > D 
n 5) E (179) 


7 Tiol H) T AT tot 


Note that an arbitrary non-SM contribution AT,,; to the total decay width is introduced. We observe 
Ttot(H) becomes the SM total decay width when g7, z =l, inf = 0, Juwwnzz = l, and Ao = 
AP1547 = 0. Note that the LO relations in Eq. (177) are most reliable when higher order corrections 
to a BSM production cross section and those to the corresponding SM one are the same and so they 
are canceled out in the BSM-to-SM ratios. Otherwise they are valid at LO strictly. 


On the experimental side, we use the direct Higgs signal strength data collected at the Tevatron and 
the LHC. Specifically, we use 3 signal strengths measured at the Tevatron, see Table |8| At the LHC 
with the center-of-mass energies of 7 and 8 (7@8) TeV, the signal strengths obtained from a combined 
ATLAS and CMS analysis [291] are used, see Table|9} On the other hand, the 13 TeV data are still given 
separately by ATLAS and CMS and in different production and decay channels. |?°| Under this situation, 
to derive the combined signal strengths of various channels, we use a simple y? method assuming that 
each distribution is Gaussian. The results are shown in Table 


36For the details of the 13 TeV data sets used, see Appendix B of Ref. [258] and references therein. 
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Table 8: (Tevatron: 1.96 TeV) The signal strengths data from Tevatron (10.0 fb~! at 1.96 TeV). 


Channel Signal strength » My(GeV) Production mode Xy4 (each) 
c.v + error ggF VBF VH ttH 
Tevatron (Nov. 2012) 
Combined H + yy[289] 6.141379 125 78% 5% 17% - 2.60 
Combined H + WW [289] 0.854983 125 78% 5% 17% - 0.03 
VH tag H — bb[290] Teor 125 - - 100% - 0.67 


X24 (subtot): 3.30 


Table 9: (LHC: 768 TeV) Combined ATLAS and CMS data on signal strengths from Table 8 of 
Ref. [291]. 


Decay mode 
Production mode | H > yy H => ZZ® HoWW” H- bb H > rtr 
ggF LHe 113ta 0.84 TF - 1.075% 
VBF 1.3705 0.170% 1.276% - 1.3794 
WH 0.5575 - 1.6719 1.070% —1.4t74 
ZH 0.5430 - 5.9433 0.4404 2.2773 
ttH 9259 - 5.0418 hie -1.9137 
Xéļm(subtot): 19.93 


Table 10: (LHC: 13 TeV) Combined ATLAS and CMS (13 TeV) data on signal strengths. The 
UIE pinea (UPES inea) represents the combined signal strength for a specific decay (prorietinn) channel 
by summing all the production (decay) modes, and y2,,,, are the ee minimal y? values. In a 
VH/WH row, the production mode for H > yy and H + ZZ“) is VH while it is WH for H > WW“ 
and H > TtT: for the remaining decay mode H — bb, we combine the two signal strengths from WH 
and VH, see Table XII in Ref. [258]. 


Decay mode 

Production mode | H>yy H—>ZZ® H—=>WW® H=>bb Harte | Eome Xun) 
ggF 1.024917 1.09491 Loo 2.514263 1.064939 | 1.114807  5.42(3.15) 

VBF 123535, TSI 0.5453: g 11503, | 1:027} T5375) 
VH/WH 14a, OTI R E 1.074923 3.3918% | 1.154670  7.05(6.44) 

ZH - - 1:001 1.2049:33 1.231162 | 1.19+932  0.45(0.02) 

ttH 1.364637  0.00+6:80 - 0.911015 - 0.931954  5.96(5.86) 

ttH (excl.) 1.397545 1.594244 0.774936 Der | 1.16+922  4.17(3.62) 

USS brea Lig. L TARE T 1.207015 150 Tag | 110790 

Ku a) 6.83(5.72)  9.13(8. 88) 9.48(7. 32) 1.56(1.51)  3.58(3. 20) 30.58(27.56) 


5.2 Constraints on the couplings of Higgs boson weighing 125.5 GeV 


In this subsection, we present a few representative results obtained from performing LO analysis of 
the direct Higgs data collected at the Tevatron and the LHC by considering CP-conserving (CPC) 
scenarios only. Note that, in the most general CPC scenario, one may vary all the 7 parameters of C?, 
CS, CP,C,, ASI, AS7, and AT, while taking vanishing pseudoscalar couplings and form factors as 
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Table 11: (CPC) The best-fitted values in various CP conserving fits and the corresponding x? per 
degree of freedom (dof) and goodness of fit. The p-value for each fit hypothesis against the SM null 
hypothesis is also shown. For the SM, we obtain x? = 53.81, x?/dof = 53.81/64, and so the goodness 
of fit = 0.814. From Ref. [258]. 


Cases CPC1 CPC2 CPC4 CPCN4 
Varying AT tot ASY c3, 03, O30, 
Parameters ASI C? C, AS7, ASI 
cS 1 1 1.001+9:956 | 1.042+9:977  1,042+9:078 Lae  —1,042+9-081 
CF 1 1 0.96273 101 1 1 1 1 
CF 1 1 1024 Foes 1 1 1 1 
Co 1 1 OIG ne | 1.02710034 Lee  1.028t0:035  7.028+9:034 
As? 0 —0.2261 0-33 0 0.129703 0.120103  3.524+0-41 3:523 83 
ASI 0 0.016+9:025 0 —0.02119:057 1.341906  0.095+9:955  1.414+9:985 
AT tot (MeV) | —0.28510'78 0 0 0 0 0 0 
x?/dof 51.44/63 51.87/62 50.79/60 50.96/60 
goodness of fit 0.851 0.817 0.796 0.791 
p-value 0.124 0.379 0.554 0.583 
4 75 0 05 1 
AT ot [MeV] 
Figure 30: CPC1: Ay? from the minimum versus AT;,, with only AT;,¢ varying in the fit. From 


Ref. [258]. 


CP =CP =CP =AP1=AP9 =0. 


Our goal is to provide constraints on the couplings of the neutral Higgs boson, which was discovered 
at the LHC, without much loss of generality when it is interpreted in various frameworks beyond the 
SM. Accordingly, we consider the four CPC fits listed in Table [M] in which the second row explicitly 
shows the varying parameters of each fit. 


Referring to Ref. for more detailed explanations, we offer highlights of each fit as follows: 


e CPC1: The best-fit value for the residual total decay width AT,,¢ is Ato, = —0.285 ae MeV 
which is 1.60 below zero. At 95% confidence level (CL), on the other hand, 


AT tot = —0.28513:33 MeV, (180) 
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Figure 31: CPC4: (Upper) The confidence-level (CL) regions of the fit by varying C,, C$, Cf, and 
C$. The contour regions shown are for Ay? < 2.3 (red), 5.99 (red+green), and 11.83 (red+green+blue) 
above the minimum, which correspond to confidence levels of 68.3%, 95%, and 99.7%, respectively. The 


best-fit points are denoted by triangles. (Lower) Ax? from the minimum versus Yukawa couplings. 
From Ref. [258]. 


as shown in Fig. [80] Using the upper error as the upper limit, we obtain the constraint AT, < 0.38 
MeV simply taking the central value equal to zero which can be translated to the following 
95% CL constraint on the branching ratio of the Higgs boson decays into non-SM particles: 
B(H — nonstandard) < 8.4%. 


e CPC2: The best-fit point (AS7, AS) = (—0.226, 0.016) of the form factors indicates an increase 
of 3.4% and 2.4% in the absolute values of the scalar Hyy and Hgg form factors of |S7| and |S], 
respectively. We note that the error of AS? is +0.025, which is numerically smaller than the SM 
bottom-quark contribution of —0.043 to the real part of S9, see Eq. (153), alerting that we have 
already reached the sensitivity to probe the sign of the bottom-quark Yukawa coupling in gluon 
fusion. 


e CPCA4: We observe that the possibility for the top-quark Yukawa coupling to be negative has 
been entirely ruled out as shown clearly in the left upper and lower panels of Fig. And, as 
already anticipated in the CPC2 fit, the bottom-quark Yukawa coupling C? prefers the positive 
sign to the negative one, see the middle panels of Fig. It is more clear from the middle lower 
panel that the point C? = —1 has Ay? > 2 above the minimum at C? = +1. The current data 
precision is yet insufficient for showing any preference for the sign of tau- Yukawa coupling, as 
shown in the right panels of Fig. 


e CPCN4: In this fit, there are 4 degenerate minima with AS! ~ 0,+1.4 for CS ~ +1 which 
could be understood from the relation |S9| ~ |0.7C% + AS9| ~ 0.7, see Table Note that AS? 
can compensate the sign change in C$ allowing it to be about —1 with AS? ~ +3.5. This could 
be understood by noting the relation |S7| ~ |—8.3C,+1.8 C3 + AS"| ~ 6.5. We further observe 
that the negative top-quark Yukawa coupling is allowed only when there exist additional particles 


37For this reason, we have considered the minimum around Cr = +1 only in Table[11] 
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Figure 32: (Left) The allowed parameter space on the (My, |guyy|) plane from AT, < 0.38 MeV at 
95% CL. Also shown are future prospects assuming two and four times stronger constraints on AT tot. 
(Right) The same as in the left panel but for (My, |guyz|). 


running in the H-y-y loop with the size of contributions equal to two times the SM top-quark 
contribution within about 10 %. This tuning on the couplings and form factors will become more 
and more severe as more data are accumulated at the LHC. 


5.3 Implications beyond the Standard Model 


In the cxSM where the SM is extended by adding a complex SU(2), singlet, there could be a light Higgs 
boson y mainly from the singlet sector and the 125.5 GeV Higgs boson H may couple to a pair of them 
through the singlet-doublet mixing term in the potential. When kinematically allowed or 2M, < 125.5 
GeV, H decays into a pair of light scalars and, in this case, AT',o, < 0.38 MeV may provide constraints 
on the mass M, of the scalar particle and the absolute value of the coupling gy,, at 95% CL. We find 
lJHep| S 0.005 for M, < 40 GeV, see the left panel of Fig. In 2HDMs, assuming that the scalar 
y is CP odd and the Higgs boson H with My = 125.5 GeV is CP even, we consider another decay 
H — yZ* to constrain the coupling guez depending on Mọ. We find |gyyz| S 0.5 for Mọ S 50 GeV, 
see the right panel of Fig. 


In the MSSM, there still exists a room for the lightest neutralino Y? to be light with its mass two 
times smaller than My = 125.5 GeV [292]. In this case, again, AT to. < 0.38 MeV may provide 95% CL 
constraints on the mass mzo and the absolute value of the coupling Pai assuming H is purely CP 


X1X1 
even. We find Sines < 0.01 for mgo £ 45 GeV, see Fig. 


From Table we find |AS7| < 0.4 and |AS9| < 0.06 by taking the errors in the CPCN4 fit. 
Exploiting these constraints, we derive constraints on the H couplings to a pair of charged Higgs 
bosons, lighter charginos, and the lightest stops/sbottoms. For the SUSY contributions to AS? and 
ASI, we refer to Appendix [B] 


In the left panel of Fig. we show the contour lines for |AS7| = 0.4,0.2,0.1 assuming it is 
dominated by the contributions from the lighter-chargino loops and considering only them. We find 
that the mass of the lighter chargino is constrained to be mz+ Z 390 GeV by the current LHC Higgs 
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Figure 33: The allowed parameter space on the (mzo, | Iroz) plane from AT, < 0.38 MeV at 95% 
1 
CL. Also shown are future prospects assuming two and four times stronger constraints on AT ot. 


Figure 34: 
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(Left) Contour lines for |AS7| = 0.4,0.2,0.1 from left to right on the (Mgt, Duxtel) 


plane assuming the lighter-chargino-loop contributions dominate AS7. (Right) The same as in the left 
panel but on the (My+, |guH+H-|) plane now assuming AS? is dominated by the contributions from 
charged-Higgs loops. For the reference value of |AS7| = 0.4, we are taking the lo error of the CPCN4 


fit. 


data when the relevant coupling is assumed to be 1. The lower-bound constraint on m + linearly 
increases as the bound on |AS7| becomes stronger. In the right panel of Fig. we show the same 
contour lines now assuming that AS7 is dominated by the charged-Higgs loops. In the 2HDM and 
the MSSM, keeping the most significant three contributions when My, , Z v, we find that the a 
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Figure 35: Contour lines for [AS] = 0.06 , 0.03 ,0.015 from left to right on the (mz, |9;77,|) plane 
assuming AS9 is dominated by the contributions from the lighter-stop loops. For the reference value of 
|AS9| = 0.06, we are taking the 10 error of the CPCN4 fit. 


coupling is given by 


30 ( “tH 45 o (ta) _ (400GeV)* (ts —1) N| (1+ t5 \ | Re(mie’*) 
Junta- 1 TeV 30 Mog z 30 F (100 GeV)? 
400 GeV\? (1+#3\ /t\? 
apo ae i£ 181 
ae) ( 7 ie Eau Ei 


under the assumption that Hı is purely CP even state. Note that the second and the third terms 
contain contributions enhanced by the factors of tg and t and, for the tZ-enhanced contributions, 
we are taking My, ~ My, = Mz and cg. = Ve = ye +03 = V2€9 Mj, /M33 with e = 0.05 
as in the subsubsection |3.7.2| see Eq. (161). It is worthwhile to note that AS](H*) — 1/3 when 
Mis Min, > MZ, v’? , Re(mi,e*), see Eq. and the two equations following it. The constraint 
on My+ becomes two times stronger when the bound on |AS7| becomes four times stronger. 


In Fig. we show the contour lines for |AS9| = 0.06 ,0.03,0.015 assuming that it is dominated 
by the contributions from the lighter-stop loops. The current LHC Higgs data constrain the mass of 
the lighter stop as mz, Z 1 TeV only when the relevant coupling is as large as 10 and, otherwise, it is 
weaker or much weaker than several direct bounds on the lighter-stop mass [6]. 


6 Conclusions 


38See the first relation in Eq. (E.8) in Appendix 
39See also https://pdglive.lbl.gov/DataBlock.action?node=S046STP. 
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as possible which induce CP-violating mixing among the neutral Higgs bosons. We have included the 
relevant higher order corrections, which are dominated by QCD corrections, as much as possible. 


Firstly, the comprehensive analytical results presented in this review can be applied for not only 
presenting but also understanding all the details of calculating the decay widths of the neutral Higgs 
boson discovered at the LHC by allowing for nontrivial pseudoscalar as well as scalar couplings of the 
Higgs boson to fermions pairs. In the SM limit, it turned out that our numerical results solely based 
on the analytic expressions and supplemental materials provided in this review show good consistency 
with those presented by, for example, the LHC Higgs Cross Section Working Group implementing the 
state-of-the-art theoretical calculations [84]. We expect that the analytic expressions for LO decay 
widths and QCD corrections together with the SM ELW corrections presented in this review could be 
used for analyzing the decays of neutral Higgs bosons appearing in BSM models with the precision 
comparable to that achieved through the full-fledged theoretical calculations. 


The second application of our analytical results was to the neutral 2HDM Higgs sector in the 
presence of nontrivial CP-violating phases of the complex quartic couplings. Specifically, we have taken 
the scenario where the lightest Higgs boson is purely CP even while the two heavier Higgs bosons are CP- 
mixed states exhibiting maximal CP violation when they are degenerate. To illustrate the typical decay 
patterns of the CP-violating heavy neutral Higgs bosons, we have taken the type-I 2HDM and we have 
contrasted them with those of purely CP-even or CP-odd neutral heavy Higgs bosons. Incidentally, we 
also have presented the decay pattern of a charged Higgs boson appearing in the framework of 2HDMs. 


Thirdly, we have presented the constraints on the couplings of the Higgs boson weighing about 125 
GeV obtained from implementing the global fits to all the Tevatron and LHC data available up to the 
summer in 2018. The global fits were based on several scenarios of the couplings and form factors for 
the main Higgs production and decay modes. Generally, the constraints turned out to be already tight 
on the tree-level couplings with the possibility that the loop-induced couplings could deviate sizably 
from the SM predictions. 


Finally, we emphasize that, even with this comprehensive review, there remain lots of improvements 
for the QCD and electroweak corrections to be made so as to match the expected ever-increasing 
precision measurements from upgraded LHC and future high energy collider experiments. Furthermore, 
to draw a more unified picture for the genuine BSM physics, many BSM scenarios considered seriously 
at present and expected to be developed fully in a concrete form in the near future have to be included. 
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Appendices 


This section consists of six appendices. Appendix[A]is for a summary of the SM parameters used for the 
numerical estimates of the Higgs decay widths and a description of the running of the strong coupling 
constant and quark masses. Appendix |B] is for the supersymmetric contributions to the loop-induced 
couplings of the Higgs boson to two gluons, two photons and Zy. Appendix|[Clis for the QCD corrections 
to the partial width of the Higgs-boson decay to two photons. We work out the relations among the 
parameters of the most general 2HDM in Appendix D] and we apply them for deriving cubic Higgs- 
boson self-couplings in Appendix [E] Finally, in Appendix [F] we briefly introduce a few well-developed 
numerical packages for calculating precise SM and full BSM-dependent ELW corrections. 


A Standard Model Parameters 


In this appendix, we summarize the SM parameters used for the estimation of decay widths of Higgs 
bosons. And we also show the running of the strong coupling constant and quark masses. 


A.1 Input parameters 
‘The SM parameters used for the estimation of decay widths of Higgs bosons are [6] [84]: | 
e Gauge coupling strengths 


a.(Mz) = 0.118 0.0015, 
a(0) = 1/137.035999 , 
a(My) = 1/128. (A.1) 


e Electroweak parameters 
Mw (80.385 + 0.015) GeV, Tw = (2.085 + 0.042) GeV; 
Mz = (91.1876 + 0.0021) GeV, Tz = (2.4952 + 0.0023) GeV; 
Grp = 1.1663787(6) x 1075 GeV-?. (A.2) 


The vev v of the SM Higgs field is given by v = (/2Gz) ~ 246.22 GeV with Gp = V29? /8M3, 
and v = 2Mw/g. For the square of the sine of the weak mixing angle, we adopt the so-called 
on-shell scheme in which the tree-level relation sj, = 1 — Mj,/M? is promoted to define the 
renormalized sj, to all orders in perturbation theory [293]: s% = 0.22290, cy = 1—s%, = 0.77710, 
g(Mw) = e/sw œ 0.66366 and g/(My) = e/cw œ 0.35544 with e = e(My) = 2,/ra(My) = 
0.31333. 


e Lepton masses 


M, = (105.6583715 + 0.00000035) MeV, M, = (1776.82 + 0.16) MeV. (A.3) 


e Quark masses 
M; = (172.51) GeV, 
my(™m) = (4.18+ 0.03) GeV, 
Me(3 GeV) = (0.986 + 0.026) GeV. (A.4) 
Note that the pole mass is used for t quark while, for b and c quarks, MS masses are used. 


40For the precision measurement of the Fermi constant Gp, see Refs. [115] : 
“1The electron and u, d-quark masses are not included as their masses are too tiny to influence the numerical analyses 
made in this work. 
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Table 12: Running of as(u), Melu), ma (u), and ™m,(u). MZ and M,_ are introduced for decoupling 
effects from matching the (effective) theory with Np — 1 with the (full) theory with Np at the scale M,. 


p [GeV] as(H) Melu) [GeV] molu) [GeV] mulu) [GeV] 
1.52+ (M7) | 3.48 x 107} 1.19 x 10° 5.40 x 10° 3.21 x 10? 
3 2.54 x 107! 986x107! 4.47 x 10° 2.66 x 10? 

4.18 (mT) | 2.25x10-! 9.21x 107! 4.18 x 10° 2.48 x 10? 
4.93_ (M,_) | 2.13 x 107} 894x107! 4.06 x 10° 2.41 x 10? 
4.93+ (M7) | 2.14 x 107! 8.92x107! 4.06 x 10° 2.41 x 10? 
10 1.78 x 107} 8.00x 107! 3.64 x 10° 2.16 x 10? 

20 1.53 x 107} 7.32x 107! 3.33 x 10° 1.98 x 102 

30 1.42 x 107} 7.00x 107! 3.18 x 10° 1.89 x 10? 

40 1.35 x 107} 6.79x107! 3.09 x 10° 1.84 x 10? 

50 1.380 x 107} 6.65x 107! 3.02 x 10° 1.80 x 10? 

60 1.26 x 107! 654x107! 2.97 x 10° 1.77 x 10? 

70 1.23 x 107! 644x107! 2.93 x 10° 1.74 x 10? 

80 1.20 x 107} 637x107! 2.90 x 10° 1.72 x 10? 

90 1.18 x 107} 630x107! 2.87 x 10° 1.70 x 10? 

100 1.16 x 107} 625x107! 2.84 x 10° 1.69 x 10? 

110 1.15 x 107} 620x107! 2.82 x 10° 1.67 x 10? 

120 1.13 x 107} 615x107! 2.80 x 10° 1.66 x 10? 
125.5 1.13 x 107} 613x107! 2.79 x 10° 1.66 x 10? 
130 1.12 x 107} 611x107! 2.78 x 10° 1.65 x 10? 

140 1.11 x 107} 6.08 x 107} 2.76 x 10° 1.64 x 10? 

150 1.10 x 107! 604x107! 2.75 x 10° 1.63 x 10? 

160 1.09 x 107! 6.01x107! 2.73 x 10° 1.62 x 10? 

170 1.08 x 107} 5.98x107! 2.72 x 10° 1.62 x 10? 
172.5_(M,_) | 1.08 x 107} 5.98 x 107} 2.72 x 10° 1.62 x 10? 
172.57 (M7) | 1.08 x 107} 5.98 x 107}! 2.72 x 10° 1.61 x 10? 
180 1.07 x 107} 5.96 x 107! 2.71 x 10° 1.61 x 10? 

190 1.06 x 107} 5.93 x 107! 2.70 x 10° 1.60 x 10? 
200 1.06 x 107! 5.91 x107! 2.69 x 10° 1.60 x 10? 
300 1.01 x 107! 5.74x 107! 2.61 x10? 1.55 x 10? 
400 9.76 x 107? 5.62 x 107! 2.56 x 10° 1.52 x 10? 
500 9.52 x 1072? 554x107! 2.52 x 10° 1.50 x 102 
1000 8.85 x 1072 5.29x 107! 2.41 x 10° 1.43 x 10? 


A.2 Running of the strong coupling constant and quark masses 


We neglect the running of the SU(2)z and U(1)y electroweak couplings and the leptons masses, measured 
experimentally with great precision. For the running of the strong coupling strength as(u) and the MS 
quark masses 77,(j1), we use the most recent version of RunDec [294] [295] in which five-loop corrections 
of the QCD beta function and four-loop decoupling effects are included. The results are shown in Fig. 
and Table We note that a,(125.5GeV) = 0.1126 and m,(M;) = 161.5 GeV. For the pole mass of 
b quark, the three-loop conversion relation is taken to give M, = 4.93 GeV [6]. On the other hand, 
for the pole mass of c quark, we take the relation between the on-shell charm-quark and bottom-quark 
masses, giving Me = M, — 3.41 GeV = 1.52 GeV [6] [84] [296]. 
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Figure 36: Running of the strong coupling constant a,() (upper) and the MS quark masses 7™7;(j1) 
(middle) and Mmg (u) (lower). In the upper frame, the vertical line locates the position of My = 125.5 
GeV. In the middle and lower frames, we also show m;,(j) used in the calculations of '(H — yy) with 
the expressions given in Eq. (A.5). The vertical lines locate the positions of the pole masses of M, and 
M,,-. And, in the lower frame, the input values for b and c quark masses are denoted by bullets, see 
Eq. (A.4). The open circles in the middle and lower frames denote the positions m;(u = M) = M: and 
M(u = Mi) = My, respectively, with M, = 172.5 GeV and M, = 4.93 GeV. 
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Specifically for the loop-induced decay H — yy, we use the running masses m;»(u) of which the 
expressions are given at the four-loop level by [148] 297] 


ee | aii ye 1+ [ew] + cf [asta] + cf Ea ne 
me" Lae) l+cef +c EJ +c [eca | 
where the numerical values of the six dimensionless coefficients are 
œ = 1.17549, œ& = 1.50071, c& = 0.172478; 
c = 1.39796, c = 1.79348, cy = —0.683433. (A.6) 


Note that m,(M,) = M; as denoted by open circles in the middle and lower frames in Fig. 


B Supersymmetric Contributions to the Hgg, Hyy, and HZ7y 
Form Factors 


In this appendix, we present the contributions to the loop-induced Hgg, Hy y, and H Zy form factors 
from the triangle diagrams in which charginos, charged and/or colored sfermions, and/or charged Higgs 
bosons are running. 


In the minimal supersymmetric extension of the SM (MSSM), the form factors of AS%7 and AP97 
denoting new MSSM contributions to the Hgg, Hyy vertices are given by: 


v2 


g ae z 
A == FHF h Ape Fo(t,3,) > 
=h, To bi d fj 
AP? = 0; (B.1) 
2 

U 

AS? = V2qg a UA Tf wa Fost) = `> NE Q InFS, Z EO) 

Í=X1 Xa fit t2,b1 ,b2,71,72 fi 
y2 
TO naa zag) 
AP? = V29 5 Ont — wa Fos Ti) (B.2) 


=X X2 


where Tjz = M}, /4m2 with i = 1,2,3 denoting three neutral Higgs bosons and N = 3 for (s)quarks 
and NÅ = 1 for status, respectively. The form factor Fo(r) is given by 


Fo(r) =r! [—1+ tT 'f(r)| , (B.3) 


which takes the value of 1/3 when 7 = 0, see Fig. 


On the other hand, the form factors of AS77 and AP% denoting new MSSM contributions to the 
H Zy vertices may take forms of 


AS" = =y] 


S 
Mgt M+, M x zty- JH.s+y- 
o . ef ( Xj pie M ) UBS x; IHE 


Cw Sw 


NeQs 
AR, Eons RIRI V ONG Mg ma) 
j,k 


stb, WSW 
2 
gJH;H+H- V 
-2 z Ty (tint, Ag+) , 
Cw Sw Mý- 
Z ig P 
KPA = = 2 v myg | Met, Mek, Met ) Vzgr he B.4 
i V2 ; g Xj , XE? Xk 2X; Xk IBEX ? ( ) 


where T;y+ = Mẹ, /4Mẹ}+ with i = 1,2,3, Ag+ = M3/4M}+, and N£ = 3 for squarks and N = 1 for 
staus, respectively. For the explicit form of the three loop functions of f(m, M2, M2), g(mMmı, M2, M2), 
and Co(m1, M2, M2), we refer to [227]. Note that they implicitly depend on M}, and M2. For the Higgs 
couplings to SUSY particles, see subsubsection [2.4.3] and the relevant Z-boson interactions are given by 
the following Lagrangian terms: 


e Z-sfermion-sfermion 


where gz = e/(Swew) and 
Izpi = KUL UL, — Q ssi, (B.6) 
with I)” = +1/2 and 12° = -1/2. 


e Z-chargino-chargino 


Laxty- = -gz Xi 7" (aaa = aziz) Xj Zu (B.7) 
where the vector and axial—vector couplings are given by 
tags = F Cda (Cr) + (Cala (Cada) - yds. 
azgi = F (Cra (Cr) — (Cr) (Ca) - (B8) 
For completeness, we recall that the Z-boson couplings to the quarks and leptons are given by 
Leary =- GZ fy" (uzri = Az ffs) f Za, (B.9) 


with ug7, = If J2 — Qysiy and AgFf = If /2 in terms of the third component of weak isospin E 
and the electric charge Qs of each fermion f. 


C QCD Corrections to [(H — yy): Css(T) and Cps(T) 


The scaling factors Css(T) and C,;(7) for the QCD corrections to the decay width of a Higgs boson H 
into two photons might be given by [183] 


4 4 
Costs p) = CH (r) + OF (1) og silos a . Gto- as tog dios a (Cl) 
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where r = Mj,/4m7(uq) with the renormalization scale yg = My/p. As demonstrated in Section 
we take p = 2. We note again the running mass mg is normalized as m,(M,) = Mọ- 


The C¥ (r) function is given via the following relation 


200 +0+0 +6") 
3(1—6)° 


—64 Li3(—0) ln 0 + 14 Lig(6) In? 0 + 8 Lig(—8) In? 0 + + In* 0 


F,s(7)C# (T) [108 Li4(0) + 144 Lia(—0) — 64 Lis (0) Ind 


+G n9 + 1663 Ind +186 


20(1 +0) . } 
——_——— | —32 Li3(—@) + 16 Liz(—0) ln 0 — 4&2 ln 0 
sQ gr 32M6) + 16Lis(—8) nd -4G 
—2 2 1 —2 2 
807 aaa es 60 (3 ox giao 
3 (1—0) 3 (1—0) 
= 2 — 7 92 3 
40 (5 r ) tn(1 — 6) In26 + 20 (3 + 256 gi +30?) in? 4 
3 (1—0) 9(1-86) 
01—140 +0? 0? 80 (1+0 400 
ell os lad 5 ame e s (C.2) 
3(1- 6) (1 — 0) (1 — 0) 3 (1—0) 
where @ is a T-dependent function defined by 
yil=r =l 
psie E (C.3) 
Jier Fl 
The three values, Çz, Cs and C4, of the Riemann’s zeta function are given by 
T? Tî 
== = 1.202 = — C.4 
Go 6 ’ E 0 05690, C4 90 ’ ( ) 
and the polylogarithm function is defined by a power series in a complex variable z as follows 
Li,(z) = — , R= 1,2,3,>-* C.5 
=m (C5) 


For analytic continuation to the complex 7 plane, the replacement T — T + Oi is understood. The 
C#(r) function is given via [213] 


Fist) CO; a) soe" i +(r-—2)fíTt)-(Tr-1)r eal , (C.6) 
with the function f(T) defined in Eq. (128). 


On the other hand, the Cf(r) and C3(r) functions are given via the following relations [213 


(1+ 6") : 128p. 
Goto) {2 Li4(0) + 96 Li4(—0) =—— | Lis(0) + Lis(—0)| Ind 


Farer) = 25 3 


28. 16. 1 
ie Lip(0) In? 6 + = Lip(—6) In? 6 + T In‘ 


42Rqs. (10) and (12) in Ref. contain typos, see footnote 3 of Ref. for locating them. 
431t is not yet clear whether ¢3 is given in a compact form or not, unlike C24. 
44 Lin (z) has a branch cut discontinuity in the complex z plane running from 1 to oo. 
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8 32 
+362 In? 0 + 3 $3 ln + 1264} 


a a -* Lis(6) — “ Lis(—0) + 16 Liz(0) nd 
+ Lip(—0) n6 + = In(1 — 6) In? — C2 Ind + 6s 
aoa T (C.7) 
Fcio = 2r [re -r LO], (C8) 


respectively. 


D Input parameters for the most general 2HDM potential 
In this appendix, we work out the relations among the parameters needed to fully specify the most 
general 2HDM potential, the masses of neutral and charged Higgs bosons, and the mixing matrix O. 


In subsection [2.3] we demonstrate that one needs all the elements of the following set of parameters 
L= {v tg ; Imal ; Ai ; A2 ; A3 ; AA ; làs] ; BY ; |A7| , Q5 F 2E , De F E, Q7 F E ; sign|cos(¢12 + E)]}} ’ (D.1) 


to fully specify the most general 2HDM scalar potential, see Eq. (37). The set Z contains 13 parameters 
plus 1 sign with sin(¢;2 + €) being determined by the third CP-odd tadpole condition in Eq. (35). 
Alternatively to the set Z, one may use the following equivalent set 


T! = {v,ty, Re(m?c%); (D.2) 

Ar, À2 , A3 , Aa, Re(Ase”*) , Re(Age) , Re(Aze*) , Sm(Ase%) , Sm (Age) , Sm(Aze*)} . 
The above set Z’ contains 10 parameters for the real and complex quartic couplings ;_7 and any 7 
of them, in principle, can be traded with the 4 masses of charged and neutral Higgs bosons and the 3 
independent angles of the 3 x 3 orthogonal mixing matrix O by judiciously exploiting Eq. and the 


matrix relation OT M§O = diag( M}, , M} , Młp). || By choosing Sm(Ase”*), Re(Age’), and Re(Aze"*) 
as independent input parameters, one may use the following set of more physical parameters: 


P = {v,tg, Relm? c); Mz, , Mu, , My, , Mpe , {0} ; Sm(Ase”*) , Re(Age’®) , Re(Aves)}.  (D.3) 
B 3 


Explicitly, we find that the 7 quartic couplings of \1, A2, Ag, Aa, Re(Ase?4), Sm(Age“’), and Sm(Aze") 
in the set Z’ can be expressed in terms of My,,,, My+, and the elements of the mixing matrix O in 
the set P as follows: 


Sp 2 Ai O51 2 O59 2 O33 2 3tg ig t3 k 
À = wa Re(mj,e">) + Ma, + Mn oe Mn, — 7 Relee y+ J Bere ), 


3 2v2c3 2v2c3 2v2c3 
2 2 2 
Cp 2 iê Oi 2 Obn2 2 O43 2 l i£ 3 it 
ào = —— R M M5, 4 M;,, 4 Re(A — —RFe(A ; 
í 2us% ele 2u?s% Hy + 2u?s% Ho 20253 Hs At? elie”) At, ene”) 
1 2 On Ogi O$,20 02 04130403 
ee R 2 „if * M2 pı 2 M2 1 2 M2 1 2 M? 
USC ar eg ye H $ v2sgeg  ™ i U?S6Cg vsecg ™ 


45For M2, see Eq. and the two subsequent relations following it. 
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To 

àg = a Re(m?e) z ZM + 7 Mj, + 2 Mj, + ee Mi, -= Ae) = E eare") 
Re( ae) = =i Reloige®) — SEL Mẹ, — O22 Mp, — SE Mp, = T-Re( Ages) — Peoro“), 

Sm (Aso) = -TE Mp, = R Mhp ~ SO Mẹ, — tolse) 

Sm(Aze%*) = a Mi, — were Mi, — sence Mi, — (Ase) . (D.4) 


We find that our results are consistent with those presented in Refs. [298] [299]. We note that Az can be 
obtained from A; or vice versa by exchanging cg + sg, ¢1 + 2, and Ag + àz. The same observation 
could be applied for Sm(Age“*) and Sm(A ze") which are vanishing when each Higgs boson is purely 
CP-even or CP-odd state and Sm(Ase”5) = 0. About àz and Ay, we note that 3 + Aq is independent 
of Re(mize) and Mz., the neutral Higgs mass contributions to 3 (Ay) are involved with only the 
CP-even (CP-odd) components of each Higgs boson, and the contributions from Re(Age*) and Re(A7e") 
are in common. In passing, we check that Eq. for the difference between \4/2 and Re(Ase”") is 
satisfied by noting the relation O7, M}, + O7,My, + O23 Mh, = M3. 


E Cubic Higgs-boson self-couplings in 2HDMs 


In this appendix, we apply the relations among the 2HDM input parameters obtained in Appendix [D] 
to derive cubic Higgs-boson self-couplings when the lightest Higgs boson is purely CP even as assumed 
in subsubsection by taking the following O matrix: 


—Sa Caw CaSw 
O = Ca S46, Sasu . (E.1) 
0 -s,; Gy 


In this case, using Eq. (55), Eq. for the cubic self-couplings of the Higgs weak eigenstates, and 
Eq. (D.4) for the conversion relations, we find the couplings of the heavy neural Higgs bosons Hə to 
a pair of the lightest Higgs bosons are given by 


C Re(m7,e") 
Jmm OT z |4s28cg-a + 60288 8-aC3-0 > 652904 a| ( = 
538 v 
c M3 1 M2. 
E |s28c8-a + 2C285p-aCh—a = 252803_a ( 2 ) ; 
2 2 2 
SCA _ ; 3CuS8—aCh_ ; 3CuS8—-alh_ , 
Ea Ym (Azer) — So Re(Age’®) + = Rezet), 
S28 4s 4c 
B B 
s Rte(m?,e% ) 
ae -= |4s28c8-a + 60265p-aC3_a — 652904 a| ( ie 
S38 v 
2 2 
Sw 2 3 Mp 1 (Ma, 
= |s29¢9-0 + 202856-aCg_o, — 282904 a ( J2 ) 3 \ 2 
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Re(Azves).  (E.2) 


We note that gu, m m, [Co = Im m m /Su When Sm(Ase**) = 0 and My, = Mu, are taken and the 
contributions proportional to the input quartic couplings are suppressed by the factor of ae We 
further note that the couplings gp, m m, and Jy, m, m, ave vanishing in the limit where gy vy = 9i,vv = 9 
or Cg-q = 0. Otherwise, they are non-vanishing. To be more specific, as in subsubsection|3.7.2} we take 


Mu, \? Mu, 
dy = Vép (a) ; d3 = E0 & 
2 3 


together with Eq. (E.1). These specific choices of 523 fully fix all the elements of the 3 x 3 orthogonal 
mixing matrix O in terms of the masses of the neutral Higgs bosons and the parameter €o through 


(E.3) 


Sa = —Vl—ecg+Ves83, Co=V1—€S8g+ Vee, 
Ĉu = Mir y Sym Mir (E.4) 
Mh, + Ma, Mz, + Ma, 
where i 1 
C= Gi + ve) (E.5) 


Incidentally, we note sg, = /l—€e and cg_,g = Ve. Taking ceo = 0.05, My, = 125.5 GeV, and 
My, ~ Mzy,, we find the couplings of the heavy neural Higgs bosons H23 to a pair of the lightest Higgs 
bosons are given by 


M? 
keeping the two most significant contributions: the first term comes from —s,,cg_a 5> and the second 
Retna = 0.1 as 
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taken in subsubsection [3.7.2] especially for |Re(mi,c")| /s23 ~ M. Incidentally, we find 
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Re(Aze’s) , (E.7) 


which becomes the SM coupling of Mj, /2v? when cg_ = 0, see Eq. (6). 


Finally, we address the contributions from the charged-Higgs-boson loops to the decay processes of 


neutral Higgs bosons into two photons in the 2HDM and the MSSM which are mentioned in subsubsec- 


tion]3.7.2|and subsection [5.3] Using Eq. (55), Eq. (58), and Eq. (D.4), we obtain the following couplings 


46Note we are taking cu > 0 and s,, > 0. 
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of neutral Higgs bosons to a pair of charged Higgs bosons: 


M? + 4 Rte(m?,e% ) 
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S6 Cp 


In the 2HDM as well as in the MSSM, the contributions from the charged-Higgs-boson loops to the 
neutral Higgs boson decays into two photons enter through the form factor £] 


AS?(H*) = vn, (Mm E.9 
a ) = Imut- IMP, 0 4M? ” ( 3 ) 


In the infinite charged-Higgs-boson mass limit, using Fo(0) = 1/3, we find 


ASI (HF) = s6-a/3, ASZ (HF) = cucg-a/3, ASI(HF) = SwCg-a/3. (E.10) 


We note the non-decoupling feature of AST (H=) 300| which does not vanish even when all 
the heavier Higgs bosons including the charged one are decoupled. On the other hand, AS} 3(H =) are 
vanishing when the heavy neutral Higgs bosons are decoupled or when sg-a — 1 and cg-a — 0. 


F Packages for electroweak corrections 


In this appendix, we introduce various numerical packages for calculating precise SM and full BSM- 
dependent ELW corrections. 


To begin with, in Table we list the two packages of HDECAY and PROPHECY4F 
for precise SM ELW corrections. We note that HDECAY is commonly used except 
for the decay mode H + VV — 4f for which PROPHECY4F provides the complete O(a) electroweak 
corrections to the Higgs decays into four fermions through intermediate W and Z bosons. And, as 
far as we have configured, there exist no reliable theoretical calculations and/or numerical packages 
implementing the ELW corrections to the radiative H —> Zy decay. Incidentally, for the other radiative 
decays of H — gg and H — yy, we refer to Refs. for theoretical calculations. 


See, Eq. (Bp. 
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Table 13: Two packages of HDECAY and PROPHECY4F for the SM ELW corrections. In each package 
line, decay modes are checked when the SM ELW corrections to them could be calculated by use of it. 
Package \Decay Mode | H > ff | H>WW/ZZ->4f | Hogg) Hoy) HO Zy 
HDECAY v v v v 

NLO Approx. NLO NLO NLO 
V 
PROPHECY4F NLO 


Table 14: Packages for BSM-dependent ELW corrections. In each package line, BSM models are 
checked when ELW corrections in those models are implemented in the package. The last column is to 
note that PROPHECY4F is exclusively for the H > VV — 4f decays and, for other packages, to indicate 


the level of precision at which the ELW corrections are calculated. 


Package \BSM Model | cxSM | 2HDM | MSSM | NMSSM Remark 
PROPHECY4F J J H > WW/ZZ > 4f 
2HDECAY V Full one-loop 
H-COUP V NLO 
HFOLD V Full one-loop 
FeynHiggs % Full one-loop 
NMSSMCALCEW V Full one-loop 
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In Table we list various packages for ELW corrections in the BSM models considered in this 
review. Additionally, as an example of BSM models containing Higgs sectors beyond the cxSM and 
2HDM/MSSM, the Next-to-MSSM (NMSSM) [70] is included. First of all, exclusively for the 
H — VV — 4f processes, PROPHECY4F can be used for the ELW corrections in the cxSM and 2HDM. 
In passing, we note that the full NLO corrections in the real singlet extension of the SM (rxSM) 
and the SM4G (SM with the fourth generation) are also implemented 
in PROPHECY4F. 


For the calculations of the ELW corrections in 2HDM, the packages 2HDECAY and H-COUP 
310) can be used. We comment that the package H-COUP can be used for the NLO ELW corrections 
in the rxSM and inert doublet models as well. For the MSSM ELW corrections at the full one-loop 
level precision, one can use the packages HFOLD [311] and FeynHiggs [812]. On the other hand, in the 
NMSSM, the package NMSSMCALCEW [313] provides the full one-loop ELW corrections. 


Before closing this appendix, we introduce two useful webpages containing significantly extensive 
lists of numerical packages or tools for various purposes from model building to event generation, etc: 


1) Supersymmetry Les Houches Accord : http://skands.physics.monash.edu/slha/~, 


2) HEPForge Projects : https://www.hepforge.org/projects . 


References 


[1] G. Aad et al. [ATLAS], “Observation of a new particle in the search for the Standard 
Model Higgs boson with the ATLAS detector at the LHC,” Phys. Lett. B 716 (2012), 1-29 


doi:10.1016/j.physletb.2012.08.020 |arXiv:1207.7214  [hep-ex]]. 

[2] S. Chatrchyan et al. [CMS], “Observation of a New Boson at a Mass of 125 GeV with the CMS 
Experiment at the LHC,” Phys. Lett. B 716 (2012), 30-61 doi:10.1016/j.physletb.2012.08.021 
arXiv: 1207.7235) [hep-ex]]. 

3] S. L. Glashow, “Partial Symmetries of Weak Interactions,” Nucl. Phys. 22 (1961), 579-588 

doi:10.1016/0029-5582(61)90469-2. 

4] S. Weinberg, “A Model of Leptons,” Phys. Rev. Lett. 19 (1967), 1264-1266 

doi:10.1103/PhysRevLett.19.1264. 

5] A. Salam, “Weak and Electromagnetic Interactions,” Conf. Proc. C 680519 (1968), 367-377 

doi:10.1142/9789812795915_0034. 

6] P.A. Zyla et al. (Particle Data Group), to be published in Prog. Theor. Exp. Phys. 2020, 083C01 

(2020). 

7| G. Aad et al. [ATLAS], “Combined measurements of Higgs boson production and decay using up 
to 80 fb of proton-proton collision data at \/s = 13 TeV collected with the ATLAS experiment,” 


Phys. Rev. D 101 (2020) no.1, 012002 doi:10.1103/PhysRevD.101.012002 |arXiv:1909.02845) [hep- 


ex]]. 


[8] A. M. Sirunyan et al. [CMS], “Combined measurements of Higgs boson couplings in proton-proton 
collisions at \/s = 13 TeV,” Eur. Phys. J. C 79 (2019) no.5, 421 doi:10.1140/epjc/s10052-019-6909- 


y jarXiv:1809.10733, [hep-ex]]. 


[9] G. Aad et al. [ATLAS and CMS], “Combined Measurement of the Higgs Boson Mass in pp Collisions 
at \/s = 7 and 8 TeV with the ATLAS and CMS Experiments,” Phys. Rev. Lett. 114 (2015), 191803 


doi:10.1103/PhysRevLett.114.191803 |arXiv:1503.07589 [hep-ex]]. 


84 


[10] A. M. Sirunyan et al. [CMS], “A measurement of the Higgs boson mass in the diphoton decay 
channel,” Phys. Lett. B 805 (2020), 135425 doi:10.1016/j.physletb.2020.135425 |arXiv:2002.06398 
[hep-ex]]. 

[11] G. Aad et al. [ATLAS], “Evidence for the Higgs-boson Yukawa coupling to tau leptons with the 
ATLAS detector,” JHEP 04 (2015), 117 doi:10.1007/JHEP04(2015)117 jarXiv:1501.04943) [hep- 
ex]]. 

|12] A. M. Sirunyan et al. [CMS], “Observation of the Higgs boson decay to a pair of 7 leptons 
with the CMS detector,” Phys. Lett. B 779 (2018), 283-316 doi:10.1016/j.physletb.2018.02.004 


arXiv:1708.00373 [hep-ex]]. 


13] A. M. Sirunyan et al. [CMS], “Observation of Higgs boson decay to bottom quarks,” Phys. Rev. 


Lett. 121 (2018) no.12, 121801 doi:10.1103/PhysRevLett.121.121801 |arXiv:1808.08242) [hep-ex]]. 


14] M. Aaboud et al. [ATLAS], “Search for Higgs bosons produced via vector-boson fusion and decaying 
into bottom quark pairs in ys = 13 TeV pp collisions with the ATLAS detector,” Phys. Rev. D 


98 (2018) no.5, 052003 doi:10.1103/PhysRevD.98.052003 |arXiv:1807.08639) [hep-ex]]. 


15] M. Aaboud et al. [ATLAS], “Observation of Higgs boson production in association with a 
top quark pair at the LHC with the ATLAS detector,” Phys. Lett. B 784 (2018), 173-191 


doi:10.1016/j.physletb.2018.07.035 |arXiv:1806.00425) [hep-ex]]. 
16] G. Aad et al. [ATLAS], “A search for the dimuon decay of the Standard Model Higgs boson with 
the ATLAS detector,” |arXiv:2007.07830) [hep-ex]]. 


17] A. M. Sirunyan et al. [CMS], “Evidence for Higgs boson decay to a pair of muons,” 


arXiv:2009.04363 [hep-ex]]. 
18] E. Gildener and S. Weinberg, “Symmetry Breaking and Scalar Bosons,” Phys. Rev. D 13 (1976), 
3333 doi:10.1103/PhysRevD.13.3333. 


19] S. Weinberg, “Implications of Dynamical Symmetry Breaking,” Phys. Rev. D 13 (1976), 974-996 
doi:10.1103/PhysRevD.19.1277. 


20| L. Susskind, “Dynamics of Spontaneous Symmetry Breaking in the Weinberg-Salam Theory,” Phys. 
Rev. D 20 (1979), 2619-2625 doi:10.1103/PhysRevD.20.2619. 


21) Y. A. Golfand and E. P. Likhtman, “Extension of the Algebra of Poincare Group Generators and 
Violation of P Invariance,” JETP Lett. 13 (1971), 323-326. 


22| D. V. Volkov and V. P. Akulov, “Is the Neutrino a Goldstone Particle?,” Phys. Lett. B 46 (1973), 
109-110 doi:10.1016/0370-2693(73)90490-5. 


23| J. Wess and B. Zumino, “Supergauge Transformations in Four-Dimensions,” Nucl. Phys. B 70 
(1974), 39-50 doi:10.1016/0550-3213(74)90355-1. 


24) N. Arkani-Hamed, A. G. Cohen and H. Georgi, “Electroweak symmetry breaking from dimensional 
deconstruction,” Phys. Lett. B 513 (2001), 232-240 doi:10.1016/S0370-2693(01)00741-9 
ph/0105239 [hep-ph]. 

[25] N. Arkani-Hamed, A. G. Cohen, T. Gregoire and J. G. Wacker, “Phenomenology of electroweak 
symmetry breaking from theory space,” JHEP 08 (2002), 020 doi:10.1088/1126-6708 /2002/08/020 

arXiv:hep-ph/0202089) [hep-ph]. 

[26] N. Arkani-Hamed, A. G. Cohen, E. Katz, A. E. Nelson, T. Gregoire and J. G. Wacker, “The Mini- 
mal moose for a little Higgs,” JHEP 08 (2002), 021 doi:10.1088/1126-6708/2002/08/021 
ph/0206020 [hep-ph]. 

[27] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, “The Hierarchy problem and new dimensions at 

a millimeter,” Phys. Lett. B 429 (1998), 263-272 doi:10.1016/S0370-2693(98)00466-3 


ph/9803315 {hep-ph]]. 


85 


[28] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, “New dimensions at a millimeter 
to a Fermi and superstrings at a TeV,” Phys. Lett. B 436 (1998), 257-263 doi:10.1016/S0370- 


aise Ci vases 0 ) arXiv: hep- Eee aai 


[29] L. Randal | R. Sundrum, Large mass hierarchy from a small extra dimension,” Phys. Rev. 
Lett. 83 3 (1999), 3370- 3373 doi: 10. 0.1108/PhysRevLett. 83.3370 |arXiv:hep-ph/9905221| [hep-ph]. 


ion,” Phys. Rev. Lett. 83 (1999), 


4690-4693 doi:10.1103 JENN 83.4690 farXiv:hep-th/9906064 nap: SEEN! [hep-th]. 
[31] G. Jungman, M. Kamionkowski and K. Griest, “Supersymmetric dark matter,” Phys. Rept. 267 
Se) 195-373 doi; 10. ILG AOAI 1573(95)00058 5 jarXiv: OE ae jep- al 


Phys. TE 405 (2008), 279.390 doi:10. 1016 /j. physrep. 3004.08.31 arXiv:hep-ph/0404175 Ta arXiv-hep-ph/0404175 hep- 
phj]. 


76 (2003), 1 doi:10. 1103/RevModPhys, 76.1 [arXiv:hep-ph/0303065 a arXiv:hep-ph/0303065 Thep- phl]. i 
[34] D. E. Morrissey and M. J. Ramsey-Musolf, “Electroweak baryogenesis,” New J. Phys. 14 (2012), 
125003 doi:10.1088/1367- 2630/14/12/125003 |arXiv:1206.2942 arXiv:1206.2942 1206.2942 er phh. 


[35] W. Buchmuller, R. D. Peccei and T. Yanagida, “Lej asis as the origin of __ _ Ann. 
Rev. Nucl. Part. Sci. 55 (2005), 311-355 doi:10. 1146 /annurev. nucl.55.090704.151558 
[hep-ph]. 

[36] M. C. Gonzalez-Garcia and M. Maltoni, “Phenomenology with Massive Neutrinos,” Phys. Rept. 
460 (2008), 1-129 doi:10.1016/j.physrep.2007.12.004 [hep-ph]. 

[37] D. H. Lyth and A. Riotto, “Particle physics models of inflation and the cosmological density pertur- 
bation,” Phys. Rept. 314 (1999), 1-146 doi:10.1016/S0370-1573(98)00128-8 
hep-ph]. 

[38] V. Shiltsev and F. Zimmermann, “Modern and Future Colliders,” Rev. Mod. Phys. 93 (2021) no.1, 

15006 doi:10.1103/RevModPhys.93.015006 [physics.ace-ph]]. 

[39] H. M. Gray, “Future colliders for the high-energy frontier,” Rev. Phys. 6 (2021), 100053 

doi:10.1016/j.revip.2021.100053. 

[40] V. Barger, P. Langacker, M. McCaskey, M. Ramsey-Musolf and G. Shaughnessy, “Com- 

plex Singlet Extension of the Standard Model,” Phys. Rev. D 79 (2009), 015018 

doi:10.1103/PhysRevD.79.015018 [hep-ph]]. 

[41] M. Gonderinger, H. Lim and M. J. Ramsey-Musolf, “Complex Scalar Singlet Dark Matter: Vacuum 
Stability and Phenomenology,” Phys. Rev. D 86 (2012), 043511 doi:10.1103/PhysRevD.86.043511 


arXiv:1202.1316) [hep-ph]. 


[42] R. Coimbra, M. O. Sampaio and R. Santos, “ScannerS: Constraining the phase diagram of a 
complex scalar singlet at the LHC,” Eur. Phys. J. C 73 (2013), 2428 doi:10.1140/epjc/s10052-013- 
2428-4 |arXiv:1301.2599 [hep-ph]. 

[43] M. Jiang, L. Bian, W. Huang and J. Shu, “Impact of a complex singlet: Electroweak baryoge- 
nesis and dark matter,” Phys. Rev. D 93 (2016) no.6, 065032 doi:10.1103/PhysRevD.93.065032 
arXiv: 1502.07574) [hep-ph]]. 

[44] F. Sannino and J. Virkajarvi, “First Order Electroweak Phase Transition from 
(Non)Conformal Extensions of the Standard Model,” Phys. Rev. D 92 (2015) no.4, 045015 
doi:10.1103/PhysRevD.92.045015 |arXiv:1505.05872) [hep-ph]]. 

[45] R. Costa, M. Miihlleitner, M. O. P. Sampaio and R. Santos, “Singlet Extensions of the Standard 
Model at LHC Run 2: Benchmarks and Comparison with the NMSSM,” JHEP 06 (2016), 034 


doi:10.1007/JHEP06(2016)034 |arXiv:1512.05355, [hep-ph]. 


86 


[46] 


[47] 


[48] 


59 


60 


61 


62 


63 


64 


65 


M. Mühlleitner, M. O. Sampaio, R. Santos and J. Wittbrodt, “Phenomenological Comparison 
of Models with Extended Higgs Sectors,’ JHEP 08 (2017), 132 doi:10.1007/JHEP08(2017)132 


arXiv:1703.07750) [hep-ph]]. 
C. Chiang, M. J. Ramsey-Musolf and E. Senaha, “Standard Model with a Complex Scalar Singlet: 
Cosmological Implications and Theoretical Considerations,” Phys. Rev. D 97 (2018) no.1, 015005 


doi:10.1103/PhysRevD.97.015005 [hep-ph]]. 

D. Azevedo, P. Ferreira, M. M. Mühlleitner, R. Santos and J. Wittbrodt, “Models with 
extended Higgs sectors at future ete~ colliders,” Phys. Rev. D 99 (2019) no.5, 055013 
doi:10.1103/PhysRevD.99.055013 [hep-ph]]. 

T. D. Lee, “A Theory of Spontaneous T Violation,” Phys. Rev. D 8 (1973), 1226-1239 
doi:10.1103/PhysRevD.8.1226. 


T. D. Lee, “CP Nonconservation and Spontaneous Symmetry Breaking,” Phys. Rept. 9 (1974), 
143-177 doi:10.1016/0370-1573(74)90020-9. 


| R. D. Peccei and H. R. Quinn, “CP Conservation in the Presence of Instantons,” Phys. Rev. Lett. 


38 (1977), 1440-1443 doi:10.1103/PhysRevLett.38.1440. 

P. Fayet, “A Gauge Theory of Weak and Electromagnetic Interactions with Spontaneous Parity 
Breaking,” Nucl. Phys. B 78 (1974), 14-28 doi:10.1016/0550-3213(74)90113-8. 

K. Inoue, A. Kakuto, H. Komatsu and S. Takeshita, “Low-Energy Parameters and Parti- 
cle Masses in a Supersymmetric Grand Unified Model,” Prog. Theor. Phys. 67 (1982), 1889 
doi:10.1143/PTP.67.1889. 

R. A. Flores and M. Sher, “Higgs Masses in the Standard, Multi-Higgs and Supersymmetric Mod- 
els,” Annals Phys. 148 (1983), 95 doi:10.1016/0003-4916(83)90331-7. 

J. F. Gunion and H. E. Haber, “Higgs Bosons in Supersymmetric Models. 1.,” Nucl. Phys. B 272 
(1986), 1 erratum: Nucl. Phys. B 402 (1993), 567-569] doi:10.1016/0550-3213(86)90340-8. 

F. J. Botella and J. P. Silva, “Jarlskog - like invariants for theories with scalars and fermions,” Phys. 
Rev. D 51 (1995), 3870-3875 doi:10.1103/PhysRevD.51.3870 [hep-ph]. 

G. C. Branco, L. Lavoura and J. P. Silva, “CP Violation,” Int. Ser. Monogr. Phys. 103 (1999), 
1-536, Chapters 22 and 23. 

M. Carena and H. E. Haber, “Higgs Boson Theory and Phenomenology,” Prog. Part. Nucl. Phys. 
50 (2003), 63-152 doi:10.1016/S0146-6410(02)00177-1 [hep-ph]]. 

G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher and J. P. Silva, 
“Theory and phenomenology of two-Higgs-doublet models,” Phys. Rept. 516 (2012), 1-102 
doi:10.1016/j.physrep.2012.02.002 [hep-ph]. 

P. Fayet and S. Ferrara, “Supersymmetry,” Phys. Rept. 32 (1977), 249-334 doi:10.1016/0370- 
1573(77)90066-7. 

H. P. Nilles, “Supersymmetry, Supergravity and Particle Physics,” Phys. Rept. 110 (1984), 1-162 
doi:10.1016/0370-1573(84)90008-5. 

H. E. Haber and G. L. Kane, “The Search for Supersymmetry: Probing Physics Beyond the 
Standard Model,” Phys. Rept. 117 (1985), 75-263 doi:10.1016/0370-1573(85)90051-1. 

R. Barbieri, “Looking Beyond the Standard Model: The Supersymmetric Option,” Riv. Nuovo 
Cim. 11N4 (1988), 1-45 doi:10.1007/BF02725953. 

S. P. Martin, “A Supersymmetry primer,” Adv. Ser. Direct. High Energy Phys. 18 (1998), 1-98 


doi:10.1142/9789812839657_0001 |arXiv:hep-ph/9709356) [hep-ph]. 


H. E. Haber, “The Status of the minimal supersymmetric standard model and beyond,” Nucl. Phys. 
B Proc. Suppl. 62 (1998), 469-484 doi: 10.1016 /S0920-5632(97)00688-9 |arXiv:hep-ph/9709450) [hep- 
phi). 


87 


Ply. Rept 407 (2005), 1-203 


Tr 10. smu a m 08.032 |arXiv:hep-ph/0312378 Xie ph/ 0012078 [hep-ph]. 


eee TE in high nergy TE 7 Publ hed i aoe USA: World 
Scientific Gara 555 p. 


edi in: Cambridge, UK: TETA Univ. Pr. ” (2006) 556 p. 

[69] P. Binetruy, “Supersymmetry: Theory, experiment and cosmology,” Published in: Oxford, UK: 
Oxford Univ. Pr. (2006) 520 p. 

[70] M. Maniatis, “The Next-to-Minimal Supersymmetric extension of the Standard Model reviewed,” 
Int. J. Mod. Phys. A 25 (2010), 3505-3602 doi:10.1142/S0217751X 10049827 [hep- 
phi]. 

[71] U. Ellwanger, C. Hugonie and A. M. Teixeira, “The Next-to-Minimal Supersymmetric Standard 
Model,” Phys. Rept. 496 (2010), 1-77 doi:10.1016/j.physrep.2010.07.001 [hep- 
phl]. 

[72] K. Cheung, T. J. Hou, J. S. Lee and E. Senaha, “The Higgs Boson Sector of the Next-to-MSSM with 
CP Violation,” Phys. Rev. D 82 (2010), 075007 doi:10.1103/PhysRevD.82.075007 
hep-ph]. 

[73] S. Weinberg, “Gauge Theory of CP Violation,” Phys. Rev. Lett. 37 (1976), 657 

doi:10.1103/PhysRevLett.37.657. 

[74] G. C. Branco, A. J. Buras and J. M. Gerard, “CP Violation in Models With Two and Three Scalar 

Doublets,” Nucl. Phys. B 259 (1985), 306 doi:10.1016 /0550-3213(85)90638-8. 


[75] I. P. Ivanov, “Building and testing models with extended Higgs sectors,” Prog. Part. Nucl. Phys. 

95 (2017), 160-208 doi:10.1016 /j-ppnp.2017.03.001 farXiv:1702.03776| [hep-ph]]. 

[76] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, “The Higgs Hunter’s Guide,” Front. Phys. 

80 (2000), 1-404 SCIPP-89/13. 

[77] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, “Errata for the Higgs hunter’s guide,” 

arXiv:hep-ph/9302272) [hep-ph]]. 

[78] A. Djouadi, “The Anatomy of electro-weak symmetry breaking. I: The Higgs boson in the standard 
model,” Phys. Rept. 457 (2008), 1-216 doi:10.1016/j.physrep.2007.10.004 |arXiv:hep-ph/0503172 
hep-ph]. 

[79] A. Djouadi, “The Anatomy of electro-weak symmetry breaking. II. The Higgs bosons in the min- 
imal supersymmetric model,” Phys. Rept. 459 (2008), 1-241 doi:10.1016/j.physrep.2007.10.005 
arXiv:hep-ph/0503173 [hep-ph]]. 

[80] E. Accomando, A. G. Akeroyd, E. Akhmetzyanova, J. Albert, A. Alves, N. Amapane, M. Aoki, 
G. Azuelos, S. Baffioni and A. Ballestrero, et al. “Workshop on CP Studies and Non-Standard 
Higgs Physics,” doi:10.5170/CERN-2006-009 |arXiv:hep-ph/0608079 [hep-ph]. 

[81] S. Dittmaier et al. [LHC Higgs Cross Section Working Group], “Handbook of LHC Higgs Cross 
Sections: 1. Inclusive Observables,” doi:10.5170/CERN-2011-002 |arXiv:1101.0593) [hep-ph]]. 

[82] S. Dittmaier, C. Mariotti, G. Passarino, R. Tanaka, S. Alekhin, J. Alwall, E. A. Bagnaschi, A. Banfi, 
J. Blumlein and S. Bolognesi, et al. “Handbook of LHC Higgs Cross Sections: 2. Differential 
Distributions,” doi:10.5170/CERN-2012-002 |arXiv:1201.3084 [hep-ph]]. 

[83] S. Heinemeyer et al. [LHC Higgs Cross Section Working Group], “Handbook of LHC Higgs Cross 

Sections: 3. Higgs Properties,” doi:10.5170/CERN-2013-004 |arXiv:1307.1347| [hep-ph]]. 


88 


[84] 


[85] 


[86] 


87 


88 


89 


[90] 


[91] 


97 


98 


99 


D. de Florian et al. [LHC Higgs Cross Section Working Group], “Handbook of LHC Higgs Cross 
Sections: 4. Deciphering the Nature of the Higgs Sector,” doi:10.2172/1345634, 10.23731/CYRM- 
2017-002 arXiv:1610.07922 [hep-ph]. 

S. Dawson, A. Gritsan, H. Logan, J. Qian, C. Tully, R. Van Kooten, A. Ajaib, A. Anastassov, 
I. Anderson and D. Asner, et al. “Working Group Report: Higgs Boson,” |arXiv:1310.8361) [hep- 
ex]]. 

M. Spira, “QCD effects in Higgs physics,’ Fortsch. Phys. 46 (1998), 203-284 
doi:10.1002/(SICT)1521-3978(199804)46:3<203::AID-PROP203>3.0.CO;2-4 
ph/9705337 {hep-ph]]. 

M. Spira, “Higgs Boson Production and Decay at Hadron Colliders,” Prog. Part. Nucl. Phys. 95 
(2017), 98-159 doi:10.1016/j.ppnp.2017.04.001 |arXiv:1612.07651 [hep-ph]]. 

S. Dawson, C. Englert and T. Plehn, “Higgs Physics: It ain’t over till it’s over,” Phys. Rept. 816 
(2019), 1-85 doi:10.1016/j.physrep.2019.05.001 |arXiv:1808.01324 [hep-ph]]. 

A. Djouadi, J. Kalinowski and M. Spira, “HDECAY: A Program for Higgs boson decays in the 
standard model and its supersymmetric extension,” Comput. Phys. Commun. 108 (1998), 56-74 
doi:10.1016/S0010-4655(97)00123-9 |arXiv:hep-ph/9704448 [hep-ph]. 

A. Djouadi, J. Kalinowski, M. Muehlleitner and M. Spira, “HDECAY: Twenty}, years after,” 
Comput. Phys. Commun. 238 (2019), 214-231 doi:10.1016/j.cpc.2018.12.010 jarXiv:1801.09506 
fhep-ph]}. 

A. Bredenstein, A. Denner, S. Dittmaier and M. M. Weber, “Precise predictions for 
the Higgs-boson decay H — WW/ZZ — 4 leptons,” Phys. Rev. D 74 (2006) 013004 
doi:10.1103/PhysRevD.74.013004 |hep-ph/0604011). 

J. S. Lee, A. Pilaftsis, M. S. Carena, S. Y. Choi, M. Drees, J. R. Ellis and C. E. M. Wagner, “CP- 
superH: A Computational tool for Higgs phenomenology in the minimal supersymmetric standard 
model with explicit CP violation,” Comput. Phys. Commun. 156 (2004) 283 |hep-ph/0307377). 
J. S. Lee, M. Carena, J. Ellis, A. Pilaftsis and C. E. M. Wagner, “CPsuperH2.0: an Improved 


Computational Tool for Higgs Phenomenology in the MSSM with Explicit CP Violation,” Comput. 
Phys. Commun. 180 (2009) 312 |arXiv:0712.2360| [hep-ph]]. 

J. S. Lee, M. Carena, J. Ellis, A. Pilaftsis and C. E. M. Wagner, “CPsuperH2.3: an Updated Tool 
for Phenomenology in the MSSM with Explicit CP Violation,” jarXiv:1208.2212) [hep-ph]. 


M. Frank, T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, “The Higgs Boson 
Masses and Mixings of the Complex MSSM in the Feynman-Diagrammatic Approach,” JHEP 02 
(2007), 047 doi:10.1088 /1126-6708 /2007/02/047 \arXiv:hep-ph/0611326) [hep-ph]J]. 

T. Hahn, S. Heinemeyer, F. von der Pahlen, H. Rzehak and C. Schappacher, “Fully Auto- 
mated Calculations in the complex MSSM,” PoS LL2014 (2014), 080 doi:10.22323/1.211.0080 
arXiv: 1407.0231 [hep-ph]]. 

T. Ibrahim and P. Nath, “CP Violation From Standard Model to Strings,” Rev. Mod. Phys. 80 
(2008), 577-631 doi:10.1103/RevModPhys.80.577 |arXiv:0705.2008 [hep-ph]. 

R. E. Shrock and M. Suzuki, “Invisible Decays of Higgs Bosons,” Phys. Lett. B 110 (1982), 250 
doi:10.1016/0370-2693(82)91247-3. 

L. F. Li, Y. Liu and L. Wolfenstein, “HIDDEN HIGGS PARTICLES,” Phys. Lett. B 159 (1985), 
45-48 doi:10.1016/0370-2693(85)90117-0. 


100] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, T. Liu, Z. Liu, D. McKeen, J. Shelton and 


M. Strassler, et al. “Exotic decays of the 125 GeV Higgs boson,” Phys. Rev. D 90 (2014) no.7, 


075004 doi:10.1103/PhysRevD.90.075004 |arXiv:1312.4992 [hep-ph]. 


89 


101) N. D. Christensen and C. Duhr, “FeynRules - Feynman rules made easy,” Comput. Phys. Com- 
mun. 180 (2009), 1614-1641 doi:10.1016/j.cpc.2009.02.018 |arXiv:0806.4194) [hep-ph]]. 


102| C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer and T. Reiter, Comput. Phys. 

Commun. 183 (2012), 1201-1214 doi:10.1016/j.cpc.2012.01.022 |arXiv:1108.2040 [hep-ph]]. 

103] A. Alloul, N. D. Christensen, C. Degrande, C. Duhr and B. Fuks, “FeynRules 2.0 - A com- 
plete toolbox for tree-level phenomenology,” Comput. Phys. Commun. 185 (2014), 2250-2300 
doi:10.1016/j.cpc.2014.04.012 |arXiv:1310.1921 [hep-ph]]. 

104] F. Staub, “SARAH,” |arXiv:0806.0538 [hep-ph]]. 


105] F. Staub, “SARAH 3.2: Dirac Gauginos, UFO output, and more,” Comput. Phys. Commun. 184 


(2013), 1792-1809 doi:10.1016/j.cpc.2013.02.019 |arXiv:1207.0906 [hep-ph]. 


106] F. Staub, “SARAH 4: A tool for (not only SUSY) model builders,” Comput. Phys. Commun. 


185 (2014), 1773-1790 doi:10.1016/j.cpc.2014.02.018 farXiv:1309.7223} [hep-ph]]. 


107] F. Staub, “Exploring new models in all detail with SARAH,” Adv. High Energy Phys. 2015 


(2015), 840780 doi:10.1155/2015/840780 [arXiv:1503.04200) [hep-ph]. 


108] A. V. Semenov, “LanHEP: A Package for automatic generation of Feynman rules in gauge models,” 
arXiv:hep-ph/9608488) [hep-ph]. 

109) A. V. Semenov, “Automatic generation of Feynman rules from the Lagrangian by means of Lan- 
HEP package,” Nucl. Instrum. Meth. A 389 (1997), 293-294 doi:10.1016/50168-9002(97)00096-X. 


110] A. Semenov, “LanHEP: A package for automatic generation of Feynman rules from the La- 
grangian,” Comput. Phys. Commun. 115 (1998), 124-139 doi:10.1016/S0010-4655(98)00143-X. 


111] A. V. Semenov, “LanHEP: A Package for automatic generation of Feynman rules in field theory. 

Version 2.0,” |arXiv:hep-ph/0208011) [hep-ph]. 

112] A. Semenov, “LanHEP: A Package for the automatic generation of Feynman rules in field the- 
ory. Version 3.0,” Comput. Phys. Commun. 180 (2009), 431-454 doi:10.1016/j.cpc.2008.10.012 


arXiv:0805.0555 [hep-ph]. 


[113] A. Semenov, “LanHEP - a package for automatic generation of Feynman rules from the La- 


grangian. Updated version 3.1,” |arXiv:1005.1909 [hep-ph]. 


[114] A. Semenov, “LanHEP — A package for automatic generation of Feynman rules 
from the Lagrangian. Version 3.2,” Comput. Phys. Commun. 201 (2016), 167-170 


doi:10.1016/j.cpce.2016.01.003 |arXiv:1412.5016 [physics.comp-ph]]. 

[115] D. M. Webber et al. [MuLan], “Measurement of the Positive Muon Lifetime and Determina- 
tion of the Fermi Constant to Part-per-Million Precision,” Phys. Rev. Lett. 106 (2011), 041803 
doi:10.1103/PhysRevLett.106.079901 |arXiv:1010.0991 [hep-ex]]. 

[116] V. Tishchenko et al. [MuLan], “Detailed Report of the MuLan Measurement of the Positive 
Muon Lifetime and Determination of the Fermi Constant,” Phys. Rev. D 87 (2013) no.5, 052003 
doi:10.1103/PhysRevD.87.052003 |arXiv:1211.0960 [hep-ex]]. 

117| J. S. Lee and J. Park, “Vacuum Structure in the Standard Model with a Complex Scalar Singlet,” 

New Phys. Sae Mulli 70 (2020) no.7, 571-574 doi:10.3938/NPSM.70.571. 

118] K. Cheung, J. S. Lee and P. Y. Tseng, “Higgcision in the Two-Higgs Doublet Models,” JHEP 01 

(2014), 085 doi:10.1007/JHEP01(2014)085 |arXiv:1310.3937 [hep-ph]]. 


119| S. Davidson and H. E. Haber, “Basis-independent methods for the two-Higgs-doublet model,” 
Phys. Rev. D 72 (2005) 035004 [Erratum-ibid. D 72 (2005) 099902] |hep-ph/0504050). 


120| S. L. Glashow and S. Weinberg, “Natural Conservation Laws for Neutral Currents,” Phys. Rev. 
D 15 (1977) 1958. 


90 


[121] S. Choi and J. S. Lee, “Decays of the MSSM Higgs bosons with explicit CP violation,” Phys. Rev. 


D 61 (1999), 015003 doi:10.1103/PhysRevD.61.015003 {arXiv:hep-ph/9907496 [hep-ph]J]. 


[122] S. Choi, K. Hagiwara and J. S. Lee, “Higgs boson decays in the minimal supersymmetric 
standard model with radiative Higgs sector CP violation,’ Phys. Rev. D 64 (2001), 032004 


doi:10.1103/PhysRevD.64.032004 |arXiv:hep-ph/0103294| [hep-ph]]. 

[123] S. Choi, M. Drees, J. S. Lee and J. Song, “Supersymmetric Higgs boson decays in the MSSM 
with explicit CP violation,” Eur. Phys. J. C 25 (2002), 307-313 doi:10.1007/s10052-002-0997-8 
arXiv:hep-ph/0204200) [hep-ph]. 

[124] M. Carena, J. R. Ellis, S. Mrenna, A. Pilaftsis and C. Wagner, “Collider Probes of the MSSM 
Higgs Sector with Explicit CP Violation,” Nucl. Phys. B 659 (2003), 145-178 doi:10.1016/S0550- 
3213(03)00192-5 |arXiv:hep-ph/0211467) [hep-ph]. 

[125] A. Pilaftsis, “CP odd tadpole renormalization of Higgs scalar - pseudoscalar mixing,” Phys. Rev. 


D 58 (1998), 096010 doi:10.1103/PhysRevD.58.096010 [hep-ph]. 

[126] A. Pilaftsis, “Higgs scalar - pseudoscalar mixing in the minimal supersymmetric standard model,” 
Phys. Lett. B 435 (1998), 88-100 doi:10.1016 /S0370-2693(98)00771-0 [hep- 
phi). 

127| A. Pilaftsis and C. E. M. Wagner, “Higgs bosons in the minimal supersymmetric standard model 

with explicit CP violation,” Nucl. Phys. B 553 (1999) 3 |hep-ph/9902371). 


128] D. A. Demir, “Effects of the supersymmetric phases on the neutral Higgs sector,” Phys. Rev. D 


60 (1999), 055006 doi:10.1103/PhysRevD.60.055006 {arXiv:hep-ph/9901389 [hep-ph]]. 


129| S. Y. Choi, M. Drees and J. S. Lee, “Loop corrections to the neutral Higgs boson sector of 
the MSSM with explicit CP violation,” Phys. Lett. B 481 (2000), 57-66 doi:10.1016/S0370- 
2693(00)00421-4 [hep-ph]. 

[130] M. S. Carena, J. R. Ellis, A. Pilaftsis and C. E. M. Wagner, “Renormalization group improved 

effective potential for the MSSM Higgs sector with explicit CP violation,” Nucl. Phys. B 586 (2000) 

92 (hep-ph /0003180). 


131] M. S. Carena, J. R. Ellis, A. Pilaftsis and C. E. M. Wagner, “Higgs boson pole masses in the 


MSSM with explicit CP violation,” Nucl. Phys. B 625 (2002) 345 |hep-ph/0111245). 


132] R. Hempfling, “Yukawa coupling unification with supersymmetric threshold corrections,” Phys. 
Rev. D 49 (1994), 6168-6172 doi:10.1103/PhysRevD.49.6168. 


133] L. J. Hall, R. Rattazzi and U. Sarid, “The Top quark mass in supersymmetric SO(10) unification,” 
Phys. Rev. D 50 (1994), 7048-7065 doi:10.1103/PhysRevD.50.7048 [hep- 
phl]. 

[134] M. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, “Electroweak symmetry break- 
ing and bottom - top Yukawa unification,” Nucl. Phys. B 426 (1994), 269-300 doi:10.1016/0550- 
3213(94)90313-1 hep-ph) 

[135] D. M. Pierce, J. A. Bagger, K. T. Matchev and R. j. Zhang, “Precision corrections in the 
minimal supersymmetric standard model,” Nucl. Phys. B 491 (1997), 3-67 doi:10.1016/S0550- 
3213(96)00683-9 hep-ph] 

[136] F. Borzumati, G. R. Farrar, N. Polonsky and S. D. Thomas, “Soft Yukawa couplings in supersym- 
metric theories,” Nucl. Phys. B 555 (1999), 53-115 doi:10.1016/S0550-3213(99)00328-4 
{hep-ph}. 

[137] K. S. Babu and C. F. Kolda, “Signatures of supersymmetry and Yukawa unification in Higgs de- 

cays,” Phys. Lett. B 451 (1999), 77-85 doi:10.1016/S0370-2693(99)00204-X 

hep-ph]. 


91 


[138] A. Pilaftsis, “Higgs mediated electric dipole moments in the MSSM: An application to baryoge- 
nesis and Higgs searches,” Nucl. Phys. B 644 (2002), 263-289 doi:10.1016/50550-3213(02)00826-X 
hep-ph] 

[139] F. Borzumati, J. S. Lee and W. Y. Song, “Threshold corrections to m(b) and the b anti- 
b —> HO(i) production in CP violating SUSY scenarios,” Phys. Lett. B 595 (2004), 347-358 
doi:10.1016/j-physletb.2004.06.040 [hep-ph]]. 

[140] J. S. Lee, Y. Peters, A. Pilaftsis and C. Schwanenberger, “Strangephilic Higgs Bosons in the 
MSSM,” Eur. Phys. J. C 66 (2010), 261-269 doi:10.1140/epjc/s10052-010-1248-z 
hep-ph]. 

141] L. Mihaila, B. Schmidt and M. Steinhauser, “T(H — bb) to order aa,,” Phys. Lett. B 751 (2015) 

442 doi:10.1016/j.physletb.2015.10.078 [hep-ph]]. 

142| E. Braaten and J. P. Leveille, “Higgs Boson Decay and the Running Mass,” Phys. Rev. D 22 

(1980) 715. doi:10.1103/PhysRevD.22.715. 


143| N. Sakai, “Perturbative QCD Corrections to the Hadronic Decay Width of the Higgs Boson,” 
Phys. Rev. D 22 (1980) 2220. doi:10.1103/PhysRevD.22.2220. 


144) T. Inami and T. Kubota, “Renormalization Group Estimate of the Hadronic Decay Width of the 
Higgs Boson,” Nucl. Phys. B 179 (1981) 171. doi:10.1016/0550-3213(81)90253-4. 

145| M. Drees and K. i. Hikasa, “Note On Qcd Corrections To Hadronic Higgs Decay,” Phys. Lett. B 
240 (1990) 455 Erratum: [Phys. Lett. B 262 (1991) 497]. doi:10.1016/0370-2693(90)91130-4. 

146] M. Drees and K. i. Hikasa, “Heavy Quark Thresholds in Higgs Physics,” Phys. Rev. D 41 (1990) 
1547. doi:10.1103/PhysRevD.41.1547. 


147| S. G. Gorishnii, A. L. Kataev and S. A. Larin, “The Width of Higgs Boson Decay Into Hadrons: 
Three Loop Corrections of Strong Interactions,” Sov. J. Nucl. Phys. 40 (1984) 329 [Yad. Fiz. 40 
(1984) 517]. 

[148] S. G. Gorishnii, A. L. Kataev, S. A. Larin and L. R. Surguladze, “Corrected Three Loop QCD 

Correction to the Correlator of the Quark Scalar Currents and y (Tot) (H° — Hadrons),” Mod. 

Phys. Lett. A 5 (1990) 2703. doi:10.1142/S0217732390003152. 


[149] S. G. Gorishnii, A. L. Kataev, S. A. Larin and L. R. Surguladze, “Scheme dependence of the 
next to next-to-leading QCD corrections to Gamma(tot) (HO — hadrons) and the spurious QCD 
infrared fixed point,” Phys. Rev. D 43 (1991) 1633. doi:10.1103/PhysRevD.43. 1633. 


150) A. L. Kataev and V. T. Kim, “The Effects of the QCD corrections to Gamma (H0 — b anti-b),” 
Mod. Phys. Lett. A 9 (1994) 1309. doi:10.1142/S0217732394001131. 


151] L. R. Surguladze, “Quark mass effects in fermionic decays of the Higgs boson in O (alpha-s**2) per- 
turbative QCD,” Phys. Lett. B 341 (1994) 60 doi:10.1016 /0370-2693(94)01253-9 |hep-ph/9405325). 


152] K. G. Chetyrkin, “Correlator of the quark scalar currents and Gamma(tot) (H — hadrons) at O 
(alpha-s**3) in pQCD,” Phys. Lett. B 390 (1997) 309 doi:10.1016/S0370-2693(96)01368-8 


ph /9608318) 

[153] K. Melnikov, “Two loop O(N(f) alpha-s**2) correction to the decay width of the Higgs bo- 
son to two massive fermions,” Phys. Rev. D 53 (1996) 5020 doi:10.1103/PhysRevD.53.5020 
ph 9511310) 

[154] K. G. Chetyrkin and A. Kwiatkowski, “Second order QCD corrections to scalar and pseu- 
doscalar Higgs decays into massive bottom quarks,” Nucl. Phys. B 461 (1996) 3 doi:10.1016/0550- 
3213(95)00616-8 [hep-ph /9505358}. 

[155] S. A. Larin, T. van Ritbergen and J. A. M. Vermaseren, “The Large top quark mass expan- 
sion for Higgs boson decays into bottom quarks and into gluons,” Phys. Lett. B 362 (1995) 134 


doi:10.1016/0370-2693(95)01192-S |hep-ph/9506465). 


92 


156| J. Fleischer and F. Jegerlehner, “Radiative Corrections to Higgs Decays in the Extended Weinberg- 
Salam Model,” Phys. Rev. D 23 (1981) 2001. doi:10.1103/PhysRevD.23.2001. 


157| D. Y. Bardin, B. M. Vilensky and P. K. Khristova, “Calculation of the Higgs boson decay width 
into fermion pairs,” Sov. J. Nucl. Phys. 53 (1991) 152 [Yad. Fiz. 53 (1991) 240). 


158] A. Dabelstein and W. Hollik, “Electroweak corrections to the fermionic decay width of the stan- 
dard Higgs boson,” Z. Phys. C 53 (1992) 507. doi:10.1007/BF01625912. 


159] B. A. Kniehl, “Radiative corrections for H —> f anti-f (y) in the standard model,” Nucl. Phys. B 
376 (1992) 3. doi:10.1016/0550-3213(92)90065-J. 


160] A. Djouadi, D. Haidt, B. A. Kniehl, P. M. Zerwas and B. Mele, “Higgs in the standard model,” 
In *Munich/Annecy/Hamburg 1991, Proceedings, e+ e- collisions at 500-GeV, pt. A* 11-30. 


161] B. A. Kniehl and M. Spira, “Two loop O (alpha-s G(F) m(t)? correction to the H — bb decay 
rate,” Nucl. Phys. B 432 (1994) 39 doi:10.1016/0550-3213(94)90592-4 |hep-ph/9410319). 


162] A. Kwiatkowski and M. Steinhauser, “Corrections of order O(Grasm?) to the Higgs decay 
rate T(H — bb)” Phys. Lett. B 338 (1994) 66 Erratum: [Phys. Lett. B 342 (1995) 455] 


doi:10.1016/0370-2693(94)01527-J, 10.1016/0370-2693(94)91345-5 |hep-ph/9405308). 
[163] K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, “Virtual top quark effects on the 
H — b anti-b decay at next-to-leading order in QCD,” Phys. Rev. Lett. 78 (1997) 594 


doi:10.1103/PhysRevLett.78.594 [hep-ph/9610456). 

[164] A. L. Kataev, “The Order O (alpha alpha-s) and O (alpha**2) corrections to the decay width of 
the neutral Higgs boson to the anti-b b pair,” JETP Lett. 66 (1997), 327-330 doi:10.1134/1.567516 
hep-ph]. 

[165] A. Grau, G. Panchieri and R. J. N. Phillips, “Contributions of off-shell top quarks to decay 
processes,” Phys. Lett. B 251 (1990) 293. doi:10.1016/0370-2693(90)90939-4. 


[166] S. Moretti and W. J. Stirling, “Contributions of below threshold decays to MSSM Higgs branch- 
ing ratios,” Phys. Lett. B 347 (1995), 291-299 |erratum: Phys. Lett. B 366 (1996), 451] 
doi:10.1016/0370-2693(95)00088-3 [hep-ph]]. 

|167] J. Chang, K. Cheung, J. S. Lee, C. T. Lu and J. Park, “Measuring properties of a Heavy Higgs 
boson in the H > tt — bW*bW~ decay,” JHEP 05 (2018), 162 doi:10.1007/JHEP05(2018)162 


fr Xivi1803.05762, [hep-ph] 
[168] S. Dawson, A. Djouadi and M. Spira, “QCD corrections to SUSY Higgs production: The Role 
of squark loops,” Phys. Rev. Lett. 77 (1996), 16-19 doi:10.1103/PhysRevLett.77.16 
ph/9603423| [hep-ph]. 
[169] H. Eberl, K. Hidaka, S. Kraml, W. Majerotto and Y. Yamada, “Improved SUSY QCD cor- 
rections to Higgs boson decays into quarks and squarks,” Phys. Rev. D 62 (2000), 055006 


doi:10.1103/PhysRevD.62.055006 |arXiv:hep-ph/9912463) [hep-ph]]. 
[170] C. Weber, K. Kovarik, H. Eberl and W. Majerotto, “Complete one-loop corrections to decays of 
charged and CP-even neutral Higgs bosons into sfermions,” Nucl. Phys. B 776 (2007), 138-169 


doi:10.1016/j-nuclphysb.2007.03.048 [hep-ph]. 

[171] E. Accomando, G. Chachamis, F. Fugel, M. Spira and M. Walser, “Neutral Higgs Boson Decays 
to Squark Pairs reanalyzed,” Phys. Rev. D 85 (2012), 015004 doi:10.1103/PhysRevD.85.015004 
hep-ph]. 


[172] S. Heinemeyer and C. Schappacher, “Heavy Higgs Decays into Sfermions in the Complex MSSM: A 
Full One-Loop Analysis,” Eur. Phys. J. C 75 (2015) no.5, 198 doi:10.1140/epjc/s10052-015-3383-z 
hep-ph]. 

[173] R. N. Cahn, “The Higgs Boson,” Rept. Prog. Phys. 52 (1989) 389. doi:10.1088/0034- 
4885 /52/4/001. 


93 


[174] A. Djouadi, J. Kalinowski and P. M. Zerwas, “Two and three-body decay modes of SUSY Higgs 
particles,” Z. Phys. C 70 (1996), 435-448 doi:10.1007/s002880050121 [hep- 
phi]. 

[175] W. Y. Keung and W. J. Marciano, “HIGGS SCALAR DECAYS: H —> W+- X,” Phys. Rev. D 
30 (1984), 248 doi:10.1103/PhysRevD.30.248. 


[176] A. Bredenstein, A. Denner, S. Dittmaier and M. M. Weber, “Precision calculations for the 
Higgs decays H —> ZZ/WW —> 4 leptons,” Nucl. Phys. B Proc. Suppl. 160 (2006), 131-135 


doi:10.1016/j.nuclphysbps.2006.09.104 |arXiv:hep-ph/0607060| [hep-ph]. 

177] A. Bredenstein, A. Denner, S. Dittmaier and M. M. Weber, “Radiative corrections to the semilep- 
tonic and hadronic Higgs-boson decays H —> W W / Z Z —> 4 fermions,” JHEP 02 (2007), 080 
doi:10.1088/1126-6708 /2007 /02/080 |arXiv:hep-ph/0611234 [hep-ph]]. 

178] L. Altenkamp, S. Dittmaier and H. Rzehak, “Renormalization schemes for the Two-Higgs- 

Doublet Model and applications to h — WW/ZZ — 4 fermions,” JHEP 09 (2017), 134 


doi:10.1007/JHEP09(2017)134 |arXiv:1704.02645 [hep-ph]]. 

179] L. Altenkamp, S. Dittmaier and H. Rzehak, “Precision calculations for h => WW/ZZ > 

4 fermions in the Two-Higgs-Doublet Model with Prophecy4f,” JHEP 03 (2018), 110 

doi:10.1007/JHEP03(2018)110 |arXiv:1710.07598 [hep-ph]]. 

180] L. Altenkamp, M. Boggia and S. Dittmaier, “Precision calculations for h > WW/ZZ > 4 
fermions in a Singlet Extension of the Standard Model with Prophecy4f,” JHEP 04 (2018), 062 
doi:10.1007/JHEP04(2018)062 farXiv:1801.07291| [hep-ph]]. 

[181] A. Denner, S. Dittmaier and J. N. Lang, “Renormalization of mixing angles,” JHEP 11 (2018), 
104 doi:10.1007/JHEP11(2018)104 |arXiv:1808.03466 [hep-ph]]. 

[182] A. Denner, S. Heinemeyer, I. Puljak, D. Rebuzzi and M. Spira, “Standard Model Higgs-Boson 
Branching Ratios with Uncertainties,” Eur. Phys. J. C 71 (2011) 1753 doi:10.1140/epjc/s10052- 
011-1753-8 |arXiv:1107.5909 [hep-ph]]. 

183] M. Spira, A. Djouadi, D. Graudenz and P. M. Zerwas, “Higgs boson production at the LHC,” 

Nucl. Phys. B 453 (1995) 17 doi:10.1016/0550-3213(95)00379-7 |hep-ph/9504378}. 

184] A. Djouadi, M. Spira and P. M. Zerwas, “Production of Higgs bosons in proton colliders: QCD 

corrections,” Phys. Lett. B 264 (1991) 440. doi:10.1016/0370-2693(91)90375-Z. 

185] K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, “Hadronic Higgs decay to order alpha-s**4,” 

Phys. Rev. Lett. 79 (1997) 353 doi:10.1103/PhysRevLett.79.353 |hep-ph/9705240). 

186] P. A. Baikov and K. G. Chetyrkin, “Top Quark Mediated Higgs Boson Decay into Hadrons to Or- 

der a®,” Phys. Rev. Lett. 97 (2006) 061803 doi:10.1103/PhysRevLett.97.061803 |hep-ph/0604194). 

187| K. Melnikov, M. Spira and O. I. Yakovlev, “Threshold effects in two photon decays of Higgs 
particles,” Z. Phys. C 64 (1994), 401-406 doi:10.1007/BF01560100 |arXiv:hep-ph/9405301, [hep- 
phi). 

188] V. S. Fadin and V. A. Khoze, “Threshold Behavior of Heavy Top Production in e+ e- Collisions,” 

JETP Lett. 46 (1987), 525-529 LENINGRAD-87-1333. 

189] V. S. Fadin and V. A. Khoze, “Production of a pair of heavy quarks in e+ e- annihilation in the 

threshold region,” Sov. J. Nucl. Phys. 48 (1988), 309-313. 

190| M. J. Strassler and M. E. Peskin, “The Heavy top quark threshold: QCD and the Higgs,” Phys. 

Rev. D 43 (1991), 1500-1514 doi:10.1103/PhysRevD.43.1500. 

191] K. G. Chetyrkin, B. A. Kniehl, M. Steinhauser and W. A. Bardeen, “Effective QCD interac- 

tions of CP odd Higgs bosons at three loops,” Nucl. Phys. B 535 (1998) 3 doi:10.1016/S0550- 

3213(98)00594-X |hep-ph/9807241). 


94 


[192] A. Djouadi and P. Gambino, “Leading electroweak correction to Higgs boson production at proton 


colliders,” Phys. Rev. Lett. 73 (1994) 2528 doi:10.1103/PhysRevLett.73.2528 |hep-ph/9406432). 
[193] K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, “Three loop O (alpha-s**2 G(F) M(t)**2) cor- 
rections to hadronic Higgs decays,” Nucl. Phys. B 490 (1997) 19 doi:10.1016/S0550-3213(97)00051- 
5 |hep-ph/9701277]. 
[194] U. Aglietti, R. Bonciani, G. Degrassi and A. Vicini, “Two loop light fermion contribution to 
Higgs production and decays,” Phys. Lett. B 595 (2004) 432 doi:10.1016/j.physletb.2004.06.063 


hep-ph /0404071). 


195] U. Aglietti, R. Bonciani, G. Degrassi and A. Vicini, “Two-loop electroweak corrections to Higgs 


production in proton-proton collisions,” hep-ph /0610033 
196| G. Degrassi and F. Maltoni, “Two-loop electroweak corrections to Higgs production at hadron 
colliders,” Phys. Lett. B 600 (2004) 255 doi:10.1016/j.physletb.2004.09.008 |hep-ph/0407249}. 


197] S. Actis, G. Passarino, C. Sturm and S. Uccirati, “NLO Electroweak Corrections to Higgs Boson 
Production at Hadron Colliders,” Phys. Lett. B 670 (2008) 12 doi:10.1016/j.physletb.2008.10.018 


ar Xiv:0809.1301, [hep-ph]. 


[198] S. Actis, G. Passarino, C. Sturm and S. Uccirati, “NNLO Computational Techniques: The Cases H 
— gamma gamma and H > g g,” Nucl. Phys. B 811 (2009) 182 doi:10.1016/j.nuclphysb.2008.11.024 


ar Xiv:0809.3667 [hep-ph]. 


[199] S. Actis, G. Passarino, C. Sturm and S. Uccirati, “Two-Loop Threshold Singularities, Unstable 
Particles and Complex Masses,” Phys. Lett. B 669 (2008), 62-68 doi:10.1016/j.physletb.2008.09.028 


hep-ph] 
[200] G. Passarino, C. Sturm and S. Uccirati, “Complete Two-Loop Corrections to H —> gamma 
gamma,” Phys. Lett. B 655 (2007), 298-306 doi:10.1016/j.physletb.2007.09.002 
[hep-ph]. 
[201] G. Degrassi and F. Maltoni, “Two-loop electroweak corrections to the Higgs-boson decay H 
— gamma gamma,” Nucl. Phys. B 724 (2005) 183 doi:10.1016/j.nuclphysb.2005.06.027 


ph/0504137]. 


[202] J. Brod, F. Fugel and B. A. Kniehl, “Dominant two-loop electroweak corrections to the hadropro- 
duction of a pseudoscalar Higgs boson and its photonic decay,” Phys. Rev. D 78 (2008) 011303 
doi:10.1103/PhysRevD.78.011303 {arXiv:0802.0171 [hep-ph]]. 

[203] J. Brod, F. Fugel and B. A. Kniehl, “T'wo-Loop Electroweak Corrections to the AO gamma 
gamma and AO g g Couplings of the CP-Odd Higgs Boson,” Nucl. Phys. B 807 (2009) 188 
doi:10.1016/j.nuclphysb.2008.08.021 |arXiv:0807. 1008) [hep-ph]]. 

[204] M. Muhlleitner and M. Spira, “Higgs Boson Production via Gluon Fusion: Squark Loops at NLO 
QCD,” Nucl. Phys. B 790 (2008) 1 doi:10.1016/j-nuclphysb.2007.08.011 |hep-ph/0612254). 


[205] R. Bonciani, G. Degrassi and A. Vicini, “Scalar particle contribution to Higgs production via gluon 


fusion at NLO,” JHEP 0711 (2007) 095 doi:10.1088/1126-6708/2007/11/095 |arXiv:0709.4227 

hep-ph]. 

[206] H. Q. Zheng and D. D. Wu, “First order QCD corrections to the decay of the Higgs boson into 
two photons,” Phys. Rev. D 42 (1990) 3760. doi:10.1103/PhysRevD.42.3760. 

[207] A. Djouadi, M. Spira, J. J. van der Bij and P. M. Zerwas, “QCD corrections to gamma 
gamma decays of Higgs particles in the intermediate mass range,” Phys. Lett. B 257 (1991) 187. 
doi:10.1016/0370-2693(91)90879-U. 

[208] S. Dawson and R. P. Kauffman, “QCD corrections to H > gamma gamma,” Phys. Rev. D 47 
(1993) 1264. doi:10.1103/PhysRevD.47.1264. 


— 


95 


209| K. Melnikov and O. I. Yakovlev, “Higgs —> two photon decay: QCD radiative correction,” Phys. 


Lett. B 312 (1993) 179 doi:10.1016/0370-2693(93)90507-E |hep-ph/9302281). 


210] A. Djouadi, M. Spira and P. M. Zerwas, “Two photon decay widths of Higgs particles,” Phys. 
Lett. B 311 (1993) 255 doi:10.1016/0370-2693(93)90564-X |hep-ph/9305335]. 


211| M. Inoue, R. Najima, T. Oka and J. Saito, “QCD corrections to two photon decay of the Higgs 
boson and its reverse process,” Mod. Phys. Lett. A 9 (1994) 1189. doi:10.1142/S0217732394001003. 


212| J. Fleischer, O. V. Tarasov and V. O. Tarasov, “Analytical result for the two loop QCD correction 
to the decay H —> 2 gamma,” Phys. Lett. B 584 (2004) 294 doi:10.1016/j.physletb.2004.01.063 


hep-ph /0401090]. 


[213] R. Harlander and P. Kant, “Higgs production and decay: Analytic results at next-to-leading order 


QCD,” JHEP 0512 (2005) 015 doi:10.1088/1126-6708/2005/12/015 |hep-ph/0509189). 


[214] C. Anastasiou, S. Beerli, S. Bucherer, A. Daleo and Z. Kunszt, “T'wo-loop amplitudes and master 
integrals for the production of a Higgs boson via a massive quark and a scalar-quark loop,” JHEP 


0701 (2007) 082 doi:10.1088/1126-6708/2007/01/082 |hep-ph/0611236). 


215| U. Aglietti, R. Bonciani, G. Degrassi and A. Vicini, “Analytic Results for Virtual QCD Corrections 
to Higgs Production and Decay,” JHEP 0701 (2007) 021 doi:10.1088/1126-6708/2007/01/021 
ph/0611266). 

216] A. Djouadi, V. Driesen, W. Hollik and J. I. Ilana, “The Coupling of the lightest SUSY Higgs boson 
to two photons in the decoupling regime,” Eur. Phys. J. C 1 (1998) 149 doi:10.1007/BF01245805 
hep-ph /9612362|. 

217] J. R. Ellis, M. K. Gaillard and D. V. Nanopoulos, “A Phenomenological Profile of the Higgs 
Boson,” Nucl. Phys. B 106 (1976), 292 doi:10.1016/0550-3213(76)90382-5. 


218] M. A. Shifman, A. I. Vainshtein, M. B. Voloshin and V. I. Zakharov, “Low-Energy Theorems for 
Higgs Boson Couplings to Photons,” Sov. J. Nucl. Phys. 30 (1979), 711-716 ITEP-42-1979. 

219| B. A. Kniehl and M. Spira, “Low-energy theorems in Higgs physics,” Z. Phys. C 69 (1995), 77-88 
doi:10.1007/s002880050007 |arXiv:hep-ph/9505225) [hep-ph]. 


220| O. V. Tarasov, “Anomalous dimensions of quark masses in the three-loop approximation,” Phys. 


Part. Nucl. Lett. 17 (2020) no.2, 109-115 doi:10.1134/81547477120020223 [hep- 
phl]. 

221] S. L. Adler and W. A. Bardeen, “Absence of higher order corrections in the anomalous axial 
vector divergence equation,” Phys. Rev. 182 (1969), 1517-1536 doi:10.1103/PhysRev.182.1517. 
222| S. L. Adler, “Axial vector vertex in spinor electrodynamics,” Phys. Rev. 177 (1969), 2426-2438 
doi:10.1103/PhysRev.177.2426. 

223] J. S. Bell and R. Jackiw, “A PCAC puzzle: 7° — yy in the ø model,” Nuovo Cim. A 60 (1969), 
47-61 doi:10.1007/BF02823296. 


224) A. Djouadi, P. Gambino and B. A. Kniehl, “Two loop electroweak heavy fermion corrections to 

Higgs boson production and decay,” Nucl. Phys. B 523 (1998) 17 doi:10.1016/S0550-3213(98)00147- 

3 |hep-ph/9712330]. 

225| J. F. Gunion and H. E. Haber, “The CP conserving two Higgs doublet model: The Approach to the 
decoupling limit,” Phys. Rev. D 67 (2003), 075019 doi:10.1103/PhysRevD.67.075019 
ph/0207010) [hep-ph]. 

[226] H. E. Haber, “Decoupling and the radiatively corrected MSSM Higgs sector,” Nucl. Phys. B Proc. 

Suppl. 116 (2003), 291-295 doi:10.1016/S0920-5632(03)80186-X |arXiv:hep-ph/0212010) [hep-ph]. 


[227] A. Djouadi, V. Driesen, W. Hollik and A. Kraft, “The Higgs photon - Z boson coupling revisited,” 
Eur. Phys. J. C 1 (1998) 163 |hep-ph/9701342). 


96 


[228] R. N. Cahn, M. S. Chanowitz and N. Fleishon, “Higgs Particle Production by Z => H Gamma,” 
Phys. Lett. 82B (1979) 113. doi:10.1016 /0370-2693(79)90438-6. 

[229] L. Bergstrom and G. Hulth, “Induced Higgs Couplings to Neutral Bosons in e*e~ Collisions,” 
Nucl. Phys. B 259 (1985) 137 Erratum: [Nucl. Phys. B 276 (1986) 744]. doi:10.1016/0550- 
3213(86)90074-X, 10.1016/0550-3213(85)90302-5. 

230| G. Gamberini, G. F. Giudice and G. Ridolfi, “Supersymmetric Higgs Boson Production in Z 

Decays,” Nucl. Phys. B 292 (1987) 237 doi:10.1016/0550-3213(87)90644-4. 

231] M. Spira, A. Djouadi and P. M. Zerwas, “QCD corrections to the H Z gamma coupling,” Phys. 

Lett. B 276 (1992) 350. doi:10.1016/0370-2693(92)90331-W. 


232| T. Gehrmann, S. Guns and D. Kara, “The rare decay H —> Zy in perturbative QCD,” JHEP 

1509 (2015) 038 doi:10.1007/JHEP09(2015)038 |arXiv:1505.00561) [hep-ph]]. 

233] R. Bonciani, V. Del Duca, H. Frellesvig, J. M. Henn, F. Moriello and V. A. Smirnov, “Next- 
to-leading order QCD corrections to the decay width H —> Zy,” JHEP 1508 (2015) 108 
doi:10.1007/JHEP08(2015)108 |arXiv:1505.00567 [hep-ph]]. 

[234] A. Abbasabadi, D. Bowser-Chao, D. A. Dicus and W. W. Repko, “Radiative Higgs bo- 
son decays H —> fermion anti-fermion gamma,” Phys. Rev. D 55 (1997), 5647-5656 
doi:10.1103/PhysRevD.55.5647 |arXiv:hep-ph/9611209 [hep-ph]]. 

[235] A. Abbasabadi and W. W. Repko, “A Note on the rare decay of a Higgs boson into photons and a Z 
boson,” Phys. Rev. D 71 (2005), 017304 doi:10.1103/PhysRevD.71.017304 |arXiv:hep-ph/0411152 
hep-ph]. 

[236] A. Abbasabadi and W. W. Repko, “Higgs boson decay into Z bosons and a photon,” JHEP 08 
(2006), 048 doi:10.1088/1126-6708/2006/08/048 |arXiv:hep-ph/0602087 [hep-ph]]. 

[237] L. B. Chen, C. F. Qiao and R. L. Zhu, “Reconstructing the 125 GeV SM Higgs Boson 
Through 44y,” Phys. Lett. B 726 (2013), 306-311 erratum: Phys. Lett. B 808 (2020), 135629] 
doi:10.1016/j.physletb.2013.08.050 |arXiv:1211.6058 [hep-ph]]. 

238] D. A. Dicus and W. W. Repko, “Calculation of the decay H — eey,” Phys. Rev. D 87 (2013) 

no.7, 077301 doi:10.1103/PhysRevD.87.077301 |arXiv:1302.2159 [hep-ph]]. 

239| Y. Sun, H. R. Chang and D. N. Gao, “Higgs decays to gamma l+ l- in the standard model,” 

JHEP 05 (2013), 061 doi:10.1007/JHEP05(2013)061 arXiv: 1303.2230| [hep-ph]]. 


240| G. Passarino, “Higgs Boson Production and Decay: Dalitz Sector,” Phys. Lett. B 727 (2013), 
424-431 doi:10.1016/j.physletb.2013.10.052 [arXiv:1308.0422| [hep-ph]]. 


241] G. Aad et al. [ATLAS], “Evidence for Higgs boson decays to a low-mass dilepton system and a 

photon in pp collisions at ys = 13 TeV with the ATLAS detector,” |arXiv:2103.10322) [hep-ex]]. 

242] H. E. Haber, M. J. Herrero, H. E. Logan, S. Penaranda, S. Rigolin and D. Temes, “SUSY QCD 
corrections to the MSSM h0 bb vertex in the decoupling limit,” Phys. Rev. D 63 (2001), 055004 
doi:10.1103/PhysRevD.63.055004 |arXiv:hep-ph/0007006) [hep-ph]]. 

[243] P. M. Ferreira, J. F. Gunion, H. E. Haber and R. Santos, “Probing wrong-sign Yukawa 

couplings at the LHC and a future linear collider,” Phys. Rev. D 89 (2014) no.11, 115003 


doi:10.1103/PhysRevD.89.115003 {arXiv:1403.4736 [hep-ph]]. 
[244] S. Y. Choi, J. S. Lee and J. Park, “Alignment of Yukawa couplings in two Higgs doublet models,” 
arXiv:2011.04978 {hep-ph]]. 


[245] J. R. Ellis, J. S. Lee and A. Pilaftsis, “CERN LHC signatures of resonant CP violation in a minimal 
supersymmetric Higgs sector,” Phys. Rev. D 70 (2004), 075010 doi:10.1103/PhysRevD.70.075010 


arXiv:hep-ph/0404167] [hep-ph]. 


97 


[246] S. Y. Choi, J. Kalinowski, Y. Liao and P. M. Zerwas, “H/A Higgs mixing in CP-noninvariant 
supersymmetric theories,” Eur. Phys. J. C 40 (2005), 555-564 doi:10.1140/epjc/s2005-02145-9 
arXiv:hep-ph/0407347| [hep-ph]. 

[247] M. Carena, J. Ellis, J. S. Lee, A. Pilaftsis and C. E. M. Wagner, “CP Violation in Heavy MSSM 
Higgs Scenarios,” JHEP 02 (2016), 123 doi:10.1007/JHEP02(2016)123 farXiv:1512.00437] [hep-ph]. 

[248] A. G. Akeroyd, M. Aoki, A. Arhrib, L. Basso, I. F. Ginzburg, R. Guedes, J. Hernandez-Sanchez, 
K. Huitu, T. Hurth and M. Kadastik, et al. “Prospects for charged Higgs searches at the LHC,” Eur. 
Phys. J. C 77 (2017) no.5, 276 doi:10.1140/epjc/s10052-017-4829-2 |arXiv:1607.01320 [hep-ph]. 

[249] A. Arhrib, E. Christova, H. Eberl and E. Ginina, “CP violation in charged Higgs pro- 
duction and decays in the Complex Two Higgs Doublet Model,” JHEP 04 (2011), 089 
doi:10.1007/JHEP04(2011)089 |arXiv:1011.6560) [hep-ph]]. 

[250] A. Arbey, F. Mahmoudi, O. Stal and T. Stefaniak, “Status of the Charged Higgs Boson in Two 
Higgs Doublet Models,” Eur. Phys. J. C 78 (2018) no.3, 182 doi:10.1140/epjc/s10052-018-5651-1 
arXiv:1706.07414 {hep-ph]]. 

[251] A. Arhrib, R. Benbrik, H. Harouiz, S. Moretti and A. Rouchad, “A Guidebook to Hunting Charged 
Higgs Bosons at the LHC,” doi:10.3389/fphy.2020.00039 farXiv:1810.09106} [hep-ph]. 

[252] A. G. Akeroyd, A. Arhrib and E. M. Naimi, “Yukawa coupling corrections to the decay Ht + W+ 
AO,” Eur. Phys. J. C 12 (2000), 451-460 doi:10.1007/s100529900266 |arXiv:hep-ph/9811431) [hep- 
phj], Eur. Phys. J. C 14 (2000) 371 (erratum). 

[253] A. Akeroyd, A. Arhrib and E. Naimi, “Radiative corrections to the decay Ht + W* A0,” Eur. 
Phys. J. C 20 (2001), 51-62 doi:10.1007/s100520100653 |arXiv:hep-ph/0002288 [hep-ph]]. 

[254] M. Krause, R. Lorenz, M. Muhlleitner, R. Santos and H. Ziesche, “Gauge-independent Renor- 
malization of the 2-Higgs-Doublet Model,” JHEP 09 (2016), 143 doi:10.1007/JHEP09(2016)143 
arXiv:1605.04853 [hep-ph]]. 

[255] M. Krause, M. Miihlleitner and M. Spira, “2HDECAY —A program for the calculation 
of electroweak one-loop corrections to Higgs decays in the T'wo-Higgs-Doublet Model in- 
cluding state-of-the-art QCD corrections,” Comput. Phys. Commun. 246 (2020), 106852 


doi:10.1016/j.cpce.2019.08.003 |arXiv:1810.00768 [hep-ph]]. 

256| K. Cheung, J. S. Lee and P. Y. Tseng, “Higgs Precision (Higgcision) Era begins,” JHEP 1305 

(2013) 134 doi:10.1007/JHEP05(2013) 134 jarXiv:1302.3794 [hep-ph]. 

257| K. Cheung, J. S. Lee and P. Y. Tseng, “Higgs precision analysis updates 2014,” Phys. Rev. D 90 

(2014), 095009 doi:10.1103/PhysRevD.90.095009 |arXiv:1407.8236 [hep-ph]]. 

258] K. Cheung, J. S. Lee and P. Y. Tseng, “New Emerging Results in Higgs Precision Analysis 

Jpdates 2018 after Establishment of Third-Generation Yukawa Couplings,” JHEP 09 (2019), 098 

doi:10.1007/JHEP09(2019)098 |arXiv:1810.02521 [hep-ph]]. 

[259] M. Diithrssen, “Prospects for the measurement of Higgs boson coupling parameters in the mass 

range from 110 - 190 GeV,” ATL-PHYS-2003-030. 

[260] M. Duhrssen, S. Heinemeyer, H. Logan, D. Rainwater, G. Weiglein and D. Zeppenfeld, “Ex- 
tracting Higgs boson couplings from CERN LHC data,” Phys. Rev. D 70 (2004), 113009 
doi:10.1103/PhysRevD.70.113009 |arXiv:hep-ph/0406323) [hep-ph]]. 

[261] R. Lafaye, T. Plehn, M. Rauch, D. Zerwas and M. Duhrssen, “Measuring the Higgs Sector,” JHEP 
08 (2009), 009 doi:10.1088/1126-6708/2009/08/009 |arXiv:0904.3866) [hep-ph]. 

[262] J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott, “Fingerprinting Higgs Suspects at the 
LHC,” JHEP 05 (2012), 097 doi:10.1007/JHEP05(2012)097 |arXiv:1202.3697 [hep-ph]]. 


x 


98 


[263] A. Azatov, R. Contino and J. Galloway, “Model-Independent Bounds on a Light Higgs,” JHEP 
04 (2012), 127 [erratum: JHEP 04 (2013), 140] doi:10.1007/JHEP04(2012)127 
hep-ph]. 

[264] D. Carmi, A. Falkowski, E. Kuflik and T. Volansky, “Interpreting LHC Higgs Results from Natural 
New Physics Perspective,” JHEP 07 (2012), 136 doi:10.1007/JHEP07(2012)136 
[hep-ph]. 

265] A. David et al. [LHC Higgs Cross Section Working Group], “LHC HXSWG interim recommenda- 

tions to explore the coupling structure of a Higgs-like particle,” [hep-ph]. 


266] M. Gonzalez-Alonso, A. Greljo, G. Isidori and D. Marzocca, “Pseudo-observables in Higgs decays,” 


Eur. Phys. J. C 75 (2015), 128 doi:10.1140/epjc/s10052-015-3345-5 |arXiv:1412.6038) [hep-ph]]. 


267| A. Greljo, G. Isidori, J. M. Lindert and D. Marzocca, “Pseudo-observables in electroweak 
Higgs production,” Eur. Phys. J. C 76 (2016) no.3, 158 doi:10.1140/epjc/s10052-016-4000-5 
arXiv:1512.06135 [hep-ph]]. 

[268] G. F. Giudice, C. Grojean, A. Pomarol and R. Rattazzi, “The Strongly-Interacting Light Higgs,” 
JHEP 06 (2007), 045 doi:10.1088/1126-6708/2007/06/045 |arXiv:hep-ph/0703164) [hep-ph]. 

[269] S. Willenbrock and C. Zhang, “Effective Field Theory Beyond the Standard Model,” Ann. Rev. 
Nucl. Part. Sci. 64 (2014), 83-100 doi:10.1146 /annurev-nucl-102313-025623 |arXiv:1401.0470) [hep- 
phi). 

[270] I. Brivio and M. Trott, “The Standard Model as an Effective Field Theory,” Phys. Rept. 793 
(2019), 1-98 doi:10.1016/j.physrep.2018.11.002 |arXiv:1706.08945 [hep-ph]]. 

[271] B. Grzadkowski, M. Iskrzynski, M. Misiak and J. Rosiek, “Dimension-Six Terms in the Standard 
Model Lagrangian,” JHEP 10 (2010), 085 doi:10.1007/JHEP10(2010)085 |arXiv:1008.4884) [hep- 
phi). 

[272] M. Klute, R. Lafaye, T. Plehn, M. Rauch and D. Zerwas, “Measuring Higgs Couplings from LHC 
Data,” Phys. Rev. Lett. 109 (2012), 101801 doi:10.1103/PhysRevLett.109.101801 |arXiv:1205.2699 
[hep-ph]. 

273| D. Carmi, A. Falkowski, E. Kuflik and T. Volansky, “Interpreting the 125 GeV Higgs,” Nuovo 

Cim. C 035 (2012) no.06, 315-322 doi:10.1393/nec/i2012-11386-2 |arXiv:1206.4201 [hep-ph]]. 

274| I. Low, J. Lykken and G. Shaughnessy, “Have We Observed the Higgs (Imposter)?,” Phys. Rev. 

D 86 (2012), 093012 doi:10.1103/PhysRevD.86.093012 |arXiv:1207.1093 [hep-ph]. 

275| J. Ellis and T. You, “Global Analysis of the Higgs Candidate with Mass ~ 125 GeV,” JHEP 09 

(2012), 123 doi:10.1007/JHEP09(2012)123 |arXiv:1207.1693) [hep-ph]]. 

276| J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott, “First Glimpses at Higgs’ face,” JHEP 

12 (2012), 045 doi:10.1007/JHEP12(2012)045 jarXiv:1207.1717 [hep-ph]. 

277| D. Carmi, A. Falkowski, E. Kuflik, T. Volansky and J. Zupan, “Higgs After the Discovery: A 

Status Report,” JHEP 10 (2012), 196 doi:10.1007/JHEP10(2012)196 |arXiv:1207.1718 [hep-ph]]. 


278| S. Banerjee, S. Mukhopadhyay and B. Mukhopadhyaya, “New Higgs interactions and recent 
data from the LHC and the Tevatron,” JHEP 10 (2012), 062 doi:10.1007/JHEP10(2012)062 


arXiv:1207.3588 [hep-ph]]. 


[279] F. Bonnet, T. Ota, M. Rauch and W. Winter, “Interpretation of precision tests in the Higgs 
sector in terms of physics beyond the Standard Model,” Phys. Rev. D 86 (2012), 093014 


doi:10.1103/PhysRevD.86.093014 {arXiv:1207.4599 {hep-ph]]. 


[280] T. Plehn and M. Rauch, “Higgs Couplings after the Discovery,” EPL 100 (2012) no.1, 11002 


doi:10.1209/0295-5075/100/11002 |arXiv:1207.6108 [hep-ph]. 


99 


[281] A. Djouadi, “Precision Higgs coupling measurements at the LHC through ratios of pro- 
duction cross sections,” Eur. Phys. J. C 73 (2013), 2498 doi:10.1140/epjc/s10052-013-2498-3 


arXiv: 1208.3436 [hep-ph]]. 


[282] B. A. Dobrescu and J. D. Lykken, “Coupling spans of the Higgs-like boson,” JHEP 02 (2013), 
073 doi:10.1007/JHEP02(2013)073 jarXiv:1210.3342) [hep-ph]]. 

[283] G. Cacciapaglia, A. Deandrea, G. Drieu La Rochelle and J. B. Flament, “Higgs couplings beyond 
the Standard Model,” JHEP 03 (2013), 029 doi:10.1007/JHEP03(2013)029 |arXiv:1210.8120 [hep- 
phi). 

284| G. Belanger, B. Dumont, U. Ellwanger, J. F. Gunion and S. Kraml, “Higgs Couplings at the End 

of 2012,” JHEP 02 (2013), 053 doi:10.1007/JHEP02(2013)053 |arXiv:1212.5244 [hep-ph]. 

285| G. Moreau, “Constraining extra-fermion(s) from the Higgs boson data,” Phys. Rev. D 87 (2013) 

no.1, 015027 doi:10.1103/PhysRevD.87.015027 |arXiv:1210.3977| [hep-ph]J]. 

286| T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile and M. C. Gonzalez-Garcia, “Constraining 
anomalous Higgs interactions,” Phys. Rev. D 86 (2012), 075013 doi:10.1103/PhysRevD.86.075013 
arXiv:1207.1344 [hep-ph]]. 

[287] T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile and M. C. Gonzalez-Garcia, “Robust De- 
termination of the Higgs Couplings: Power to the Data,” Phys. Rev. D 87 (2013), 015022 
doi:10.1103/PhysRevD.87.015022 |arXiv:1211.4580| [hep-ph]]. 

[288] E. Massó and V. Sanz, “Limits on anomalous couplings of the Higgs boson to elec- 
troweak gauge bosons from LEP and the LHC,” Phys. Rev. D 87 (2013) no.3, 033001 
doi:10.1103/PhysRevD.87.033001 |arXiv:1211.1320) [hep-ph]]. 


[289] A. Juste, in Proceedings of HCP2012, 15 November 2012, Kyoto, Japan, 


[290] Talk k by K K. Herner, in Proceedings of ICHEP 2014, Spain. 


[291] G. Aad et al. [ATLAS and CMS Collaborations], “Measurements of the Higgs boson production 
and decay rates and constraints on its couplings from a combined ATLAS and CMS analysis of the 
LHC pp collision data at ys = 7 and 8 TeV,” JHEP 1608, 045 (2016), |arXiv:1606.02266 [hep-ex]]. 

292| R. Kumar Barman, G. Belanger and R. M. Godbole, “Status of low mass LSP in SUSY,” 

arXiv:2010.11674 {hep-ph]]. 

293] A. Sirlin, “Radiative Corrections in the SU(2)-L x U(1) Theory: A Simple Renormalization 

Framework,” Phys. Rev. D 22 (1980), 971-981 doi:10.1103/PhysRevD.22.971. 

294] K. G. Chetyrkin, J. H. Kuhn and M. Steinhauser, “RunDec: A Mathematica package for running 
and decoupling of the strong coupling and quark masses,” Comput. Phys. Commun. 133 (2000), 
43-65 doi:10.1016/S0010-4655(00)00155-7 |arXiv:hep-ph/0004189 [hep-ph]]. 

[295] F. Herren and M. Steinhauser, “Version 3 of RunDec and CRunDec,” Comput. Phys. Commun. 
224 (2018), 333-345 doi:10.1016/j.cpe.2017.11.014 |arXiv:1703.03751 [hep-ph]]. 

[296] C. W. Bauer, Z. Ligeti, M. Luke, A. V. Manohar and M. Trott, “Global analysis of inclusive B 
decays,” Phys. Rev. D 70 (2004), 094017 doi:10.1103/PhysRevD.70.094017 |arXiv:hep-ph/0408002 
[hep-ph]. 

[297] J. A. M. Vermaseren, S. A. Larin and T. van Ritbergen, “The four loop quark mass anomalous 
dimension and the invariant quark mass,” Phys. Lett. B 405 (1997), 327-333 doi:10.1016/S0370- 
2693(97)00660-6 |arXiv:hep-ph/9703284| [hep-ph]. 

[298] B. Grzadkowski, O. M. Ogreid and P. Osland, “Measuring CP violation in Two-Higgs-Doublet 
models in light of the LHC Higgs data,” JHEP 11 (2014), 084 doi:10.1007/JHEP11(2014)084 


arXiv:1409.7265 [hep-ph]]. 


100 


[299] B. Grzadkowski, H. E. Haber, O. M. Ogreid and P. Osland, “Heavy Higgs boson decays 
in the alignment limit of the 2HDM,” JHEP 12 (2018), 056 doi:10.1007/JHEP12(2018)056 


arXiv:1808.01472 [hep-ph]]. 
[300] A. Arhrib, M. Capdequi Peyranere, W. Hollik and S. Penaranda, “Higgs decays in the two Higgs 
doublet model: Large quantum effects in the decoupling regime,” Phys. Lett. B 579 (2004), 361-370 


doi:10.1016/j.physletb.2003.10.006 |arXiv:hep-ph/0307391 [hep-ph]. 


[301] G. Bhattacharyya, D. Das, P. B. Pal and M. N. Rebelo, “Scalar sector properties of two-Higgs- 
doublet models with a global U(1) symmetry,” JHEP 10 (2013), 081 doi:10.1007/JHEP10(2013)081 


arXiv:1308.4297 [hep-ph]]. 


[302] J. McDonald, “Gauge singlet scalars as cold dark matter,” Phys. Rev. D 50 (1994), 3637-3649 


doi:10.1103/PhysRevD.50.3637 |arXiv:hep-ph/0702143  [hep-ph]]. 

[303] C. P. Burgess, M. Pospelov and T. ter Veldhuis, “The Minimal model of nonbaryonic dark mat- 
ter: A Singlet scalar,” Nucl. Phys. B 619 (2001), 709-728 doi:10.1016/S0550-3213(01)00513-2 
arXiv:hep-ph/0011335) [hep-ph]. 

[304] D. O’Connell, M. J. Ramsey-Musolf and M. B. Wise, “Minimal Extension of the Standard Model 
Scalar Sector,” Phys. Rev. D 75 (2007), 037701 doi:10.1103/PhysRevD.75.037701 
ph/0611014 [hep-ph]. 

[305] O. Bahat-Treidel, Y. Grossman and Y. Rozen, “Hiding the Higgs at the LHC,” JHEP 05 (2007), 
022 doi:10.1088/1126-6708/2007/05/022 |arXiv:hep-ph/0611162) [hep-ph]. 

[306] V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf and G. Shaughnessy, “LHC Phe- 
nomenology of an Extended Standard Model with a Real Scalar Singlet,” Phys. Rev. D 77 (2008), 


035005 doi:10.1103/PhysRevD.77.035005 [hep-ph]. 
307| X. G. He, T. Li, X. Q. Li and H. C. Tsai, “Scalar dark matter effects in Higgs and top quark 
decays,” Mod. Phys. Lett. A 22 (2007), 2121-2129 doi:10.1142/S0217732307025376 
[hep-ph]. 
308| H. Davoudiasl, R. Kitano, T. Li and H. Murayama, “The New minimal standard model,” Phys. 


Lett. B 609 (2005), 117-123 doi:10.1016/j.physletb.2005.01.026 \arXiv:hep-ph/0405097 [hep-ph]]. 


309| S. Kanemura, M. Kikuchi, K. Sakurai and K. Yagyu, “H-COUP: a program for one-loop corrected 
Higgs boson couplings in non-minimal Higgs sectors,” Comput. Phys. Commun. 233 (2018), 134- 


144 doi:10.1016/j.cpc.2018.06.012 |arXiv:1710.04603 [hep-ph]]. 
[310] S. Kanemura, M. Kikuchi, K. Mawatari, K. Sakurai and K. Yagyu, “H-COUP Version 2: a 
program for one-loop corrected Higgs boson decays in non-minimal Higgs sectors,” Comput. Phys. 


Commun. 257 (2020), 107512 doi:10.1016/j.cpc.2020.107512 |arXiv:1910.12769 [hep-ph]. 

[311] W. Frisch, H. Eberl and H. Hlucha, “HFOLD - a program package for calculating two-body 
MSSM Higgs decays at full one-loop level,” Comput. Phys. Commun. 182 (2011), 2219-2226 
doi:10.1016/j.cpe.2011.05.007 |arXiv:1012.5025 [hep-ph]. 

[312] H. Bahl, T. Hahn, S. Heinemeyer, W. Hollik, S. Pa®ehr, H. Rzehak and G. Weiglein, “Precision 
calculations in the MSSM Higgs-boson sector with FeynHiggs 2.14,” Comput. Phys. Commun. 249 
(2020), 107099 doi:10.1016/j.cpce.2019.107099 |arXiv:1811.09073 [hep-ph]. 

[313] J. Baglio, T. N. Dao and M. Miihlleitner, “One-Loop Corrections to the Two-Body Decays 
of the Neutral Higgs Bosons in the Complex NMSSM,” Eur. Phys. J. C 80 (2020) no.10, 960 


doi:10.1140/epjc/s10052-020-08520-8 jarXiv:1907.12060) [hep-ph]. 


101 


